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SYNOPSIS 


This study is concerned with the problem o£ estimating 
the parameters and testing the equality of location parameters 
©i Ci = 1#2, of K(> 2) two-parameter e3q)onential distribu- 

tions E(0^^a) based on K independent type II censored samples* 
■Heare the common scale parameter a is assumed to be untaiown. 

Type II censored sample is an ordered sample in which a Tcnown 
number of smallest (left) observations anchor largest (right) 
observations are missing* An oirdered sample is obtained by 
rearranging the variates in an ascending order of magnitude# 

Most of the work in this field is done when complete or 
right censored samples are available# However, there are 

ii 

situations when some smallest observations are also not availably 
In the present work, main attention has been paid to type II [ 
doubly censored samples# j 

4^12 4^3^ (i = 1,2,...,K) be K I 

independent ordered samples from E(0^,a) with r^ > 0, Sj^ > 0 
and rj+1 < n.-— s . • Based on these observations, the Least Square | 
(LS) and Maximum Likelihood (ML) estimators of 0^ (i =5 1,2, 
and o are obtained# For K = 2, the LS and ML estimators are 
derived under the assumption 0^ “ ®2 “ Some distribution 
theory results regarding these estimators are also obtained* 

A brief comparison of the LS and ML estimators is made by using 
the mean square error criterion# 



SYNOPSIS 


This study is concerned with the problem of estimating 
the parameters and testing the equal ity of location parameters 
0 ( i = 1# 2/ • • -/K) of K(> 2) two-parameter exponential distribu- 

tions E(0^, a) based on K independent type II censored samples* 
Her^ the common scale parameter (7 is assumed to be unknown* 

Type II censored sample is an ordered sample in which a known 
number of smallest (left) observations anchor largest (right) 
observations are missing* An ordered sample is obtained by 
rearranging the variates in an ascending order of magnitude# 


Most of the work in this field is done wiien complete or 
right censored samples are available# However, there are 
situations when some smallest observations are also not available# 
In the present work, main attention has been paid to type II 
doubly censored samples# 

Let ***'^n^s (i = 1/2, ...,K) be K 

independent ordered samples from Eid^go) with r^ > 0, > O 

and Based on these observations, the Least Square 

(LS) and Maximum Likelihood (ML) estimators of 6^ (i = 1, 2, 
and a are obtained# For K = 2, the LS and ML estimators are 

Some distribution 
theory results regarding these estimators are also obtained# 

A brief comparison of the LS and ML estimators is made by using 
the mean square error criterion# 


derived under the assumption 0^ “ ®2 “ 



For K = 2/ the test statistics based on LS and ML 
estimators are proposed for testing the null hypothesis 6^ = 
against one-sided alternatives* These statistics are equivalent 
to 

^ - ^\+l > 

where 


= F^/d is the pooled estimator of Oj 


K ^i~®i 
P- = S { S 

i=l j=rj^+l 

K 


^(i) , ^ ^(i) , . ^^(i) , 


d = S (nj_~r£-.SjL-l ) and K = 2( 
i=l 


Against the alternative G^ i> © 2 / large values of T lead to the 
rejection of the null hypothesis* The null and non-^null 
distributions of T are derived# Some exact and approximate 


critical points of T are tabulated# Power values along with 
a normal approximation are also tabulated# From this^ it is 
concluded that the test T is more sensitive for the left 
censoring than for the right censoring# 


The test statistics ^^^^2 and K based on LS estimators/ 
ML estimators and LR test procedures respectively/ are proposed 
for testing 0^ = ©2 against 0^^ © 2 - These are given by 

= I T-q^l 


(1 = 1 , 2 ), 



and for Y > 0, 


X = 


Y ^ ^ (d* a^+fY-.P)^ exp(-p/a ) 


° { d"'^ a^+( r^+f ) Y~P} { P-d* + r^Y} 


where 52 ~ 

= ■ E -(n^-^j+l) ^ (i = ly2), 

j=l 


= log {n^/(nj^~r^) } (i = 1,2)^ 
a = P^/d^ ^ d* = d+2^ 
f = n3_+n2 - r^-^r2, 

2 n. -(n^^r. ) 

R = £ n ri^^(n^~rj_) ^ } exp (d*)/ Y = 

i==1 


2 


- = E £ S ' -(ni-ri) 

i=l j=ri+l J "" ''i ®i "" "■ 

and 0Q is the solution of the equation 


V^o 

e = {1 + 


^2^ 


P - d*a 




r^Y 


d* cTq + fY-P 


The expression for X given above is for r^^ > 0/ r 
Considerably siinpler e3?pressions are obtained if r 
r 2 s= O* For Y < O, X is obtained by replacing n^,: 
’^2'^!' ^2'' ^ above expression# 


.(2) ^^(1) 
■r2+l r^+1 


> 0 -- 

= O an^or 
,/r^,r 2 by 



The exact critical points and power values of the testi 
and ^2 are evaluated. Due to the complex nature of X , only 
simulated ciritical points and power values of X are tabulated# 

On the basis of power calculations, it is concluded that the test 
V2 is somewhat more biased than and there is very little 
difference between and X# Since the statistic X is far more 
complicated than use of is recommended* 


For testing K(> 3) populations, two test statistics are 
proposed* We first consider the case of right censoring, so that 
the smallest observation of each sample is available# For 
simplicity of notations let =s (i = 1 , 2 , #.#,k) and 

X^l) = min(X^,X2, • * • The statistic given by 

= {X^~mir(X2,Xg, ..-,Xj^)}/a* 


is proposed for testing : 
^ max(02/ 


=©2 = = 0 
and the statistic 


against 


! 


T^ = { max (X. ) - min (x. )}/a* 
l<i<K l<i<K 

is proposed for testing against H2 : at least one 0 ^ is 
different from 0 • The exact critical points and power values 
of the tests T^ and T2 are evaluated# For K = 3 case, the 
performance of statistic T^ is studied for different combinations 
of n^,n2#n3 and d# It is observed that, the test T^ is more 
sensitive for changes in n^^ than changes in ■a2t'n.-^ and d# For 
K = 3 and equal sani>le size case, the performance of T^ and T2 



along with that of and is studied* The statistics and 
are given by 


and 


2 


K 

U4 = n^(Xj^-X(^))/{ (K~l)a*} . 


The test is actually the LR test and has been proposed by 
other v^ithors* On the basis of power calculations carried out# 
is recommended for testing against and the LR test 
statistic is recommended for testing against ^ 2 * 

If the smallest observation is missing in each sample of 
size n# then the following test statistics are studied ; 


V, = {X^^^ ~ min (X^^^)}/a* 

2<i<K ^ 

for testing against 

V, = £ max (X^^^) - min (X^^^)}/ff* # 
l<i<K l<i<K 

and ” 

V- = [ max £ (X^^^-X^^ ^)# ( X^^ ) } ]/a* 

2 <i<K z z 

for testing against H2* For K = 3 # exact critical points 
and simulated power values of these tests are tabulatejd* On the 
basis of these calculations is recommended for testing 
against but nothing can be said regarding the pireference of 
V2 over as in some cases performs better and in some cases 
V3 perform^ better * ' 



CHAPTER I 


INTRODUCTION AND SUMMARY 


1*1 Scope * 

Exponential distribution is often proposed for modelling 
the lifetime distributions of items like electronic components^ 
mechanical breakdowns^ light bulbs etc* [ Davis (1952) J Epstein 
(1958)/ Proschan(l963)/ Nelson (l975)J* In a two-parameter 
exponential distritautionj the location parameter is interpreted 
as the minimum (or guarantee) time/ before which no failures 
occur/ and the scale parameter/ as the mean life measured from 
the location parameter as the starting point* 

There are several situations/ where the complete 
sample is neither available nor desirable* Since life-testing 
experiments are usually destructive/ this limits the number of 
items to be tested (Sinha and Kale 1980, p» 18)* Moreover/ in 
the ordered samples frequently found in biological data either 
some smallest and/or some largest observations are not available 
[ Ipsen (1949)]* An ordered sample is obtained by rearranging 
the variates in an ascending order of magnitude* In an ordered 
sample# if a known number of smaillest (left) values or largest 
( right) values or both are missing# then such a sample constitutes 
a type II censored sample* 

Several problems dealihg with the estimation and testing 
of hypothesis of a two-parameter esqjonential distribution are 



discussed in several books, for example see Sarhan and 
Greenberg (1962), Mann, Schafer and Singpurwalla (1974), Bain 
(1978), Sinha and Kale (1980)* Complete and right censored 
samples have been considered by many authors (Walsh (195-0), 
Halperin (1952), Epstein and Sobel (1953), Hogg and Tanis (1963), 
Grubbs (1971), Kumar and Patel (1971), Dubey (1973), Weinman et 
(1973), Khatri (1974), Pemg (1978), Mathai (1979), Regal (1960), 
Bhattacharyya and Mehrotra (l98l), Hsieh (l98l), Goria (1982), 
Mehrotra and Bhattacharya (1982), Singh (1983), Singh and 
Narayan (198 3) etc*) • 

Although the left censored samples have not been 
considered so thoroughly as the right censored samples, there 
are situations in v^ich some smallest observations are not 
available* For exaitple, in experimental biology, n animals 
are tested for antibodies after a certain period of time* Only 
(n— r) of these samples contain measurable amounts while r of 
the animals develop the antigen at a level too low for measurem^t 
by the prevailing technique (Ipsen, 1949) * This gives rise 
to a censored sample from left* Another example where the 
smallest order statistic is difficult to observe is the failure 
time of human kidney It is not easy to tell the failure time 
of one kidney since noticeable symptoms occur only when both, 
kidneys fail* Similarly* in a transistor set four battery 
cells may be used* Failure of only one cell may not affect 
the performance, and hence the first failxxre time may go 
unreportedi; But ^en two cells fail, the transistor may not 



3 


work properly and the second smallest order statistic becomes 
the first available observation* 

Based on doubly censored sart^les, the Least Square (LS) 
estimators and Maximum Likelihood (ML) estimators of the 
parameters of a two~parameter exponential distribution were 
derived by Sarhan (1955) and Kanibo (1978) respectively* Recently, 
Tiku (l98l) and Khatri (1981) have considered the probl^i of 
testing equality of location parameters of two-parameter 
exponential distributions based on type II doubly censored 
samples* 

One of the problems arising in life-testing experiments 
is the estimation and the comparison of minimum (or guarantee) 
time of the items* Related to this probl€3n, the present work 
is conceimed with the problem of estimating the parameters and 
testing the equality of location parameters of K(> 2) two- 
parameter exponential distributions* Here the scale parameters 
are assumed to be equal taut unknown* In this thesis, main 
attention has been paid to type II doubly censored samples* 

The following topics are studied in this thesis : 

1* Comparative study of the ML and LS estimators of the 
parameters * 

2* Testing equality of location parameters against one-sided 
alternatives * 

a* Testing equality of location parameters against two-sided 
alternative* 



4 * 


4 


Generalized statistics for K right censored samples* 

5* Generalized statistics in equal saitple case, when one 
observation is missing on the left* 

The results obtained are compared with the existing 
results* Suitable graphs and tables aire provided to support 
the theory, wherever necessary* The notations v/hich are used 
consistently in the text arc given in the next section and the 
subsequent sections describe briefly the above mentioned topics* 

1*2* Notations and abbreviations * 

As far as possible random variables will be designated 

by upper case letters, and their realizations (observations) by 

! 

the corresponding lower case letters* 

1th ordered sample of size 
nj^ (i a 1/2, • **,K) with superscript 
i dropped If there is only one sample* 

ith type II censored sample in which 
first rj^ and last Sj^ observations 
are missing, where r^ ^ O# ®i - 
and rjL+1 < Hj^-Sj^* 

cumulative distribution function of X 
probability density function of X 
mean of X 
variance of X 
mean square error of T 
standard error of T 
likelihood rs^io 


i) A 
"i ®i 

Fj^(x),P(x) a P [x < x] 
. f^Cx)/ p(x) 

E(X) 

Var(x) 

MSECt) 

sect) 

LR 
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4 * Generalized statistics for K right censored samples* 

5* Generalized statistics in equal sanple case^ when one 
observation is missing on the left* 

The results obtained are compared with the existing 
results* Suitable graphs and tables are provided to support 
the theory, wherever necessary* The notations vhich are used 
consistently in the text arc given in the next section and the 
subsequent sections describe briefly the above mentioned topics* 

1*2* Notations and abbreviations * 

As far as possible random variables will be designated 
t'Y npper case letters, and their realizations (observations) by 
the corresponding lower case letters* 

ith ordered sample of size 
n^ (i ss 1/2/#**,K) with superscript 
i dropped if there is only one saitple* 

i th type II censored sample in which 
first r^ and last Sj^^ obseirvations 
are missing, where r^ >0, > 0, 

and r^+1 < 

cumulative distribution function of X 
probability density function of X 
mean of X 
Variance of X 
mean square error of T 
standard error of T 
likelihood r»tio 


<•'•< 

„( i) •) 

ni-Si 

Px(x),P(x) = P [x < x] 
• f^Cx), p(x) 

E(X) 

Var(x) 

mse(t) 

sect) 

LR 



KP test 


5 


LS 

pdf 

jp<af 

cdf 


Kumar and Patel (1971) test 
least square 
maxirnum likelJLhood 

independent identically distributed 
probability density function 
joint pdf 

cumulative distribution function 


w»r«to 


with respect to 


MVU 

DP 


Me, a) 


h ) 


V 


a*,8* 


^ ^ A A 

cj /© /O .e 
o' o' O'' o 



minimum variance unbiased 

degrees of freedom 
two~parameter exponential 
distribution with location 
parameter 0 and scale parameter a 

“normal distribution with mean ju 
and variance 

“central chi-square distribution 
with V DP 

“LS estimators of a and 0 
respectively 

LS and ML estimators under the 
hypothesis 0^ - ^2 = 

is distributed according to « 

is approximately equal to 


- is identically equal to 

B(a,b) = X t^~^(l-t)’^^ dt,a > 0, b > 0 
o 

TTa) = X t^”^ dt, a > 0 
o 
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QpCxI 


s) = 


^ ^ dt/(p-l)i , X > 0, p = 1^2, for s > -1 


X 


p-1 


Lp(xis) 


S e3qp £-x(l+s)}{x(l+s)} Vui (l+s)^} 
j=o 

X 

= / ^t/(p-l)! , X > O, p * 1,2,... 


xP/pI , s = 1 


l/(l~s)^ ~ S eXp{-J-x(l-s))x^(l-s)^'’^/j| , s 1 
j=0 
K 

E (n^-ri-Si) 

K 

E (n^-r^-s.-l) = d*~ K 
i=l ^ 

E £ E ^ X^^| } 

i=l j=r^+l J ^ ^ ^ 


_ n.-s, 

2 i i 


(i) 


i=l j=rj^+l J 1 nj^ 


- (“i-^i) 4 J!i5 


1 *3* Comparative study of the ML and LS estimators of the 
parameters » 

Let ^2-4-2' ‘'^n~s ^ censored sample 

from an exponential distribution E(0,a) with pdf 


f(x?6gCj) = a exp {— (x-^)/ct}, O < a < «>^ G < x < <»# 

Sukhatme (19 36) obtained the best unbiased estimators of 0 and 
a based on a conplete sanple (r = 0, s s= 0 case) of size n. 
Lloyd (1952) discussed the technique of estimating 0 and a by 
applying general LS theory to an ordered sample. By applying 


} 



/ 


Lloyd's method Sarhan (1954/1955) and Greenberg and Sarhan (1962) 
have obtained the best linear unbiased estimators of mean and 
standard deviation for dcxible and middle censoring case as well* 

Let ^r^il'^r^+a' “ 1/2/..-/K) be K 

independent type II doubly censored samples from E(0^/Cf)» The 
LS estimators 0^ arid d* of 8j[(l- * 1/2/*»*/K) and o are obtained 
in Section 2*2# These are given by 


® i = - V* < i = 

^i"^^ -1 

where b. = I (n^-j+l) 

0=1 


1/2/. ../K) and 0 * = P^/d, 

K 

(i s= 1/2/-../K)/ d =: S ( n .-r .-s .~l) 

1=1 ^ ^ ^ 


and P. = E £ E + s^ 

^ i=l j=ri+l J ”i ®i ^ ^ 


In Section 2*3/ it is shown that (i = 1/2/.../K) and a* 

^2 

are independent, and 2d0*/a has ^ ^2d distribution * 


ML estimators of the parameters for a complete and 
right censored sample were derived by Sukhatme (1936) and 
Epstein and Sobel (1954) respectively. For a doubly censored 
sample, ML estimators were discussed by Tiku (1967). He has 
obtained the ML equations and has mentioned that, those equations 
do not have explicit solutions. Hence, he obtained a modified 
ML estimator of the parameters. Using these equations, 

Kambo (l978) has obtained esqjlicit expressions for the ML 

A /V 

estimators* In Section 2*4/ the ML estimators 0^ and a of 
0j^ (i = 1/2/*../K) and O are obtained. These are given by 



8 


* ( t) A 

©i = + a log (l-r^/nj,) (i 


where d = d+K« 


= 1,2, 0 *^,k) and a * 




Kambo (1978) conpared the minimum variance unbiased 
(MVU) estimators with the ML estimators, 'hJien a doubly censored 
sample is available* He has shown that for a single sample 
MSE(a*) > MSE(cr) and verified mmierically that MSE(6*) can be 
greater or less than MSE(e)* By this, he concluded that, some 
times ML estimators are better than MW estimators* Similar 
type of comparison is done in Section 2 *5 and it is concluded 
that in general for K > 3, MSECa^) is less than MSE(a)» However, 
for K < 2, MSECa"^) is greater than MSE(a), 


Epstein. and Tsao (1953) derived the ML estimators 
under the hypothesis 0^^ ~ ^2 ~ ® right censored samples* 

In Section 2*6, the LS estimators and ML estimators for doubly 
censored samples are discussed under the same hypothesis* A 
brief comparison has been done between LS estimators and ML 
estimators in this case also* 


1 *4 * Testing ecaialitv of location parameters against one">sided 
altematlves • 

Eor the right censored samples from two populations 
Kumar and Patel (1971) have proposed a test statistic for testing 
H^ : 0^ s= ©2 against H 2 : 0^ ^ © 2 * Weinman ^ ^* (1973) 
extended it for testing H^ against one-sided alternative 
H' : < © 2 * Their statistic is W = where 



^ ~ i® tine pooled estimatoi’ of Cf • They obtained, tbe 

critical point c as 

d[{nj_/(n3_+n2)a}^/^-l]/n2 if < an2/(l-a) 

G ss 

_ d [ l-{ n2/ ( nj^ +n2 ) ( 1-a ) } otherwise, 

where a is the chosen level of significance. Note that c > O 
for n^ > an^ iX-Qc) • They obtained the power function P(<P) of 
W for c > 0 and <P = (e^-ej^Va > 0 as 

P(^p) = l-{ V ( n. tn- ) } [(n.+n«) 2 {6.<p/c)^/it 

^ 1=0 

•'{ n^/d+rioc/d) i S {q)(d+n-,c)/c} I ] 

i=0 ^ 

-[rue ^ E {<P ( d-n- c)/g} yi I ] 

^ i=0 

. l/{ ( nj^+-n 2 ) ( l-Hji c/d)^} , 

provided that ^ d* For the case n^c = d, the last term in 
P ( ^ ) becomes 

"■n-tp d ' 

— n2e (d 9 /c) /(n^+n2)dl» 

For c < 0^ P(<P ) i.s given by 

P(<P) = L-n2 e3qp(-n^g)) Cl— n^c/d) /(n^tn2)* 

In Chapter III, a test statistic based on ML and LS 
estimators is proposed for testing against : 0^ > 02 * 
This test is equl-valeht to 
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The distributions of the statistic T under and are 
derived. Approximations for null distribution in terms of 
Student^ s t and normal distributions are studied. Appxoxlinate 
critical points obtained from above approximate null distributions 
are also tabulated along with the exact critical points. It is 
observed that the normal approximation is better if r^ > r 2 t 
otherwise Student's t approximation is better. Some exact and 
normally approximated power values are tabulated. Also/ the 
variation in power due to different combination of r^^ and r 2 
is plotted. With power function as a base^ it is concluded 
that the test T is unbiased and it is more sensitive for the 
variations in r^^ compared to the variations in r 2 . Further^ 
the performance of the test is not seriously affected for 
variations in right truncations for fixed values of 

^ 2 - 

1 *5 . Testing equality of location parameters against twD~sided 
alternative » 

Epstein and Tsao (1953) discussed the LR test 
procedure for testing various types of hypotheses based on right 
censored samples. Kumar and Patel (l97l) proposed a test based 
on 1 ^ )/a^ 1 for testing : 6^^ = ©2 against H 2 : ©^/402* 

Dubey (1973) and Weinman et al « (1973) derived the power function 
of KP test. Weinman ^ coitpared the KP test with LR test 

and they concluded that in general LR test performs better than 
KP test. The power function P(<SP) of KP test with critical 
point c was obtained by Weinman et. al. yor ^ * (02~6j_)/c ^ 0 



it is given by 


d-1 

P((p) = l-exp(-d<p/c) E (d(p/c)/H 

i=l 

+ { n 2 /(nj^+n 2 )} exp(~nj^<p){ (l+n^c/d) ( 1-nj^c/d) ^ 

( 3.***1 

+ { n^/(n^+n2 )} (l+n2c/d)"’'^ S {(p(n2+<^c)} ^/±l 

IssO 

+ {n2/(n^+n2)} Cl-nji^c/d)"'*^ 

for 1-n^c/d 0* However/ for 1-nj^c/d = O the third and last 

terms of P(<p) become respectively 

P."* { n2/(nj^+n2)} e3^(-n^(p)and ~n2(<i<P/c) e3^(-nj^<p)/{d I (n^+nj)}* 

The power P(<P) for (p < 0 is obtained by interchanging n 2 and n^ 
and evaluating the above expression for 1 <pl • 

Recently/ Tilcu (1981 ) considered the problem of 
testing : 8^ =02 against ^2 ' ^ ®2 ’^^^sed on doubly 

censored samples* He prc^osed the test statistic t 

^ 

for testing H^ against H 2 and obtained its null distribution as 
^2 . r 

f^(u) = H {] E (~l)^( j^) B( f+j/r^^+l ){ l+h 2 ( j )u} ^ 

j=0 

^1 . r 

+ E (- 1 )^( >)B(f+j,,r 2 +l){l+h^(j>u>’'‘^“’^] / u > 0/ 
j=£) ^ 

2 .1 

where H = 7~l {B(nj— rj/r.+l}, / ^i^J^ ~ (nj^— r^^+j )/d (i = 1/2.)/ 



Khatri (l98l) has derived the non-null distribution of U as 
follows : 

1 j r • d j • ^ - 

g^(u) = H E ?)B(f+j,r 5 +l) E (d(p/u) { 1+h. ( j )u} 

j =0 J -1=0, ^ 

^2 . r 

. e3{p(-d<})/u)/il + H E .^)B(f+j,ri+l){l+h«(j)u}~^"’^ 

j=0 J • ^ 

* exp {-h2(j)d(f>} - H ^ (-l)^( j^)B(f+j,rj^+l){ l-h2( J)u}”^'’^ 
d . ■ . 

. exp{*-d<P/u} [ E (d(j)/u)^{l-h^(j)u} yin , u > O, 
i=l ^ 

where <p = (0^-02)/<7 >0* This final expression is slightly 

wrong due to some integration errors • Puirther^ he has not .taken 

into account the singularities at u =* l/h 2 (j) for j = 0,1/ 2/ r 2 

(for <P > O), as has been done by Weinman ^ (197 3) for 

right censored samples* We give correct form of this expression 

in Section 4 *2 • 

In Chapter IV, two statistics defined by 
= |T-qj_i and V 2 = I T-q 2 1 based on lS and ML estimators 
respectively are proposed for testing against H 2 , where 

^ == ” ^r2+2^/^*' % = ^1~^2' ^2 * ^-^ 2 * 

^i+1 ' 

b£ =r E (nj_~j+l)~^ (i = 1 , 2 ) and m^ = log{nj/(n^-rj_)} ( i=l, 2) • 
a — 1 

The null and the non-null distributions of those statistics aare 
derived* As a special ease, the noi>-null distribution of Tiku's 
st.atistic(XJ) is also obtained* 



The LR test statistic X is derived for testing 
against For r^ > 0 , r2 > O, Y > 0 , it is given by 


. A ^ 

= const. [7-]'^ 


r. +r„ ^ . 

1 2/ ,* 


(d* a^+fY-P)-^ exp(-P/0Q) 
{ d*aQ+(rj^+f)Y~P} ^"^^{P-d* 


-(n.-r. ) “(n^-r,) ^ a jl 

where const = ^2 ^ ^ ^ e3^(d*)/0 = P^/di 


2 ^i^^i 


^ i=l j=r,+l J ^ ^i-®i ^ ^i+l-* 


Y = , P = P. +(n„— r^)Y and *_ is the solution of 

JL2'T'i JTji^'T'X X ^ O 

the equation 


A 

Y/a, 


1 + 


^2^ 


•W' A 

P-d* 0 


'/ 


1 + 




d* S„ + £Y-P i 


For Y < 0 , X is obtained by replacing nj^/n2/r^/r2 by n2#nj^/r2/r^ 
respectively in the above es^ression* 


The critical points of the tests ^ 

tabulated. The corrparative performance of all the four test 

/ 

statistics is studied. Since the distribution of X is not 
easy to derive, only simulated critical points and power values 
are used. Using the power function as a base, it is concluded 
that the test U is more biased than and V2/ is more biased 
than ,and there is very little difference in the power values 
of X and . Since the statistic X is very complicated, while 





the statistic is considerably simple^ v;e strongly recommend 
the use of the test statistic in such situations* 

1 *6 * Generalized statistics for K right censored samples * 


Khatri (1974) derived the LR test and two test 

statistics and using union intersection principle for 

testing : 6^=0 2 =^ ® against : at least one 0 ^ 

is different from 0, based on K independent right censored 

samples xj (i = from E(0 .,ct), For 

1 1 

C ' "S 

simplicity of notations^ let X^ = X^^ (i = 1 ^ 2 , ••*,k) and 
^(l) ~ niin(Xj^,X 2 / * • • Then Khatri' s statistics are given 

by 


u 


2 


[ max {nj(Xj~X/- Oll/db* 
l<i<K 1 VJ--' J 


and U- = [ max £ n. (Xj~X- ), n- (X^-^X. )} ]/da** 

2<i« X i -L i 

He has obtained the null distributions of ^3* ^'^rther, 

he has discussed their non-null distributions without carrying 
out any power calculations* Singh (1983) also discussed the 
LR procedure for testing against B. 2 * He obtained the LR test 
U 4 as 

K 

U 4 = E n^(x^-x^^))/{ (K-l)a*} = dU^/(K-l) 

and has shown that, the null distribution of is ^ 2 (K-l), 2 d ' 
However, he has not studied the power function of • 



It. Is reasonable to study a test for testing against 

the alternative : 6^ > max (6j)» For this purpose, in 

^ 2 <i<K 

Chapter V a test statistic given by 
Tl = {Xj^-minCx^^Xg, 

is proposed# Although there are several tests (as mentioned 
above) for testing against H2/ yet we propose another test 
based on 

- { max (X. ) — min (X. )}/a*# 
l<i^ 

Chapter V is mainly devotcid for studying the performance of 
for K = 3 and different combinations of n^,n2^n^ and d# 

For comparison purposes the performance of statistics 
and T2 for K = 3 and n^^ = ^2 - is also studied t Note 

that for equal sample sizes = nT2/d* The necessary critical 
points of these tests are tabulated* Some exact power values 
of the tests T^ and T2 are tabulated for points in the 
parametric space satisfying 6^ > ©2 ^ 63* Siimalated power 
values of T^,T2#U3 and are also tabulated, since the 
expressions for power functions of and provided by 
Khatri ( 1974 ) arc extremely complicated* From these calculations 
of power values, it is concluded that 

( i) the test is more sensitive to changes in 
compared to arKi d, 

(ii) the power of tt^ test T^ is oonsiderably higher than 
that of other three statistics. 



(ill) the test T2 performs slightly better than if 
is small i otherwise reverse is the case^ 

( iv) the test T2 performs better than if 
largO/ otherwise reverse is the case* 

Finally, the LR test (or equivalently U^) is 
recommended for general alternative hypothesis since its 

critical points are easy to evaluate from the F— distribution. For 
the specific alternatives like the statistic is recommended, 
since its critical points are available in a compact form and 
its power is considerably higher than that of other tests* 

1 * 7 . Generalized statistics in equal samp le case, when one 
observation is missing on the left * 

In Chapter VI, tests for the equality of location 

parameters of K populations are considered, when the smallest 

observation is missing and atleast second smallest observation 

is available in each sample of equal size n* Here the test 

defined by 

V- = - min 

^ ^ 2 <i<K 

is proposed for testing i 6^ = ©2 " ®K 

max (02/63 / '• -/©k^ • Similar to the previous section, 

the statistics V2 and V3 given by 



and 



are proposed for testing against H2 s atleast one 0^^ is 
different from 0 « Compared to two— sample case (K ,= 2 )^ the 
LR test is much more complicated even for K = 3 # Consequently^ 
this has not been studied at all* 

For K = 3 , the exact critical points of all the three 
tests are tabulated for scxne selected values of n and d« Since 
it does not appear simple to evaluate the non— null distributions 
of these statistics/ the power of these tests are calculated by 
Monte-Carlo techniques for 0 ^^ > ©2 ^ 03* From this study, we 
have recommended for testing against a specified alternative 
like Hj^/and for testing against the alternative H2/ the test 
is recommcaided if Q^roected to be very small, otherwise 

V.^ is recommended * 



CHAPTER II 


ESTIMATION OF THE PARAMETERS AND BASIC DISTRIBUTION THEORY 

Introduction » 

In this chapter, derivation and the cc»nparison of 
the Least Square (LS) and Maximum Likelihood (ML) estimators 
of the location and the scale parameters of two— parameter 
exponential distributions are discussed* Some of the results 
are established, which are used in later chapters* 

Let X ,,X ..«,X be a type II censored sample. 

r+1' r+2' n~s 

from an ejcponential distribution E(6,cr) with pdf 

(2-1*1) f(x7e,a) = cr"'^ exp{-(x-e )/£J} , 0<a<«>, e<x<«>. 

Many researchers have investigatcjd the problem of 
estimation of the scale parameter a and the location parameter 
G * For a complete sample (case r = s «= o) of size n, 

Sukhatrae (19 36) obtained the best unbiased estimators of 0 
and a as 

n 

(2*1*2) 0* = x.-a*/n^ a* = ( E x.-nX^)/d, 

3 = 1 '^ 

where d = n— 2 >'S— 1 = n— 1 for r = s = O* 

Lloyd (1952) discussed the technique of estimating 0 
and a by applying genera], LS theory to an ordered sample* 

Type II censoring on ^ the right (case r = O) was considered by 



Sarhan (1954), and he obtained the LS estimators of 6 and a 


as 

( 2 . 1 . 3 ) = X^-(7*/n, a* = ( Xj+sX^_g-nXj^)/d. 

The LS estjjnators of 6 and a based on a type II doubly censored 
sample was discussed by Sarhan (1955), and is given by 


( 2 #1 . 4 ) 0* = X «-cJ* 


r+1 . n~s 


Generalization for K > 2 type II doubly cotisored seunnples from 

are discussed in Section 2.2. The necessary distribution 
theory results are obtained in Section 2.3. 


The ML estimators of 0 and a for a complete sample 
was obtained by Sukhatme (1936) as 


/ \ ^ A ^ 

(2«1,5) 0 = X^, a = ( E Xj.-nXj^)/d* , 

where d* = n-r-s = n for r =s s = O. For type II censoring on 
the right, the ML estimators was discussed by Epstein and Sobel 
(1954). This is given by 


( 2 # 1 . 6 ) 


^ A 

0 = X^, <7 


n-s 

( S X^+sX„ -nX. )/d^ 


Tiku (l967) has derived the ML estimators of the parameters 
based on a type II doubly censored sample. He has obtained the 
following ML equations : 



(2-1-7) i [l-- 
or ^ n 

and 


(2-1-8) g [-(l-f 


f ( z ) 


ns 


^ ,2) ^ £ p 


^ j=r+l J 


"*' n ^n-s~ n ITz 


_ f ( z - ) 
r r+l _ -1 

“ 7. V+lJ 


r+l 


= 0, 


where f(z) =s exp(— z) and F(z) = l“*exp(~z) are the pdf and cdf of 
Z = (X-0)/cy respectively# According to Tiku (1967), ML equations 
(2-1-7) and (2-1-8) do not have explicit solutions (for r > o) 
duG to the presence of the term f(z)/F(z)- He therefore obtained 
modified ML estimators and of 0 and a, which are 


given by 


r+l 


e„od = ^r+l-^mod .f, 

3 


n-s 


'’mod = f 

j =r-fl 

In tsquations (2-1-7) and (2-1-8), Kajnbo(l978) 
eliminated ^ ( V^( ) and obtained explicit expressions for 

the ML estimators- These are given by 


(2-1.10) 0 = ^i-+i + ^ (1-r/n), a = 

Note that, for r = 0, equation (2-1-10) reduces to equation 
( 2 -1 -6 ) which in turn reduces to equatim ( 2 *1 -5 ) for s = 0 • 

Prom equations (2-1-4) and (2-1-10), it is easy to see 

that 


(2.1.11) 0* = - d* S 


r+l 


r+l - % 

1 - 2 9 . 

' ^ - d 





In Section 2 *4, ML estimators for K(> 2) samples are 
derived* A brief study of MSE of ML estimators and MVU estimators 
is made in Section 2*5. 

Epstein and Tsao (1953) considered the problem of 
testing equality of two exponential distributions based on 
right censored sarrples* They derived the ML estimators under 
the hypothesis 0^^ = ~ These are given by 


( 2 . 1 * 12 ) 


6^ = min 






2 

where d* = E (n^-rj^-s^) with r^^ - ^2 ~ doubly censored 

samples, the LS estimators and the ML estimators xinder the 
hypothesis 0^ = 02 = 0 are discussed in Section 2«6* 

2 • 2 • LS estimators of the parameters for K(^ 2) samples case * 

In this section, the LS estimators for K samples are 
derived* Let (i = 1 , 2 , ***,K) be K 

ri+l ®1 

independent type II censored samples from E(0j^,a)* Denote 
E(x) = A r and Var(x) = a^D, 

fs> «Nt ^ 


where 


E(Xj ^) = ©i + Cj 

Var(X^^^) = 

= S l/(n.-g+l), 
^ g=l ^ 


S l/{n^-g+l)% 
g=l 



X = 

fNJ 


x(l) 

X(2) 

« 

X(K) 


, X(i) = 


„(i) 

r^+l 

„(i) 

^(i) 

"i-=i 


r = 


e 


> -© = 


e. 


e. 


e 


K 


A = 

ro 


1 

0 ^ m m 

ro 

0 

CO 

c(l) 

CO 



- -n 

'^r.+l 

0 

«v> 

# 

• 

1 • ♦ • 

ho 

0 

CO 

c(2) 

CO 


, C(i) = 

fO 

(i) 

ri-i-2 

# 

# 

0 

0 #«• 

fsJ 

1 

CO 

c(k) 

CO 

'd*X(K+l) 


"l-'i 


and rankCA) = K+i; D is the dispersion matrix of X/a/ and is 
given by 




D. 0 

ho 

c c c 

0 

CO 


D 

CO 


0 ■ D- 

• 

c 

0 9 9 

0 

CO 

i 



0 0 

CO <o 

0 0 0 

Sk 

Mii 

d*xd* 



^r^+l 

^r^tl 

*•• %+l 

Si 


^(i) 

^r^+l 

^r^+2 

,( 1 ) 

• • • %+2 



m 

9. 

c 






,(i) 

^rj^+l 

r^+2 

Ji) 

0 0 0 OL 

n£~s 

d = s dl”. 

4- 

= n^-r^-s 

i' J 

r^+1, i^ji+2, • • • 


d^yd^ 


and 1 and 0 are the matrices of suitable orders with all entries 
as 1 and 0 respectively* 


i 
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Following Lloyd (1952), the LS estimator of Y in 
this setup is given by 


{ 2 . 2 - 1 ) 


r 


e 

#N> 




= (A' Q A)"'^(A' 0 x), 

4S> rs> «s» ro fM 


where Q = D 

rsa ro 


-1 


That is. 


Q- O ... O 


0 Q Q V <0 4^ O 

A fsj 


<SJ 

m 

0 

.■m 


2 2 ••• gK 


and = 


Q . 

X 


~ +(ni-ri-l) 


- ( njL-rj^-l ) 




0 ... 0 

2 


( nj|^“*r^~l ) +( n^*“rj|^~ 2 ) ■“(nj|^~rj^~2) *.»» O 


0 


o 


-( s^+i)n 

> (s^tl)^ 


Hence, A' Q A 

|v> ro ^ 


— 

1/a^ 

0 

0 

:V“i 

o 

1/3.2- #«» 

0 

iV®2 

f 

j 

# 

0 

0 ... 

1/aj^ 

f — 

"b," ■ 



I 

! Q 
) 

_ 

^2 



(K+1)X(K+1) 


bj = c 


(i) 


ri+l 


■'i = "r,+l (ni-j+l)~^ , 


j=l 


.(i) 


r^+1 


a_. = aL~’ = S (n,~j+l)"'^, 
j=l 


K 


K 


Q = d+ S 

1=1 


, d = 2 and dj = <^5-1 for 1=1,2,— K 

1=1 


where 


Now, for a partioned . matrix 


R S 

rs> 


s' T 


(where R is non~singalar) , 


the inverse is (see, for example Rao 1973, p« 33) 


+ P F' 

ro 


ro oo 


-F E 

fo ro 


.-1 


Tn/ 


-E F 

rs> CO 


where E = T-S' r"^ S and F = rT^ S. 


.-1 


CO <s> ro fO 

Applying this we have 


ro (TO CO 


(2.2.2) (A'Q A) 

ro ro <s» 


-1 _ 1 
” d 


a^d+t^ 

b^bj 

0 0 0 

1 ~bH 

bjbj 

0 

a2d+b^ 

• 

... b2bj^ 

’ -^2 

1 

, 0 

0 

0 

0 

0 0 0 

t # 

• 

t 

j - - 

t 

-^1 

-b2 

... -bj^ 

f 1 

t 


Now 


(A' Q X) = 

• ro CO c<o 




/a 

^rj^+l/^K 


K ^i~®i / . •) f .) a) 

where P. = S [ S +Si^n -s ^ ^i^^i^ ^r +J .Substituting 

^ 1=1 i=rj+l i i , i 


^ (!) 

^1+ 2 +l/^ 

1 ^ i rj^+1 . 

,^( i) . - -u-( i) 


±=1 j =r^+l 


for (A'QA)“^ ted <A'Q xj in equation (2.2*1) we obtain for 

ro rOi rO' ,irO' 


^-1 

X s: X/2/'*® K 


< = 4 i!rV* ='* = . 


( 2.2.3) 



where 

( 2 . 2 - 4 ) 


j=l 

Note that for K = equation (2^2#'^) reduces to equation (2-1*4) 
obtained by Sarhan (1955). 

The variance covariance matrix of r* is (A'q A)**^cj^, 

tva «v» fs> 

From equation (2.2*2) we have 


( 2 . 2 . 5 ) Var(ep = [aj^+bJ/d] (i = 1, 2 * . ..,K), 


( 2 . 2 . 6 ) Var(0*) » aV<i 

and 

( 2 - 2 . 7 ) Covar(eJ,a‘‘^) = -bj^ 0 V<i/Covar(e J,et)*b^b^aVd( 2 / . * »^k). 

Hence - 

( 2 . 2 . 8 ) Var(e*-eJ) = [aj+a2+(bj^-b2) 

2 - 3 . Distribution theory . 


In this section, the distribution of the statistics are 
discussed, which are main tool for deriving the distribution of 
the LS estimators * 


Theorem 2.?.1 . +1 ' ' ' ^r^+l ^ * 

distributed. 

Proof. The jpdf of X(1 ),x( 2) , . - .,X(K) is 


are independently 


f ( x( 1 ) , x( 2 ) , . 



i=:l 


-,x(K)/©,,e2,.-t,8j^,o) 

^ X ^ 

Uj-t [l-e3^f^(x^^ 


r^f a 


^r^r®i 


^i 

—x 




. exp{-sj_(x^^2 -©iVcJ - S (x5>^-^e.)/a} 

i i j=r^+l J 


4i!l i®! (i = 1-2^ 


rK) 




For the sample (i s= l,2/«*»/K), considering the 

transformations. 


( 2.3-1) 




= (n^-j+1) (Xj^^-Xjf;^)/a (j = r^+2/ . . •/ «^-i> 


• • •/ n^-"S j ) 


.(i) 


and noting that > 0, the corresponding inverse transformation 

is 

+ Si' 


^ ^ ^ c? / ^ ® *{ ^ •! • 


Hence, 


n.-s. 

^ S Ji) I /-.(i) 


( 1 ) 


3 =r^+l 3 =r^+l i i 

Following standard methods, the jacobian of the transformatlcxi 


K sA a 

J = 71 


ni~r^-Si 


i— 1 Cn^^r^)5 


is 





The jpdf of 's (j s= r^+l, » . »#n£~s^/ i =1/2, is then 

obtained as 


, y(2) *.>/ y(K)) 


K 

= n 


■V* n • • 

i , ^ ..(i). 


i=l ^i* '‘’^i"'^i 




J=r^+1 


( 2.3*2) 


" il b7r(^"r7n 


Hi-Si 


. 71 [exp{-y!.^l] for y;.^^ > O. 

j=r.+2 ^ ^ 


( * 

This shows that 's are independently distributed with the 
following density functions : 


*^(1) 'y’ Cl-exp{-y/(n^-rj^)l]''^ e3q=(-y),Y.>;0 


ri+l 


( 2 *3*4) f / j ■\(y) = e3<p (-y) / y > 0 


Vi) 

j 


for j = rj|^+2/ rj^+3, . • *, n^^-s and i = l/2/».»/K« 


Note that, 


K ^i-^i 


^ J- X ( 4 ) ^ 

(2*3*5) 2 £ y\ = dp*/a, 

±-l j=ri+2 ^ 


which is free from yi , ( i = 1/2/^* #/K)# 

ri-t-x 

From eguations (2*3*2) to (2*3*5) and trensforroation (2*3*1), 
it is easy to see that ^r^+l * ** ** 

independently distributed* This completes the proof of the 


theorem* 



Corollary 2<>3>1* 2da*/a has a chi-square dis-trltaution with 2d 
degrees of freedom (DP) i 

Proof * By making the transformation = 2y\^‘^ ( j=r.-f-2/ — .-s./ 

J J ^ ^ 1 i 

( i) 

i = 1/2>*-*,K) in equation (2-?# 3), we see that Z\ 's are 

iid chi— squaare variates with 2'DFi Consequently/ 

K ^i**®! , ^ 

S 2 z'-. = 2 da*/a 

i=l j=r^+2 ^ 

has a chi-square distribution with 2d DF. 

In particular/ the pdf of := dcJ^/cT is given by 


(2 *3 *6) f(w) = w^"*^ T(d) / w > 0» 


Corollary 2«3«2« The pdf of (i = 1, 2/**»/K) is given by 

^r. 

[l-exp{-(x-e . )/cr} ] 

(2.3*7) f (x) = — exp{~{n/-r.)(x~ej^)/a} 

^rl+l B(ri+l/ni-ri)0 

for X > 6^/ where B(rj^+1/ n^— r^^) = (n^-r^-l) i/n^^l* 

P roof * The proof follows directly from the marginal pdf of 

or from equation (2*3^3) and the transformation given by 
equation (2.3*l)» 

The distribution of the LS estimator 0* of 0 i=l/ 2/ « * ./K) 
is given in the following corollary ; 

C orollary 2*3*3 * The pdf of = 1/2/***/K) is 

# Y < 0 

P2(y)/ f I 


( 2 *3*8 ) fv Cy) = 
^i 





where 


r . 
1 




, . - {-(n.-r.+j+d/b. )y}‘ 

(y) = G E (--l)^( .^)expCdy/b.} v i x i 

j =:0 J ^ 


h =0 h { (n^-r^+j+d/ti^^* 


P 


2 ( 7 ) = G E^(~l)^ ( j^)exp{-(n^~r^+j )y}/(n£-r^+j+(Vbj^)‘^y 


r ^+1 

G =: (d/b^)^/B(rj_+l,nj|^~r^), ^ (nj^~j+l)~^ 

j -1 


ja 

and 6 jj^ is given in equation ( 2 * 2 «> 3 )« 


Proof » The proof follows immediately on making suitable 
transformation and using Theorem 2*3*1 along with equations 
( 2 .3*6) and (2*3*7) * 

2 *4 • ml estimators of the parameters for K(> 2 ) samples case* 


Let X(i) (i = l/2/*..yK) be (K > 2) independent type II 
doubly censored samples from E( 0 ^/ 0 )* Then the likelihood 
function is 


L(02,6 2/ • lx(l)f x( 2 )/ •**/X(K) ) 


-rd. 


.(i) 


K 

TT 


i_l rj^! 


— [l-e ] e3?)[~ -ei> 

X X 


n^-s . 

XI (i) 


.+ E for > e^(i=l, 2 ,*..^K),a > O, 

j=r^+l ^ i"*" 

where d* = nj^^r^-s^^* For i = 1,2#.**,K make the substitution 
j(i) _ (x4^^-0^)/d (j = r^tl# » • n^— 


Zj “ ^ J = •^i"’ 



then we have 


»• l^l)^z(2)/»»«/ zCk) ) 


,* 


~d7 ( i) n _s 

K n . I 0 r.+l r. /■) i i r 

n [l-e ^ ] exp [-Si4,-3.- S. 4 ] 


i=l 


j=rj_+l 


for >0(1 = l/2|r«»«,K)« 

r±+l - 

Differentiating log L(©j^,62/ • • 1 2)/ » • z(K) ) w«r»to 

0 £ (i = 1/2 #*«j»/K) and a we get the following (K+1) likelihood 
eouations : 

z^i^ 

rj+1 . 

(2.4.1) r/[e ^ -1] - (n^-rj^) =0 (i = 1^2 ,.. .,k) 


and 


K 


(2 *4 *2) E [ E ^ - d* = 0. 

i=l j =r j^+1 i i 3- 

Equation (2.4.1) simplifies to 


z 


.(i) 

r.+l 


®i = 4!‘ii ^ ^ 

X ( i) 

(2.4.3) = - a mj., 

where = log {ni/(n^-rj^)} (i = 1, 2, .../K). 

Substitution of (2.4.1) in equation (2.4.2) gives 

(2.4.4) 0 = dhVd* = 

where 0* is givean in equation (2-2*3). 

The distributions of O and O^can be easily obtained by 
applying Theorem 2'*3^1 and its corollaries* 
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2 *5 • Comparison of ML and LS est±niat.ors # 

Epstein and Sobel (l954) showed that/ the LS estimators 
given in equation (2il*3) are the minimum variance unbiased 
(MVU) estimators of 6 and a respectively* It follows that 
6* (i = 1/2/***/K) and a* given in equation (2#2*3) are the 
MVU estimatojTS of (i = 1/2/*»*/K) and o respectively 
[see/Sarhan and Greenberg 196 2, p* 368 ]» 

Then from equations (2*2*"^) and (2*4*4) we have 

(2-5*1) a* = d''''^/d , = x^^l.-b-a"^* 

Using the equations (2*2*5)/ (2*2*6) and (2*5*1)/ the 
means, variances and mean square errors (MSE) of ML estimator 

A 

a are respectively 
(2-5*2) E(a) = d0/d* 

(2-5*3) Var(a) = dh^d*^ 

(2-5*4) MSE(a) = (d+K^)CTV(<i+K)^ 

Prom relations (2*2*5) and (2*5*3), it can be seen that 
Var(a) < Var(a*) • 

The relative efficiency E of O w«r*to d* is given by 

E = MSE(cJ*)/MSE(a) = (d+K) V(<i^+<iK^), 

Note that/ E is less than 1 whenever d > Since in 

general d is large, E is less than 1 for most values of K > 3- 


I 
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This shows that in general for K > 3, (7* is a better estimator 
than a • However, for K < 2^ the value of E is greater than 1 

ft ^ 

and the ML estimator a even though biased, is better than <7 * 
This agrees with the conclusions drawn by Kambo (1978) for 
K = 1 • Further, for fixed K, E tends to 1 as d tends to 
infinity. Table 2.5.1 gives E for K s 1(1)6(2)10 and 
d = 2(2)10(5)30. This table shows that for K > 3 and moderate 

A 

values of d, the estimator (7 is considerably better than c* 
The equation (2.4»'') also gives 

t 

E( 0 ^) = 0 j_+a(bjpA^) and Var( 6 j^) = {a^+A^/d}a^. 

Hence , 

(2*5.5) MSE(0j|^) s= { (aj,+A^/d) + (bj^— A^) \ O / 

where A^ = dm^/d* ( i = 1 , 2 , **.,K) • 

A 

In the general case direct corrparison of the MSB of 6 ^^ 
and 0* is difficult. But from equations (2*5.5) and (2.2*5), 
it is clear that, for r^ = 0 (i = 1,2, ....k) MSE(0*) < MSE(0^) . 

For the two. sarrples case, aj, = 5 , r^ = 0 , ^^2 ~ 

r =12 and s« = 1, the MSE of estimators are calculated from 
equations (2.2*5) and (2*5.5), and are given by 

MSE(0*) = 0.05 330^ < MSE(0^) = 0*800(7^ 

A 2 

MSE(05) = 1*43240^ > MSE(02) = 1-36810 . 


arid 



This shows that nothing can be concluded about the relative 
efficiencies of these estimators* The proper estimator out 
of these two is the one with the smaller MSE, Similar 
conclusions arc drawn by Kambo (1978) for the case of one 
population* 


A A 

For later us3/ wo now evaluate the variance of * 


From equations (2*4*3) and (2*4*4), we have 


®1-S2 = 

where m = m 2 « From Theorem 2»3*1 and equation (2*2*6) we 

get 


A A 


_ 2 


( 2 * 5 * 6 ) Var( 65 ^- 02 ) = {aj^ 4 a 2 +m d/d*'^}a 


2*6* Estimators of the parameters under the hypothesis 0^=0 2 =®' 


The LS estimators and the ML estimators of the 
parameters 0 and a under the hypothesis based on 

two independent type II censored samples are discussed in tnis 
sGction* 


2 *6 *1 * LS estimators of the parameters • In this subsection. 


all unspecified notations are as given in Section 2*2 with K 
Denote E(x) =8 T and Var(x) = cr^, 

«s> cvi ro »>> 


B = 

1 C(l) 

and T = 

i. e 


1 C(2) 

ra 

0 


^ isa ' ^ 


^ — 





where 
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Then, the LS estimator of 7 is given by 


( 2 * 6 . 1 ) 


2^ 

cvjO 


0" 


o 

o 


= (B' Q B)"^ (B' Q X), 

<s9 fo tN> ro «<s> 


where ( B' Q b) 

CO ro 


l/a^^ + l/a2 

b/a^ + b/a^ Q 


C 2.6.2) (B'qb)”^ = -1- 


CO CO ro 


Q' 


Q 


-(b^/a^+b^/a^) 


where Q* = IB' Q Bi = {(bj^-b^)^ + d( a^+a2^)/^1^2‘ 
Further^ 


(3' Qx) = 


CO <0 CO 


421/^2 \ 

Simplification of equation (2.6.1) gives 


(2.6.3) e* = i[Q{-3i 


b 

-•r.. 4..,-T-x b. b 2 .ir.fl 

1 • -1-} - (^ + ~){Pi+ S — 

^2 ^1 ^2 ^ i=l X 


+1 ^ r .,+1 




and 


• (1) „(2) 


K b ^ r.,+1 r^+1 

(2.6.4) a? = X- [-(^1 4 5^)C-^ 4 -Ir) 


^ ^2 '^l 


2 / . > 

Substituting for (P^ + .+l/^i^ ^ equation (2.6.3) 

by (2*6.4), and simplifying, we get the unbiased LS estimators of 


0 and a as 



J . 


( 2.6 .5) 

CD 

0 

H 

and 


( 2 .6 .6 ) 



(2.6.5) ej= 


( 2 ) 

2"* 


* = {(l:i-b2){x[J’j,-4^|^)+<S(ai+a2)a*!/{a(3l+a2)+(bj-b2^}. 


The variance covariance matrix of 2!^ is 


The equation (2 •6*2) now gives 


(2.6.7) Var(eJ) = Qa^Q*/ var(a*) = (a^+a^^O^^ 


and- 


CovarCeQ/a^) = -(a^b2+a2^1 * 

2 .6 . 2 . ML estimators of the parameters . 

The likelihood function of X{1) and X(2) is given by 


-d 




_2 nj^! a 


L(e,aix(l),x(2)) = JT^ . ^^-1 - 1— [l-exp{x^^jj_~e}/g] ^ 


. exp [- i{s.(x^^2c 2 (x^^^-0)}],x^^^ >e(i=i,2),a >. 0. 

a X “n^ jrrr^+l ^ ^i^^ 

( 1 ) ( 2 ) 

Without loss of generality, let + 1 / otherwise, we 

1 2 




can simply relabel the samples. Now y 


(2.6.8) L( 0 ,aix(l) ,x( 2) ) = Const [l-exp{~(x^\.,-e )/a}] 

<s9 ro ^ 

, [l-oxp{- 2 exp 

4Jil > e, a > 0. 


for 




where P = E [ e ^ 1 *="l+“2-^l-^2- 

1=1 3=r^+l 1.1 1 

According to zero or non-zero values of r^ and r 2 t we have four 
different cases* 

Case (l): r^#r2=0. The likelihood function (2*6*3) becomes, 

L(6,a lx(l) ,x( 2 )) = Const* esg) [- ^{ P+f(xj^^ ^-8 )}]/cr^ 

for > 0, a > O* 

It is clear that L is maximum for 6^ = and 0 ^ = P/d*» This 

has been also obtained by Epstein and Tsao (1953)* 

Case ( II) : r^^ > O, r 2 = 0* Now, the likelihood function (2*6»8) 
simplifies to 


( 2 -6 .9 ) 


L(e,aiz(l),z( 2 )) 

V.+l r. 


~ ( 1 ) 


Const 


- [l-c ^ ] ^exp[- E -s ^ ' 

i=l ^ ®i j=r^+l 


.(i) 




for > 0, 0 > O, 


where = (xj^^— 6)/0/ j = r^^+l/ » • and ± — 1,2* 

The maximizing equations for L(G,0 i z(l ) , z( 2) ) are 


(A ) 

(2*6.10) rj_/[exp - f = O 

and 

(2.6,11) [e^ 4 ^>i}-l] -f 4 iil-*“/‘’ + = O. 

Solving the equations (2-»6^10) and (2^6«H)/ we get 


(2.6.12) e . 


^r^+1 ” = P/d*. 


^ f 1 ) 

Not.G that e < 
o — 


-V\ * 

\+l' 

Case (III) ; := 0, r 2 ^ C. In this case the likelihood 


function (2.6*8) redaces to 


^ IT 

(2.6.13) L(e,aix(l),x( 2)) = Const. {1-qw} ^ 

CO 

for 0<w<l, 0<q<l/ 

where q = exp(-y/a) and w = exp {-(x^^^-e )/<?} • We maximize 
it w«r. to 0 first/ which is equivalent to maximising it w.r* to 
w- Unlike the Case (II), the maxltnum is not necessarily at a 
point which is less than could be at also. 

Towards this end, consider the function 

^2 f 

g(w) = (l~qw) w / O < w < 1/q/ 
which is plotted in Figure 2.6.1. 



FIGURE 2.6.1. A function related to likeldLhood function. 



r 

Note that l/q > 1/ g(0) = O, g(l) = (1-q) and g(l/q) = 0. 
Differentiating log g(vr) w»r» to w and equating it to zero, 
we see that the inai-ciinum of g(w) in the range 0 to l/q occurs 
at where 

( 2»6 *14) Wq = f/{q(r2+f)} • 

Note that may be greater than 1* 

Consequently, the differentiation of likelihood 
function w»r* to 0 gives the ML estimator, only if < 1# If 
Wq > 1, then the likelihood function will attain the maximum' 
w*r> to 6 at a point where Wq = 1, that is at expC-Cx^^^-© )/0’} =* 1* 

A ( i ^ 

Thus, the maximum in this case is at 6^ = x^ • We thus have 
two possibilities* 

Case (i) : w^ < 1# With same notations as used in equation 
( 2 *6 *9 ) , we have 



Solving these equations we obtain 


(2.#6.16) 0^ = log (-^)/ = (P-fY)/d*. 

2t 

Case ( ii) J Wq > 1 • In this case^ as shown above 0^ = 
Substituting this in equation (2*6 *13)/ we have to maximize 

(2*6 *17) L(0,0 |x(l),x(2)) = Const* a ^ ^ 

L/ <sa CO 

(for y > o, a > O,) 

w»r« to cr • For this, the maximizing equation is 
( 2 *6 -13 ) d*g + r2y/{ 1} — P = 0» 

The solution of this equation gives the ML estimator of (7- 
This can be solved by using Newton— Raphson method with P/d* 
as an initial value- 

In Case (III), the procedure of choosing proper ML 
estimator is as f ollov/s : 

Calculate the quantities a = y/log il+r^f) and 
b = (P-fy)/d*- Depending on the values of 'a' and 'b' we have 
three possibilities- 

Case (l) ; a < b* In this case, the relevant ML estimators are 
given in Case (i)- Since, the likelihood function attains its 
maximum within the pertinent range, w^ < 1 - Equivalently 

% > y/lof 

A 

where q ?= b- 



Case ( 2 ) 5 ■ a > In this case, the ML estimators are given by 

'' ( 1 ) A 

©o * and which is the solution of the equation (2»6»18)» 

Note that, for 


ficj) = d^'a + r 2 y/(e^^*^-'l) ~ P, we have 
f( 0 ) = ~P < 0 , 
f(a) = d*(a-b) > C 

and f'(a) = d^+rjY^ e^'^‘^/£a^( e^'^^-1 ) ^} > 0, 

hence equation (2*6«18) has a unique solution in (0, a) • 

Case (?) t a = b » In this boundary case, the ML estimators 
given in Case (i) and Case (il) turn out to be same with 

(1 ) A 

0o = ^ ^ = a =b. 

Case (IV) : r r 2 > 0. Using the same notations as in 

oppression (2«6-9), the likelihood function (2«6*8) can bo 
rewritten as 


( 2 *6 »19 ) L(e,a I z(l) , z( 2 )) = Const* a’’^ [l-e3<p{-z^^^^} ] ^ 


to 


• [l~e3p{- ^ ^ I “ ^ 


( 1 ) 
1 ' 
( 1 ) 




,( 1 ) 


for z"-^' > 0, <7 > 0< 

rj^+i - 


Consequently, 


( 1 ) 

log L(e,ai z(l),z{ 2)) = consti-d* log 0+r- log [l-exj){-z . 

#0 ca 1 


}] 


+ Z '2 log [i-exp{-*y/a-«j^j|j^} ] - P/o-fz^^^^^ 


.( 1 ) 



The corresponding equations are given by 








+ 1 = 0 


and 


=^24211 “=^<-42li* p. _ *4J 


.( 1 ) 
,+l 


On simplification/ those equations rodu-ce to 


(2.6.20) r^/[exp{ + r^[e3q5{ -1] - f = 0 

1 2 


.( 2 ) 


and 

( 2.6 . 21 ) 


( 1 ) 

^l^r^+1 


A2) 


m ^ m — - -*4^li-| + 'i* = O' 


Eliminating [® 3 cp{ -1] in (2*6.2l) by using (2.6.20)/ we 

get 

.( 2 ) 


(2*6.22) e = ^^2+1 ^ {l+r2y/(P--d ff)}. 


( 0^ 

Similarly/ eliminating [exp { z^ 4 . 1 ^“^} obtain 


(2.6 .23) e = - o log {l+r 3 _y/(d*a + fy-P)} • 

Equating (2*6.22) and (2.6.23)/ we have 


(2.6.24) = [l + — ^]/[l + 




P-d*a 


;]- 


d%+fy~P' 

This equation can be rewritt^ in the follov/ing alternative 


forms i 



( 2,6 .25) a = P/d*-r 2 y /[d*{ l+r^V(d Wy-P)} 
or 

C2.6.26) a = (P-£y)/d* + rjy/[d*{ l+r 2 y/( P-a*a) J - d*]. 

>1 * Equation (2.6.24) has unique solution for < 7 , and 
the solution lies in the interval ((P-fY)/d* P/d*). To this aid, 
let 

f(a) = [n-r2y/(P~d*a)]/[l+rj^y/(d*a+fy-P) ] . 

Note that/ 

f(P/d*) := -oo and f((P-fy)/d^) = exp {yd*/(P-'fy)} > 0. 

Hence/ the equation f(a) = 0 has atleast one root in this 
interval. Now, for concluding it has unique solution, it is 
sufficient to show that f(a) is monotonically decreasing# 
Differentiation of f(a) w«r- to a gives 

O 

f'(a) = ~exp (y/cr)y/g'^ 

^2^^* 1 rii -ir 1 

( d*a+fY^P ) ( P-d*cJ ) ^ "^ ( P-d*ci ) "^ ^ d*a+f y~P ) 

[l+r^y/ (d*a+fy-P ) ] ^ 

< O for every <7 e ((P-fy)/d*, P/d*). 

Hence, f(cj) is a monotonically decreasing function of cr in the 
above interval. 

Now, the ML estimator of 0 is given by one of the 
equations (2.6.22) or ( 2*6 .23), and frcroa expressicai (2#6#23)/ 



it is easy to see that 6^ < since elation (2*6.26) gives 

> (P-.fY)/d’''. The ML estimator of a is the solution of the 
equation (2*6 .24) » It is difficult to show by second- derivative^ 
that these values really give the maximum of the likelihood 
function • But some numerical galculations/ show that the 

A A 

likelihood function is really maximum at 0^ and <7^# 

Remark 2.6 > 2 « Eventhough r^ > 0^ r 2 0 in Case (W) , we can 
o?otain some results for other cases fiom this case by substituting 
r^^ and/or r 2 equal to zero. 

(i) If r^=:r2=0 then equation (2.6.25) gives = P/d , as in 
Case ( l) . 

(ii) If r3_>0^r2=0 then equation (2-6.25) gives * P/d , as ir. 
Case ( II) . 

( iii) If r^=0, jtr^O then equation (2*6.26) gives <7^ = (P-fy)/d , 
as in equation (2-6.16). 

(iv) If r^=0, r2>0 and 0^ = then equation (2.6.24) reduces 

to equation ( 2 »6 .18 ) • 

Remark 2.6.3 . A computer program for evaluating the ML 
estimators for all four cases is given in Appendix A » 


2 .6 .3 * 


the estimators 


In general the explicit form of the 146 estimators 
under the hypothesis is complicated. Hence the coitparison become 
difficult. However, for r^*r 2 = 0 the ML estimators are given 
in equation (2.1.12). This can be rewritten as 



2 

(2.6.27) Sq = min = [d0*+^S n^Cx^^^-S^) ]/d 


* 


Before deriving the MSE of the ML estimators, we first prove 
the following theorem- 

Theorem 2-6.1. Let U = n^^X^^ min( X^^ ^ , x[ ^ ^ 

and V = min(X^^\x^^^), 
where N = nj|^+n 2 , then 

( i) 2U/a has a di s t r ilxit ion , 

(ii) 2N(V-e)/a has a distribution, 

( iii) The random variables U and V are independently 
'distributed. 

Proof . The jpdf of X^^^ and X^^^ is given by 

= n^n 2 exp [- fc n^4^ 

for > 6, O > O. 

T.-t n n^x[^^+n 2 xj^^-N min{ x^^ \ xj^ ^ ^ and v = min{ - 


This is not a one to one transformation* Two inverse 
transformations are 

1) = V, x[^^ =: xV'b2 + ^ 

and 

2) 3^^^ = V, xj_^^ = vi/n^Wi 

where u > O and v > 0. The respective Jacobians of transforsa- 
tions are l/n 2 and %/xi^* The jpdf of U and V is then 



n- 

f(u^v) = — 2 63^ [- ^ {n^v+u+n 2 V-NG} ] 


2 1 

+ ^ C~ a C ^ u > o^ 

(2.6 #28) == {exp(~u/CT)/a} [n e 3 ^{ -N( v-e )/(7}/a], u> 

Then the marginal pdf of U and V are 


V > 0 

Oy V > © . 


( 2.6 .29) 
and 

( 2.6.30) 


f(u) = ^ e for u > O 
_ H(v~e) 

f(v) = ^ e ^ for V > 0 


respectively. The equation (2-6.29) shows that 2U/a has a 
distribution. Epstein and Tsao (l953) also obtained the equation 
(2.6*29) by using a direct argument. Similarly/ 2N(V-6)/a has 
a distribution. This can be derived by direct argument as 

well. From equations (2.6.^)/ (2*6.29) and (2'.6*30)/ it follows 
that U and V are independently distributed. This completes the 
proof of the theorem. 

\ A f \ 

C orollary 2.6.1* With and 0^ as defined in equa.tion (2*6.27)/ 
( i) 2N(0^-0)/cr has a distribution/ 

and 

( ii) 2d^a^a has a X^ „ distribution. 

2(d''^-l) 

Proof. The ML estimators given in equation (2*6.27) can be 
written as 


(2.6.31) ©Q = V, = l26a*/a + 2U/0] d/2d* / 

where U and V ana as defined in Theorem 2*6^. Now/ part (i) 



of the corollary follows from Theorem 2i6»l* By Corallary 2*3*1/ 
^ 2 

2da /a has a distribution* Since U is a function of 

( 1 ) (2) 

and X^ , on applying a result similar to Theorem 2#'3#1^ 
we see that a* and U are independently distributed* Consequently, 
2d^a^0 has chi-square distribution with 2d-{-2 = 2(d*-l)DF and 
the corollary follows • 

This corollary gives the following results : 

E r2N(e^-6)/a] = 2 implying E(6 q) = G-kt/N, 

Var [2N(G^-6)/a] = 4 implying VarCG^) * cr^/N^, 

E [2d*ajy'a] 2d*-2 Implying E(q^) = (d*-l)a/d* 

and Var [ad^'a^/a] = 4(d*-l) implying VarCa^) = (d*~l )a 

By these relations we have 

(2.<.6.32) MSE(0^) = 2aVN^ 

and 

(2.6.33) MSE(a^) = 

In equation (2.6.7) for r^ = r 2 = O, the variance of the LS 
e st imat ors are 

( 2 .6 .34 ) Var(e5) = n^+n|)^N^} 

and 

( 2 .6 .35 ) Var (a^^) = a Vf d*-N V ( nj+n^ ) } * 

Eventhough is a biased estimator of a, from 
equ.ations (2.6.33) and (2.6.35) it follows that MSE(a^)<MSE( a J), 
where 0 * is the unbiased estimator of d given by equation (2.6.6) 
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Further, as ej<pected MSE(( 7 *) < MSE(a^) where a * is given by 
equation (2»2#3)* This can be seen from the equations (2<*2»6) 
and (2*6 *35) and the fact that N^/Cn^+n^) <2* 

Direct comparison of MSE(0*) and MSE(0^) secans to be 
difficult, but for n^^ = n 2 = n, the equations (2*6*32) and 
( 2*6 *34) sinplify to 

MSE(e^) = a^/2n^ < iMSE(0*) = d*crV[(^*-2)2n^ ] . 

Also for nj^ = 5, n 2 = 15, s^ = 0 and S 2 = 0 the equations 
( 2 *6 *32;) and ( 2 *6 *34 ) give 

MSE(0^) = 0.005 > MSE(0*) = 0*00435 

This shows that, nothing can be said regarding the preference 

^ A 

of the one estimator over the other among 6^ and 0^* However, 
one can choose tho estimator with the smaller MSE* 


\ 



TABLE 2*5#1* The relative efficiency of 
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CHAPTER III 


TESTING OP HYPOTHESIS ABOUT LOCATION PARAi-iETERS AGA-ItTST 

ON &- SIDED ALTERNATIVES 

3#1* Introduction and test statis-bics » 

In this chapter a test statistic frdm type II 
censored samplesis proposed to test the equality of location 
parameters of two exponential distributions against one-sided 
alteimatives- The common scale parameter is assumed to be 
unknown* The null and the non-null distributions of the 
proposed test statistic are obtained* Some critical points 
and some values of power are tabulated* An approximation 
in terms of Student's t distribution for the null case is 
studied* 

Let be two independent sairjjles 

ri+l t ^+2 ^±"^1 

from EiQ^fO) (i = 1,2)* Consider the problem of testing 

Hq : = ®2 one-sided alternative hypothesis > © 

or H' : ®2' ^ unknown. 

For the right censored samples (case r^ “ ^2 ~ * 

Kumar and Patel (l9 7l) have proposed a test statistic (KP test) 

for testing H^ against H^ i 0 2 * Weinman et (1973) 

extended it for testing H^ against one-sided alteraatives* 

o 

Against their test statistic is/ 

= (X^^^ ~X^^"^ )/a* , 
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where a ' is the pooled estimator . of o, given in equation (2-2*3) 
They had obtained the null and non— null distribution of as 


( 3-1 -l) P [Wq < clH ] = 


n^(l-n 2 C/d)"' /n for c < 0 
l**n 2 ( 1+nj^c/d)** /n for c > 0 


n^ eT<p(-n2'!P)(l~n2c/'i) /i^/ c < 0 

(3-1-2) P [Wq < c|<pj = Q^(c^l0)-n2e''^’’’ Q^{ I n^c/d}/n 

+ n^^ esp(-n29>)L^{ c^in2c/«3l}/n/ c > C 
respectively^ where n = n^+n 2 / d = ’^l'^^2"’^l*'^2"'®l”*®2"‘^'' 


OO ^ 

Q (xls) = 7 y^^^e dy/(p— l)l, x > 0, p = 1/2/---/S > ~1 

^ X 

X - 

and L (xls) = / y^"'^e”^'*’^^dy/(p-l)I , x>0^ p = l/2,»-» - 
^ o 

The expressions for the critical points obtained by 
Weinman et al» (given in Section J -4) are actually intearchanged, 
although his tabulated values are correct- By using the 
notations as in Section 1 -4# the correct critical points of W 



if n^(l~a) < 

\ 

if- n'j(l*-a) > 



All the foti^olhg statistics are based on intuitive 
giX)unds • On the basis of LS and ML estimators/ we propose the 
following test statistics for testing against ; 

"^LSE ~ of SE(G*— e*) 

and 

'^MLE = of SECe^^-e^)/ 

’¥t ^ 

whore and 6^^ denote the LS estimator and the ML estimator 
of 0^ (i = 1/2) respectively* Using expressions for 9^ and 0^ 
given in equations (2 *2 *3) and (2*4*3)/ and expressions for the 
variances of (©^-S*) and ( 8 ^~ 82 ^ given in equations (2*2*8) and 
(2*5*6)/ these statistics reduce tO/ 

i 

(3.1.3) IlSE = 

and 

1 

( 3.1 .4) T^e = . 

where a* is the LS estimator of a given in equation (2*2*3)/ 

r^+l 

b. = S (n/-j+l)'"^/a^ = E (nz-j+l)"^' and 

j=l J=1 

= log {n^/(n^-r£)} for i = 1/2* 

Both of thtse statistics against one-sided alternative are 
equivalent to the st^istic ■ 
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Against very large values of T lead to the rejection of 
that is/ we reject against if T > c^f where c^^ is 
determined so that p[t > = tt»and oc is the chosen level 

of significance of the test# Similarly^ against very small 
values of T lead to the rejection of H^# We thus need the null 
distribution of T for finding the criticeil point c^, which is 
derived in Section 3 #2* The non-null distribution is obtained 
in Section 3 #3# In Section 3*4/ the exact and approximated 
critical points are evaluated# Power function and its approxima- 
tion are studied in Section 3*5# 

3*2* Null distribution of the statistic T • 

We first prove the following lemma : 

Lemma 3#2#1 « Suppose and Z 2 are two independent random 
variates with pdf of as 

r. 

(3*2.1) iz) = {l-exp(-z)} exp£-(n^-r^)z}/B(r^+l/n^-rj^) /Z>0 . 

i 

Then the pdf of Z » ^i“^2 giveai by 

(3*2.2) g( z) = 

2 

where H = JJ £ 1/B( f 

■ , ^ - • 

B(p,q) » / u^^(l-u)^^ p > O, q > O. 
o 


2 . (n2--r«+j)z 

H E /)B(r-+l,f+j)e "" r z < 


^1 . r- -(n--rj+j)z „ > n 

H S (-1)^( ?‘)B(r2+lrf+j)e ' ' 

j=0 ^ 





Proof. The jpdf of and ±s 


f(z^,Z2) = H(l-e ^( 1 -e ^ expC-Cn^^-r^^) 2^-(n2~r2) 22? 

for 0 < Z2 < 

Malcing a transformation z a: Zj|^— Z2 and 22 = we see that the 
range of Z2 and Z are ^2 ^ max( 0 , — z) and < z < 00 respectively 
The inverse transformation is = 2+22 and Z2 = Z2* 

Jacobian of the transformation is 1, and the jpdf of Z and Z2 is 

X* x* 

(3.2 .3) g(z,Z2) = H [ l-exp{ -'( Z2+2)} ] ^ { l-expC— 22)} ^ 

• exp [-f22-(nj^-r^) z ] / 22^ max(0,-z) /-<»< 2 < ^ 
The marginal pdf of Z is thus 

r- 00 

/ g(z/22)d22 for 2 < 0 

-2 *" 


( 3 *2 .4) g( 2) = 


S g(z, Z2)dz2 for z > 0 . 


L o 


For 2 << 0 ^ on putting y = 2+22# g( 2) becomes 


-(hi-ri)z « ~fz„ -(25+2) r. 
gC2) = He ^ ^ ^ ® ^ ^ ^ ^^2 


— 2 


= He 




o 


Now expanding by binomial expansion and simplifying, 

we get „ , , . ' ' 

g( 2 ) = H E f) e ^ ^ e (1-G^) dy. 

j=0 ' ® 



Substituting u = and simplifying, we obtain equation (■?i2*2) 
for z < 0, For z > 0, we expand [ l~e ^ ^ of equation 

(3i>2*3) as a binomial sum and get the corresponding result of 
equation (3*2;2). 


-C orollary 3>2>1 . If Y = Z + g), then the pdf of Y is 


( 3 • 2 *5 ) g( y) = 


r„ (n2-r2+j)(y-4^ 

Jr At,/- .. for y < <p 


H E /)B(r.+l,f+j), 

js=0 J 


1 . i ^1 ^ - ( n- -r- +j ) ( y- 

H E (— 1 ) ( . )B(r^+l,f+j )e for y > (p 

j=0 . J ^ 


Proof * The corollary follows from Lemma 3#2*1. 


Theorem 3 «2 •! > The null distribution of T [v/here T is defined 
by equation (3*1 -5)] is 


' ^2 . r 

H E .^)B(ri+l,f+j){l~h,( j)cr^/dh,( j) for c < 0 

, j=0 J 1 2 2- 

. r 

1-H S .^)B(r 2 +l,f+j)(l+h^( for c > 0 

j=0 ^ 


where hj^(j) = (n^-rj^+jVd (i = 1,3, d = f-Sj^-S2-2, H and f are 
as in Lemma 3-2*l* 


Proof * Note that given 

in equation (3.2^) . Then under : 6 ^ = 62 , we have T * d^W, 
where Z = %-Z2 ^ =? dJ^cx* By liieorCTi 2 #3*1, Z and W are 

independent with pdf of Z giv®i in Leanma 3*2*1, and that of W 
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*2*2 * With the notations used in Theorom 3*2-1, we 
have the identity 


( 3.2*8 ) 



B( r‘2+1/ f+j ) 
“■"dhJIljl 


+ s i-i)H 

■ j=0 



B(r^+1, f+j ) 

— ah"(T) 3 =!• 


Proof. On using the fact that f(t) given in ejquation (3*2»7) 


is a pdf, wa get the required result. 


This corollary can be rewritten as an interesting 
Gominatorial idvsntity, namely 


1 

j r:0 



B(r2+l/f+j ) 



B ( r^i^ +1 , f+j ) 


2 

= TT B(rj+l,ni-r-0 , 

i=l ^ ^ 

wherfi f = n^+n2“r]_~r2. 

3*3# Non~null distribution of the statistic T» 


We next derive the distribution of T under the 
alternative hypothesis * ®1 ^ ®2* needed for 

studying the power function of the test* 


Theorem 3 *3 .1 . The non-ntill cdf of T for f = (0j^--0 2 V‘^ ^ ® 
is given by 


(3.3.1) P [T < cif] = 


P- (c| <^) f G < O 




where 


=H S r)B{rj^+l,f+j).e3?>I-h2( .£l->h2(3)c} /dh^j), 

^ j=0 ^ 



F 2 (c|(p) = Qd^C3_|0)~H^^(-l)J(^l)B(r2+l,f+j) .exp£h^(j)d«))} i 

j=0 J 1 ! 

. ej(p{-h2(j)d<p} #L^{aj^|h2(j>c}/dh2(j) 

and the remaining notations are given in eguation (3*1 *2) and 
in Lertttna 3 *2 *1 * 


Proof. Let (i = 1/2). Then follows the 

distribution given in equation (3.2*1)/ and Y = Z,-Z^+(}) 

X M 

( 1 ) ( 2 ) 

= C^r +l"‘^r +1^/*^ the pdf given in Corollary 3*2.1, Clearly^ 

1 2 

T can be written as dY/W, Using similar results as in the null 
case, we now have the jpdf of Y and W as 


where 
£ 3 ^ ( Y/ w) 


f(y/w) 


fj^(y/w) for w > y S ^ 
±2(7,^) for w > 0 / y > <p. 


H 


^2 

2 

j =0 


(- 1 )^( 


r^ (n^-ro+j )(y-<p) 1 

,^ ) B( r^ +1 , f +j ) e ^ ^ e\r^" V ( <3-1 ) I 


and 

^1 

f,(y/w) = H 2 (£ 
j=0 

Making the transformation t = dy/w and w = W/ the corresponding 
inverse transformation is y = wt/d/ w = w. The iregion w > O, 
y < <p is transformed to w > 0/ wt < d^. Similarly, the region 
w > 0/ y > is transformed to w ^ O, wt > dfi Note that, the 


■f -(n-,-ri+j)(y-^) „ j 1 

1)^( ^^)B(r2+l/f+j)e ^ ^ e'^w^ V(<a“l) !• 
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boundary y = <J> is "transformed to the hyperbola wt = d<!P» (see. 
Figure -3«l)* The jacobian of this transformation is w/d, and 
the jpdf of T and W is 

gj^(t,w) for w > 0, wt < d^ 


g(t/w) = 


g 2 ('t,w) for w > 0, wt > d^^ , 


where 


2 r. d 


and 


g^(t,w) = H E (“l)'^( . ) B( rj^+1, f+j ) e 
j =0 ^ 


e w /dl 


A r- -h- ( j ) ( wt-d a 

g2(t,w) = H E (-irC ?)B(r2+l/f+j)e ^ &^\r /d I . 

j ssO ^ 



FIGURE Joint pdf of (T,w) for ^5 O# 

The cdf of T upto the point c is the integral of joint 

density of (T,w) over the shaded region sis shoim in Figure 


This is given by 

<0 C 

P [T < Cl = / (/ 


and 


o 

00 0 


= (say) for c < 


dp/w 


o 


P Ct < cl <p] »i (/ g^ dt)dw +/ u dt> dw 

- A df/c o 


o 
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d<P/c c oo G 

+ / 91 dt)dw + / (/ g- dt)dw 

® ® d^/c dg/w ^ 


^2'^^3'**^4'^^5 for c > 0 * 


Simplification of (i = 1 / 2 ^ 3 # 4 / 5 ) involves lengthy cypress Ions, 

although the method is straight forward* Simplification of one 


factor I. is illustrated below : 


1 c 

14 =/ (/ dt) dw 


o o 


H E ,^)B(r-+l,f+j) 

j =0 J ^ 


^^h 2 (j)d«l> c^ , 


rTTT ^ ^ ® 


2(1 )wc ^ ^ 

rd-D'T 


^ i ^9 

= H E (-1)J( r) B( r. +1, f+j ) 
j =0 J 


-h,( j )d<}) 


=1 d-l -Mtl-h 2 (J)c) c 

dw - ^ TdiTTf 


^2 r -■h 2 (j)dg> 

H E (-1)^( ?)B(r 3 ^+l,f+j) ^ ^h ' rn ’ 

j =0 ^ . 2 J 


Similarly, 


L = H E (~l)^( . )B(ri+l,f+j)e 
j =:0 ^ 


“h2(j)d<P {i-.h„(j)c}~‘Vdh,(j), 


I 2 = H E (~1)^( 


-h 2 (j)d(P/^ 


■^2 . r^ “h5(j)d^ 

= H E j )B(r^+l/ftj) £l~e } Q^(c^l0)/dh2( j 5 


■ j .=0 


and' 



H 


'1 . r. 

2 h 

i=o j 


“^TTjl — r^d^c^io) 


- exp {h 2 _(j)d<J>} Q^{c^|h^(j)c} ] . 


Combining these expressions and using Corollary 3.2 .2, the cdf 
of T as given in equation (3.3.1) is obtained* 

I'or calculation purposes/ we need some further 
simplifications- For this, note that 


Q^( xl s) 
and 


2 ejp{-x(l+s)}{ x(l+s)} Vt Jl(l+s)‘S / s > -1 

j=o 


X /d I for s = 1 

d-1 . 

[l- 2 exp{-x(l-s)}{x(l~s)} VJ jy(i-s)*^ for s /t 1 

j=0 


The proof of these results and a computer program for evaluating 
these functions are given in Appendix B- 


Substituting these results in expression (3.3.1) we 
obtain the cdf of T as 


(3.3.2) P [t < c|<p] = P-(c|9>) = 


Gj^( cl 9) for c < 0 


GjCcl 


for c > 0/ 


wnore 


G. (clf ) = H S ( 
2=D 


2 . r„ B(r-+l,f+j) esq) {-h,(j)d^} 

) — : —3 


<ah2(j)Cl-h2(j)c} 
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^ i=o i« 1 “ j^o ^ j' ^ " - d t q n i — 

d~l cj e5^(-c ) ^2 r- B(r,+l,f+j) 

i=0 i!{l+h^( j=0 j dh2( j){l~h2( j)c}^ 

~c. d~l 

• [exp{-h (j)d<p} - e ^ E Ti ] 

i=0 ^ ‘ ^ 

^ ^ (j = 0,1, 2/ • • •,r 2 )' However, for h 2 (j)c = 1, 
the ( j +1 ) th term of the last factor in simply becomes 

1? 

H(~l)^ i B(rj^+l,f+j) e>p {-h 2 ( j)d<p}c^/{dh 2 ( j)dl} • 

Note that, for <p = 0, equations (3#3*1) and (3;»3*2) reduce to 
the null distribution of T., given by eqpiation (3»2#€)# It is 
easy to show that, for r^i^ = r 2 = 0, the equation (3*3*2) reduces 
to equation (3«1»2), the non-null cdf of the test statistic 
considered by Weinman (1973)» 

Usually, the null hypothesis is also taken as one-sided 
against one-sided alternatives* Thus, if we take H* : < ©2 

and test it against > 0 2 ^ then we also require the power 

function for negative values of (p • The following corollary 

gives the necessary distribution theory results* 

% 

Corollarv 3 *3 -I* The non-null cdf of T for <p < O is 


(3*3*3) P [t < c|<p] = P^Ccif) = 


M^(clf), c < 0 

M«(cl<p), c > 
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where 


Ml(c|q.) = 1-Q^(cil0)+H S (-l)^(^2>B(r.+l,f+j) 6xot-h,(j)d<« 

j aO J * ~ ^ 

• r, 

•Q<i{Cj^l~h 2 (j)c}/dh 2 (j)~H 2 (-1)^( *)B(r„+l,f+j) 

^ j=0 J 2 

.exp{h^( L^£c^l-h^(j)c}/dh^(j), 

M 2 (cl<p) = 1-H s (-l)^( j )B(r 2 +l,f+j )exp{hj^( j )d^P}{ I4hj( j )c}*' /dhj^(j 
j =0 I 

and the remaining notations are same as in Theorem 3 * 3 *!• 


Proof * The proof of the corollary follows on the same lines as 
that of Theorem 3 •3 *1, but the cdf of T upto the point c is 
the integral of joint density of (T,V7) over the shaded regicsi 
as shown in Figure 3*3*2- instead of Figure 3-*3*l* 



This is given by 


P [t < cicp] = 


00 oo 1 Q 

(-T g2d.'fc)<ivr + / (/ g«d.t)clw 
o o o c ^ 

“ O oo d(f)/w 

+f (/ g 2 dt)dw +X (/ g|^dt)dw] for c < O 


Cl c 




oo OO 

1-X (/ g 2 dt)dw for c > 0» 
o c 


On simplification, this gives the equation (3»3»3)» 

Rgntark 3 *3 »1 » Aitaough due to the complicated expressions 
involved, it is not possible to show that the test based on T 
is unbiased, yet we feel that this test is unbiased* Limited 
calculations carxied out in Section 3*5 strengthen this feeling- 


Remark 3 * 3 *2 * Critical points and the power values for testing 
Ho I = ®2 H^ : ©i < ©2 obtained by relabelling 

the samples* 

3*4* Exact and approximate critical points of T * 

3 *4 *1 • Exact critical points * The critical point is 
obtained by solving the equation 


(3.4*1) P [t > %|Ho] = 

whore a is the chosen level of significance of the test* Prom 
equation (3 *2 *6) it is clear that Cj^(< O) is the solution of 

^2 ■ r 

1-H E (-1)^( /) B(ri+l/f+j)fl-h2(j)cQ^} /dh2(j) = a, 
j =0 ^ 

if Po ^ l~a, where 
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^2 

= P tT < 0|H„] = H E (-l)^(7)B(r,+l>f+j)/ah,(j). 

j=rO J -*• Z 

Similarly, if then o^(> 0) is the solution of 

r 

Note that, if = 0,then 


= 1*"B( r2+l/ f )/B( r2+l/n2~r2) ♦ 

Now, if Pq < l~a, then is given by 

% = d - l]/n^ > o. 

Similarly, if r 2 = 0 and 


Pq = B(rj^+l,f )/B(rj^+l,n^-rj^) > l~a, 
then is given by 

= d [l-{pQ/(l-a)}^'^'^yn2 < Oi 

In particular, if r^^ = r 2 = 0, then P^ = n^/(n^-m 2 ^ 


'a 




n- 


1/d 


n. 


d_ Tr 

“i 


- N 


i€ < n2(l~ix)* 


Other than these cases, c^. can be obtained by evaluating P^ 
and then applying Newton — Raphson method to the pertineit 
equation* Note that if < r 2 tha:^ it raay be more convenient 
to calculate P [t > critical points are tabulated 
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in Table 3*4-l and Table 3*4-2 for a = 0*05. These are tabulated 
for following combinations of sanple sizes and censoring patterns? 

(i) for n3^=n2=10,(r3L/r2)«(0,o),(0,l),(0,2),(l,0),(l,l),(2,0) 
and d = 4(2)16 in Table 3*4*1, and 
( ii) for n3_=10,n2=6(2)24,(ri,r2)=(0,0),(0,l),(0,2),(l,0),(l,l), 
(2,0) and d = 12 in Table 3*4*2. 

Note that and S 2 appear only throigh d* Due to many variables 
involved, it is not possible to give a large table of critical 
points* In these tables, we have provided mostly those valnes 
which were later used for studying the power function of the test- 


initial value for solving equation ( 3 *4 *1 ) can be taken 
as an approximate critical point, given in Section 3*4*2* 


3 *4 *2 » Student* s t. approximation* The asymptotic distribution of 

sample quantiles ^^+1 vrall known (for example, see 

1 2 

David 1981, p. 255)* We however use the exact mean and variance 

of Z = under H • Prom the results given in Section 2*2, 

we have 


( 3 . 4 * 2 ) E(z) = BC, Var(z) = Ac^, 

r^+l 

where B = bj^-b 2 , A = a^4a2/ 

3 “1 

bi+l 

bi = 2 (n^-j+l)"^ . (i = 1«2)- 

j =1 

Using the asymptotic normality we see 


(n^-j+l) 


that' 


-2 

* 
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t f AA 2 AN(0, 1) i 

NOW/ following Tiku (l98l)> we can approximate the null 
distribution of 

(3.4-3) = (T-B)//A 

as Student's t distribution with d DP. Thereto re, the exact 
critical value is approximately given by 

(3-4.4) c^ = c^Va + B=!C^/ ^ 

where c^^ is the upper ath percentile point of Student's t | 

distribution with d DP. ' 

i'' 

E 

3 . 4 - 3 . Normal approximation » We now consider a normal | 

approximation for the null distribution. For this we need the | 

i 

mean and variance of T. [ 

From Corollary 2-3*1/ W = 2d0*/<7 has a H<asco 

( 3 - 4 . 5 ) E [l/W^j = (d-m-l)l/[2"‘(d-l)l]/ m = 1/ 2, , . (d-l) - 
Now/ T can be written as 

( 3 . 4 . 6 ) T = 2dV(cfW). 

.From the independaace of Z and W, and using eguation (3.4.2’)/ 
we get 

( 3 . 4 . 7 ) E(t) = dB/(d-l)/ VariT) = d^f (d-l)A+B^}/{ (d-1 ) (d-2)} . 
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We now use the fact that 

{T-E(T)}/{ Var(T)}V^ S AN(0,1) , 
where E(T) and Var(T) are given in equation (3*4»7)* Hence, 


(3*4-8) Cq_ = c^^{Var(T)}^/^ + E(t) = c*, 

V I#, 

where is the upper ath percentile point of normal 

d i st r ibut ion • 


The exact values and the approximated values c* and c* 
obtained from the equations (3-4-4) and (3-4-8) are tabulated 
in Table 3-4-3, for a = 0«05 and for some selected values of 
2 * ^ 1 ' ^2 basis of this and some other 

calculations, it is observed that, in general normal approxlmaticxi 
value c* is bettor than cj for > r 2 , otherwise cj is better 
than c 2 - 


3 -5 * Power function and its approximation- 



The power of the test T for testing H* : 

against 

Hf : > 02 is given by 





for 9 > ^ 0 

( 3 .5 -1 ) 

P [T > %I<P] = 


for < 0, 


< e 


2 


where is the exact critical point, P^(Cq^I^) and p 2 (c^l ^) 
are given in equations (3-3-2) and (3-3-3) respectively^ For 
a = 0.05, tho power of the test T is tabulated for various 
values of 9 in Table Table 3^fl2 and Table 3.-5 -3 for 



w 


the following combinations of sample sizes and censoring 
patterns : 

(i) for n]_=n 2 =l 0 ^d=l 6 and (rj^,r 2 )=( 0 , 0 ) , ( 0 ,l)^ (O^ 2 )^ (l, 0 ), ( 1 , 1 ), 
( 2 , 0 ) in, Table 3-5.1, 

( ii) for n^=n2=10,d=:4( 2)16 and r^=r 2 =l in Table 3*5.2, and 
( iii) for n^= 10 ,n 2 = 6 ( 2 )l 6 and rj^=r 2 =l in Table 3*5*3. 

The variation in power due to different combination . 

of r^i^ and r 2 is showed in the Figure 3*5*1. From this and some 
other calculations the following points emerge* 

(a) Table 3.5*1 and Figure 3*5.1 show that the test T is more 

sensitive for variatioi' in r- compared to variation in j 

r 2 • Consequently, it is desirable to take greater care 

in handling first sample so that censoring on the left is j 
reduced to a minimum for this sair^jle* i 

(b) It is clBar froxti the Tablo 3 • 5*2 that fof £xxed | 

f 

and r -5 the power of the test is relatively less affected | 

by variation in right truncation which is repres^ted by ^ 

variations in d- | 

(c) Table ^.5*3 shows that for fixed r^,r 2 and d, the power 

of test is not much affected by increasing the sample size* 

Since the power function given in equation (3*5*1) is 
very complicated, for large d and moderate critical points (c^) 
the following normal approximation is suggested. 

Let U = ^ 

normal approximation of chi s^are aistrltation yields 
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•jf ^ 2 

0 -^N(c7^(7 /d) and foir luodeirate value of the effect of o^Cf^ 

vXr U 

is not negligible in the linear combination defining Then 

from Theorem and Corollary 2«3*1^ we see that under 

E(u) = (<P+B-CQ^)a and Var(u) = (A+c^/d)0^, 
where <p = ( 0 ^— 62)/0» Therefore^ 

p [T > =«'»]= p ’/<’*> c^if] 

= p [u > Ol(pJ 

( 3 - 5 . 2 ) = 1 - f(b*), 

where b* = -(<P+B-Cq. )/(A+ c^/d)^'^^ cand f(x) = / dy/v ( 27i ) , 

mmCO 

The exact power [left hand side of equation ( 3 - 5 * 2 ) j and the 
approximated power [right hand side of equation (3 -5 *2)] are 
tabulated in Table 3 - 5-4 for a = 0 . 05 . As can be seen from this 
table, the normal approximation is quite satisfactory in this 
case- However, some calculations performed for small values of 
show that this approximation is not that good for such 


situations 
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TABLE Exact upper critical points c„ of the test 

- I I r- ' ' " - vir 



statistic T 

for a = 

M. 

□ .05 and n^^ 

=n2=10. 

— 









(0,0) 

(0,1)^ 

(0i2) 

(1,0) 

(1,1) 

(2,0) 


4 

0.3113 

0.1901 

0.0818 

0.6401 

0.4780 

1 .0100 


6 

0.2807 

0.1776 

0 ,079 3 

0.5575 

0»A'2B0 

0^614 


8 

0.2668 

0.1717 

0 .0780 

0.5209 

0.4054 

0.7960 


10 

0 .2589 

0 .168 3 

0 .077 3 

0.5002 

0.3925 

0 .759 3 


12 

0.2538 

0.1661 

0 .0768 

0.4870 

0.3842 

0.7358 


14 

0 .2503 

0 .1645 

0.0765 

0.477S 

0.3734 

0.7195 

( 

16 

0.2477 

0.1633 

0.0762 

0.4711 

0.3742 

0 .7075 

i 

1 

TABLE 3 ^4 

•2 . Exact upper 

critical 

points ( 



Df the test 


statistic T 

for a = 0.05, n^=lO 

and d = 

12. 







. 

"2 

(0,0) 

(0,1) 

(0,2) 

(1,0) 

(1,1) 

(2,0) 


6 

0.2194 

0.0952 

-1 .4357 

0.4436 

0.2928 

0.6863 


8 

0.2396 

0.1^70 

0.0321 

0.4693 

0 .3474 

0 .7154 


10 

0 .2538 

0.1661 

0.0768 

0.4870 

0.3842 

0.7358 


12 

0.2644 

0.1876 

0.1097 

0.5000 

0 .4110 

0 .7506 


14 

0.2726 

0.2043 

0.1351 

0.5099 

0.4315 

0.7618 


16 

0,2792 

0.2176 

0.1553 

0.5178 

0.4476 

0 .7707 


18 

0.^46 

0.2285 

0.1718 

0.5242 

0 #4606 

0.7779 


20 

0.2891 

0.2375 

0 .1856 

0.5295 

0.4714 

0 .78 36 


22 

0.29 29 

0.2453 

0.1972 

0.5339 

0.4805 

0.7835 ; 


24 

i 

0.2962 

.0 .2519 

0 .2072 

O .5377 

0.488 2 

0.7927 
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TABLE 3 *4 *3* Comparison of exact Student"' s t approximation 

( and normal approximation (c*) critical 
points of the test T for a U!)»05> 


ni=^2= 


.2519 
.1652 
.0766 
.48 20 
.3811 
.7 270 


10, d = 


.2505 
.20 74 
.1518 
.4296 
*3744 
.*6240 


nj^=n2=15, d 


22 


a 


-2632 
.2197 
.1712 
.4604 
.39 35 
.682.7 


.1618 
*1088 
.0548 
.2939 
.239 3 
.4358 


*1620 

.1318 

.0941 

.2746 

.2374 

.3908 


_!L 

.1665 
*1359 
.1002 
*28 55 
*2440 
*4111 







TABLE 3«5*1« Power of the test T for testing against 



for a 

= 0.05, 

n^=n2=10 

and d = 16 

* 









(0,0) 

(0,l) 

(0,2) 

(1,0) 

(1,1) 

(2,0) 

-0 .50 

.0003 

.0003 

#0003 

.0009 

.0009 

.0017 

1 

o 

« 

o 

.0009 

,0009 

#0009 

.0020 

#0020 

.0035 

-0 .30 

.0025 

.0025 

#0025 

.0046 

.0046 

.0071 

~0 #20 

.0068 

.0068 

.0068 

-0104 

4)105 

.0139 

-0.10 

.0184 

.0184 

.0184 

.0231 

.0232 

#0267 

-0.05 

.0303 

-0303 

.0303 

.0341 

-0342 

#0367 

0.00 

.0500 

.0500 

.0500 

.0500 

#0500 

^5CX3 

0.05 

.0824 

.0824 

.08 24 

.0726 

.0725 

•0675 

0.10 

.1 358 

.1358 

.1355 

.1042 

.1038 

4)903 

0.15 

.2231 

.2218 

.2146 

.1476 

#1467 

.1194 

0.20 

•3540 

.3416 

.3136 

' #2054 

.2035 

.1559 

0.25 

.5137 

.4758 

-4205 

.2792 

.2754 

-2006 

0.30 

.666 2 

.6016 

#5249 

.3683 

.3613 

.2539 

0.35 

.7851 

*1010 

■ .6198 

.4684 

*4566 

.3155 

0.40 

.8665 

#789 5 

.7020 

.5721 

.5543 

.3843 

0.45 

.9183 

.8514 

.7703 

.6709 

*6412 

#4581 

0 *50 

i9504 

.3966 

.8 256 

*1511 

-7294 

-5341 

0.60 

mn 

.9514 

*90'B 

.3329 

-3528 

#6803 

0.70 

#9933 

.9773 

*9419 

.9498 

.9 259 

^016 

i 

0 .80 

.9975 

#9901 

.9723 

#9301 

.9646 

^381 

0.90 

♦9991 

.9956 

.9862 

.9924 

.9836 

#9421 

b 

o 

.9997 

.9931 

.9931 

x9912 

^926 

*9122 


-..ui ^ 
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TABLE 3 0.2. Power of the test T for testing against 

gpr g O ^OS; and r, = 1 / ^ 










d 


6 



... 

Mi 


<p 

‘jr 

3 

10 

12 

14 

16 

-0.50 

-0010 

.0009 

.0009 

.0009 

.0009 

.0009 

'iocog 

-0.40 

.002 2 

.0021 

.0021 

.00 21 

•0021 

.0020 

.0020 

-0.30 

.0050 

.0043 

.0043 

.0047 

.0047 

.0047 

.0046 

•^0 0 20 

.0110 

.0103 

.0107 

.0106 

.0105 

.0105 

.0105 

o 

• 

o 

r 

.0239 

.0236 

.0234 

.0233 

•0233 

-.0232 

<^2^2 ■ 

t 

in 

O 

D 

0 

1 

.0^47 

#0^45 

.0344 

.0343 

.0343 

.0342 

.0342 j 

Q 0 OO 

.0500 

.0500 

.0500 

.0500 

.0500 

.0500 

..0500 i 

I 

0.05 

-0710 

.0716 

.0720 

•0122 

.0723 

.0724 

.0725 

0 .10 

.0992 

.1012 

.1022 

.1029 

•1033 

4037 

.1033 

0 .15 

.1357 

.1403 

.1423 

•1444 

•1454 

.1462 

.1467 1 

0.20 

•IB 03 

.1901 

.1952 

.1935 

•2007 

.2023 

.2035 ' 

{ 

0.25 

.2341 

.2504 

.2597 

.2653 

.2700 

.2732 

.2754 j 

0.30 

.2942 

.3193 

.3343 

.3449 

.3520 

.3573 

♦3613 

0.35 

.359 2 

.3955 

.4173 

-4320 

.4425 

.4505 

.4566 

0 .40 

•4266 

,.4740 

.5026 

.5219 

.5353 

.5463 

.5543 

0.45 

.4942 

.5519 

.5363 

.6093 

•6256 

•6373 

.6472 

0 .50 

.5600 

•6260 

.6644 

• .6394 

.7069 

.7197 

.7:©-i 

0 .60 

.6795 

.7539 

.7937 

.3173 

.6337 

.3440 

.3528 

0 .70 

.7772 

.8437 

.30 23 

.9016 

.9131 

.9207 

•9259 

0 

0 

•3512 

^123 

.9376 

.9501 

.9572 

•9616 

•9646 

0.90 

.9041 

.9517 

*9635 

-^9759 

^793 

-49321 

.9336 

0 

0 

H . 

.9401 

.9745 

.9347 

.0337 

.9907 

.9919 

.99 26 


■ , 
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TABLE 3 #5 .3. Po\^r of the test T for testing agatost H , 

'd’= 12 and r 3 ^=r 2 =1l- 


<P 

6 

8 

10 

12 

14 

16 

-0 #50 

.0009 

.0009 

.0009 

.0009 

.0009 

.0009 

-0 *40 

-00 21 

.00 21 

.0021 

.0020 

•0020 

.0020 

-0 .30 

.0048 

.0047 

.0047 

.0047 

♦0047 

.0046 

~0 .20 

.0107 

.0106 

.0105 

.0105 

.0105 

.0105 

*^0 #1 0 

.0235 

.0233 

.0233 

.0232, 

.0232 

.0232 

-0 .05 

1 

•0344 

.0343 

♦0343 

♦0342 

.0342 

.0342 

0 .00 

.0500 

..0500 

.0500 

.0500 

.0500 

.0500 

m 

O 

• 

o 

.0718 

•0721 

♦0723 

.0724 

.0724 

' .0725 

o 

T-l 

o 

-1018 

.1028 

.1033 

.1036 

.1037 

.1038 

0.15 

•1418 

..1442 

♦1454 

-1461 

.1464 

.1466 

0 .20 

.19 29 

..1981 

.2007 

.2021 

.2028 

.2033 

0.25 

.2550 

.2650 

.2700 

.2727 

.2742 

.2751 

0.30 

*3256 

.3430 

.3520 

.3568 

.3595 

.361 2 

0.35 

.4010 

.428 3 

.4425 

.4504 

.4548 

.4576 

0 .40 

.4771 

.5153 

.5358 

.5472 

.5538 

.5578 

0 .45 

.5502 

♦5990 

.6256 

.6405 

.6492 

.6544 

0.50 

-6179 

.6754 

.7069 

.7247 

.7349 

.7411 

0.60 

.7323 

.7983 

.8337 

♦8535 

.0647 

.8713 

0.70 

.8179 

.8811 

.9131 

♦9 301 

.9394 

.9447 

0.90 

.8787 

.9325 

♦9572 

.9693 

.9756 

.9789 

0 .90 

.9 205 

.9627 

.9798 

•9874 

.9910 

,.99 28 

1 .00 

.9485 

..9797 . 

.9907 

.9950 

.9968 

..9977 
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TABLE 3-5.4. ^cfe,..and^srox W power of t.ha 

tor g = 0.057n.,=n^=io iSd d = 16~ ^ 



-0 *50 
-0 -40 
-0#30 
~0 -20 
- 0.10 
-0 -05 
0-00 
0-05 
0-10 
0.15 
0.20 
0,-25 
0-30 
0-35 
0 -40 
0-45 
0-50 
0-60 
0-70 
0-80 
0.90 
1 -OO ' 


Exact 

Approx. 

Exact 

Approx * 

Exact 

i^pro 

.000^ 

-0001 

.0009 

-0001 

.0017 

.0003 

•0009 

.0006 

-0020 

-0004 

-0035 

-0010 

.0025 

-0027 

.0046 

•0018 

-0071 

-0031 

-0068 

-0099 

.0105 

-0065 

-0139 

-0085 

.0184 

-0304 

.0232 

.0201 

i0267 

.0212 

.0303 

.0497 

-0342 

-0333 

-0367 

-0320 

.0500 

.0777 

»0500 

.0528 

-0500 

.0470 

-08 24 

-1164 

*0725 

-0604 

-0675 

.0672 

-1355 

-1671 

-1038 

-1178 

.0903 

-09 35 

• 2146 

-2302 

-1467 

.1661 

.1194 

.1268 

-31 36 

-3047 

-2035 

-2256 

•1559 

-1675 

-4205 

.3884 

.2754 

.2955 

.2006 

.2158 

.5 249 

.-4777 

.3613 

-3741 

.2539 

-2713 

-6198 

•5681 

.4566 

.4583 

.3155 

-3332 

.7020 

-6550 

•5543 

.5444 

.3843 

.3999 

-7703 

-7344 

.6472 

.6285 

-4581 

-4697 

-8 256 

-8034 

-7294 

.7068 

-5341 

-5405 

.9029 

.9047 

.8528 

-8357 

-6803 

-6761 

.9479 

-9611 

-9259 

-9206 

.8016 

-7916 

*9728 

.9867 

.9646 

.9673 

-8881 

-6784 

.9862 

' .9962 

-9836 

-9885 

-9421 

.9361 

-9931 

.9991 

-9926 

-9966 

.9722 

.9698 




CHAPTER IV 


TESTS OP OTPOTHSSIS FOR LOCATION PARAMETERS 
against TWO-SIDED ALTERNATIVE 

iPJtrodaction and test statistics . 

In this chapter, test statistics based on the LS and 
ML estimators are proposed for testing ^ 2 * 

the two-sided alternative non-mil 

distributions of the proposed test statistics are derived in 
Section 4*2« In Section 4 *3, the LR test statistic for testing 
against H 2 is discussed* An approximation to the critical 
point is investigated In Section 4*4. In Section 4*5, the 
performance of tests is studied and coit^arisons with Tiku's | 

(1981) test and LR test are made* 

Epstein and Tsao (1953) discussed the LR test based 

■ ' I 

on right censored samples* For same problem, Kumar and Patel i 

I 

(l97l), also proposed a test (KP test) based on Uj=| )/0'**‘j, | 

where a* is the pooled estimator of a, given by equation (2 •2*3)* | 

t 

They obtained the mil distribution and have tabulated some | 

\ 

critical points of • Dub^ (1973) and Weinman et (1973) 
derived the power function of the KP test* Weinn^n et also 
compared the KP test with LR test* They used the power function 
expression for the LR test as given by Paulson (l94i)* 

Recmtly, Tilcu (l90l) generalised the KP test for 
left censored samples by pre^osing a test statistic 
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U = 1(X 


( 1 ) 


■ r^+1 




He obtained the null distribution of U • Puirther, he studied a 
Student s t approximation of the null distribution to obtain 
approximated critical points. Khatri (1981) derived its non- 
null distribution. There are some practical limitations for 
evaluating the power function of the test by using the non— null 
pdf given by Khatri. As a particular case of our statistic, the 
power function of Tiku^s statistic (U) is giv^ in Section 402 . 


For the one-sided case we have obtained the statistics 

equationsC 3.1 .3) and (3*1.4) respectively 
Both are equivalent to statistic T given in equation (3*1 *5). 

For testing : 0^ = ^2 ^2 * ^ either 

(4-1.1) Vj. gj, = =Vj^ (say), 

or 

(4.1.2) = |(x‘^j,-x5.^|^)/(7*‘-q2l=|T-qj| = (say), 

where = ^ 1 “^ 2 ' ^2 " nij_-'ni 2 , 

tj^+l 

bj = E (nj^-j+l)"^ and mj,^ = log = 1 , 2 )- 

j=l 

Note that, q 2 = O if = r2/T^2" 

idesntify cases where q^^ = 0 except for the trivial case of 
= ^2 ^1 " ^ 2 * ^^*=230 statistics are of the form 


( 4 •! •S ) 
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where q is a suitable constant indiich is either equal to or 

equal to q 2 * Without loss of generality, we can take q > 0, 
since 

where q* = -q, and the distribution is obtained by Interchanging 
the samples* 

For q = 0, V reduces- to the Tiku^s statistic U* Further, 
for r^ = r 2 = 0 and n^ = n 2 -the test statistics ^MLE' \ 

and the LR test given by Epstein and Tsao (1953) are equivalent. = 

The test procedure is to reject against H 2 / if 
V > Cq^/ where c^^ is determined so that P [v > ^ot a. I 

is the chosen level of significance. 1 

j 

4 .2 • Distribution theory . I 


This section is devoted to deriving the null and non~null 
distributions of the statistic V and its special cases. 

Theorem 4.2*1 . For q > 0, the Cdf of V under is 


(4.2.1) P [v < ciH^] 


GiCc) , O < c < q 


G 2 (c) i q < c < o»/ 


where 

^1 • ^ ““d 

G-(c) = H S (~l)-^( .^)B(r 2 +liff-+j) [{ l+hj^( j)Cg-c)} 
j=0 J 

- {l+h^( j)(q+c)}“‘^3/<%(i>z 



c) 


1 i ^1 

jHq+c)j*’‘V<3hjC j) 

^2 r. 

S' j)(q-c)}'’‘V'3h2( j) 


and the remaining notations are same as in Theorem 3*2#1« 


Proof* Note that V = IT— q|, where T has pdf given in equation 
( 3*2*6)* For the random variable Y = T-q, the cdf is given by 


P [Y < YIHJ = 


F^(y) for y < -q 


^2 r 

where F-(y) = H E (-1 )^ ( .^)B(r* +1, f+j ){l-hQ( j )(q+y)}’*‘V^o( j ^ 

j=0 J J- ^ ^ 

,. r. 

and F^Cy) = 1-H ^E^C-l)-' ( j )B(r 2 +l, f+j ){ l+hj_( j )(q+y)} /dh^( j ) 

Since V = |Yu P [v < ciH^] = P[y< clH^] - P [Y< -clH^ ] • 
Consequently, for q > 0, 


P [V < OIH^] 


0 for c < 0 

P 2 (c)-P 2 (-c) for 0 < c < q 
F 2 ( c)-Pj|_(-c) for q < c < CO, 


On sin^lif ication/ equation (4*2*1) follows* 


The above theorem immediately gives Corollary 4*2*1 
for q = 0, which agrees with the mil distributicai derived by 
Tiku (1901) . 

Corollary 4.2*1 * For q = O, the cdf of Tilcu's test statistic U 
under is 
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for 0 < c < oo« 

For the KP test/ we have q = O and - ^2 - 
cdf is now given in Corollary 4-2.2. This agrees with the cdf 
obtained by Kumar and Patel (1971) • 

, 4 « 2 .2 j » For q 5 = 0/ r^ = r 2 = 0/ we have the cdf of KP 

test statistic as 


( 4 . 2 . 3 ) P [u^ < clH^] = l-{n 2 (l+nj^cy'd)’*‘^+nj^(l+n 2 c/d)"^/(nj^-!-n,) . 

for 0 < c < oD. 

The distribution of the statistic V under H 2 ; ©x / ©2 
can be obtained directly from the distribution of T given in 
Theorem The derivation of the non-null distribution of 

V is similar to that of the null distribution of V obtained in 
Theorem 4.2.1. This distribution is given in the following 
theorem i 

T heorem 4.2.2 . For (p = t. © <3^0/ 'th® of 

V under H 2 is given by 

( 4 . 2 . 4 ) P [V < c|(j) ] = Px(cl (Jj) 

F2(q+cl O < c < q 

F 2 (q+c| ^)-F 2 ^(q-c|f )/ q < c < «>, 
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where F^(»i(j)) and PjC ^Kp) are given in equation (3*3*l)* 

The equation (4 *2 .. 4 ) reduces to equation (4«2#l)^ the 
null cdf of V for q) = 0. For q = 0^ Theor^a 4 •2*2 gives the 
non-null cdf of the Tiku's statistic. 


Corollary 4.2.3 • The non— null cdf of the Tiku^ s statistic is 


n r 

( 4 . 2 . 5 ) P [u < cicp] = Q^(c 3 ^| 0 )-^H S .^)B(r 2 +l/f+j) 

j=0 ^ 


^2 jc 

.exp{h^( j)d<P} Q^{c 3 ^ihL(j)c}/dh^( j)4H 

.B(rj^+l,f+j ) exp{-h 2 (J )d<})} [L^{ G^lh 2 ( j)c} 

~{ l+h 2 ( j )c}"'^]/dh 2 ( j ) for O < c < «>, 
where the notations are same as in Theorem 3 *3 .1 • 

The following lemma can be used to obtain the non— null 
pdf of V as well as U • 

Lemma 4 >2.1 . Let = d<p/c, (p > 0 and Q^(xis) and L^(xis) 
be as defined in equation (3.1.2). The first derivative of 
Q^i *1 .) and L^( .1 .) functions w.r* to c aire given by 


(i) Q^(c 3 ^ 10 ) = e Vc iTd), 


d e 

( ii) Q^£ cj^ = -dh^ S — 

JL~sO 


-c^CX+hjc) 


{ c^d+hj^c)} 


(l-rti^c)^’*'^ il 


+ c^ e 


/c JTd) 


and 
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d -c (l-h^c) £c, (l-h^c)l^ 

[1 . el 1 ] 


H "* ^*1 ( 1 “* 1 ^ c ^ 

C JTd) 


for G ^ 1/h, 


Proof' proof follows at once by writing these functions as 

a sum and then differentiating term by term w»r» to c or by 
writing them as integrals and applying the Leibnitz rule for 


partial derivatives . 


Prom Corollary 4 *2 *3 and Leima 4*2 •!, we get the non-* 
null pdf of the Tilcu's statistic U for <5P > O as 

(4*2*6) fyCu) = H S f )B(r«+l,f+j)e^ 

A —O J ^ 


- E u^Cl+hj^C j)u}^’"‘^‘'Vi! + H s (-l)^’(^^)B(ri+l,f+j) 
isO j=0 ^ 

^2 . r 

* exp {"h2(j)d(^{l+h2U^^j’’^~^ + H S ?) 

3=0 ^ 


B ( r-i +1 , f +j ) -h^ ( j ) d'P — \jL, d . 

S+i ^ [ui{l-h2(J)n}]Vil] 

{l-h2(j)u} i=0 ^ 


for u > O, V. l/h 2 (j) (j = 0,1, •••^r 2 )* 
where u. s: d*^/u i the singularity points, u = 1/h (j) 

mJLi 

the pdf fyCu) is taken as zero* 

It may be noted that, Khatri's (l98l) exprtassion for 
the pdf of U as given in Section l*5f differs slightly from 
the above expr^sion di^e to sewte minor iptegration errors* Also 
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the derivation of the cdf from the pdf given in equation ( 4 * 2 * 6 ) 
is extremly difficult* Consequently, the evaluation of power 
function using Khatri^s approach is not easy, altliough one may 
resort to numerical integration techniques* But equation (4*2*5) 
yields the power function immediately* 

Corollary 4»2*'^ reduces to the following corollary for 

r^ * r2 = 0. 

Corollary 4 *2 *4 » The non— null cdf of the KP test is given by 


(4*2*7) P [ < ci<p] = *^^“’^2 ®3q>(nj^<p)Q^{ nj^c/d}/(n^+n 2 ) 


+ n^L e3?p(-n2<J)) [l^{c^ 1 n2c/d}-f 1402^/^1 

for O < c < »• 

The cdf agrees with the cdf of the KP test derived by 
Weinman at ^* (1973) and Dubey (1973) and given in Section l *5o 

4 *3 * The LR test statistic * 

Let Q = > 0} be the parameter space and 

_ 0 ^ _ 9 ^ a > 0 } be the null hypothesis subset of Q • 

Let 42il 

independent type II doubly censored samples from E(e^,a) (i = 1^2) 
with likelihood function 




(1) J2) 


( 2 ) 


-q 

2 n . !a . t , f i) \ n^i 
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n ."'S . 

1 ^ ^ (i) 




for Gj < y ^ ..( i) 


< ••• < (± = 1 / 2 ) 




'i %+l ^ 

we here discuss the case r^ > 0, r 2 > 0 In detail, Considexably 
simpler expressions hold for other cases* Prcmi equations (2*4,3) 

and (2*4*4), the ML estimators of 03^,82 and a in the parameter 
space 0 are 


i = ^r^ll " ^ log {ni/(nj^-r^)} (i = 1 , 2 ) 


A ^ ^ 

and o = respectively • Hence, the maxiraum value of the 

likelihood function under q is 

/ V ovnf-.,3*) I (nj-r.) 

( 4 . 3 ..I) L(q) = ^ n T- ~: ^ • 

A ^i • ®i' ^i 


Ad*^ 

£7 


In the null hypothesis subset , the ML estimators of 0 and a 

for [see. Section 2*6] are given by 

JL 2 


(4.3.2) ^ -r 

® d* a^+fy~P ^2 


o 








)* 


where 0 ^ is the solution of the equation 

Y/a r^Y 

e o ^ ^ ] / j-l + _ 

p-d a 

o 




d* 0^ + fY-P 


and the notations are same as in Section 2 >6 . 

Prom equation (4.3*2) we immediately get 

{- tt (4^il-8o^> = riY/{a* 5„+(f+ri)Y-Pi 
"o ^ 
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and 


1-exp {- ^ = r2Y/{P-.d* + r2Y} , 

o ^ 


Further, 


n s . 


^ ^ 1 

4“ [ S 


‘1 “i 


Then the maxiinura value of the likelihood function under u is 


A, r 2 n.i 

lCo)) =[ TT ^ 


iJl ^i^^^^d*a^+(ri+f)y-.P - "p-d'^^^+r2Y 


r^Y 


.]^1 [ 


Oi 

] 


^ A 

d*a„+fY~P 


•Hta 




d*0^+(rj_+f)Y-P 
Now/’ the likelihood ratio is given by 

A A 

X = L((o)/L(q) 


( 4 .3 .3) = const . fi-!-** Y’'l'”'2(d*5„+fY-P)* e^(-P/5„) 


_______ 

{ d*0^+( r^+f )y-.P} ^l"^^ £ P-c^0^+r2Y} 


for Y > O, 

2 n^ -(n^~r^) 

where Const* = e>p(d*)* M n^ ^^i^^i^ * 

i=l 

For Y = <0, X is obtained by replacing r^,n 2 /rj_,r 2 

r 2 tx r tX 

by n 2 #n^^ ^2'^1 equation (4 *3 •a)* The LR test 

then rejects if ^ 5. where is chos^ so that 

P [X < X^iHq] = a. 
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4#4» Critical points of the te<!-i-c t.r \t 

-cne tests V.,^Vy U and X . 

crxtxcal point o^ of V is obtained by solving the equation 
(4-4.1) P[v>c^|Hj,a, 

whejTS d is tne chosen level n-F o4/wi-p» 

significance* Fix)m equation 

(4»2*l) it is clear that, c is ei-hh^T' +-v,f» i 

* a either the solution of 

r^ 


(4.4.2) H [fi+h^CDCq.^^,, 


-d 


- {l+h^(j)(q4<^)r<^]/dj,^(j) ^ 

(case in which < q) or the solution of 
h r 

( 4 . 4 . 3 ) H ^l)B(j-2+l^f+j) {l+h^(j)(q+c^)}“‘^/dh^(j) 

""2 r 

jio ^ ^ ^®^^l+l'f+JHl-h2(j)(q~CQ^)}‘''V<ah2( j)] = a 

(case in which c^j. > q) according as > 1-a, or P^. < i«a 
respectively, where 

Pg == P [v < qlH ] = H E (-l)^(^^)B(r,+l,f+j) 

j =0 ^ 


. [l-{l+2h^(j)qr'^/dhj(j). 

The value of c^ can thus be obtained by first evaluating P and 
then applying Newtonr-Raphson method to the relevant equation* 


Unlike the one~sided test statistic disciussed in Chapter III, 



88 


here it ivS not possible to get a compact expression for even 
for = r 2 = 0 case* 

Now, the critical points tests 

V and U arc nothing but the' solution of the equatic»i C4*4*l) 

1 2 

with q = ^ respectively. It is clear from the 

equation (4-2.1)/ that the critical point c^^^ of the test U is 
simply the solution of the equation (4*4.3) with q = 0» Some 
of these critical points are tabulated in Table 4*4.1 for 
studying tho performance of these satistics. 

Since the LR test statistic K given in equation (4*3.3) 


is very complicated/ we have not attenpted to derive the 

distribution of X . However, we have obtained simulated critical 

points Xq, based on 10,000 iterations, for studying the 

porfonnance of the statistic. These are tabulated in Table 4,4.1* 

For r =: r. = 0/ the tabulated critical points are exact. 

1 2 

In section 4.5, we have oonpared the perfomance of 
these statistics. It is very difficult to ccropare the power 
functions of the test statistics. But Tables 4.5.1, 4.5.2 and 
some other numerical calculations show that, to between Vj, and 

n+-ne difference of the power values* However, 
V„ there is very little Cittereric. 

^ t- ir H^ncG we are suggesting a method 

V, is slightly better than V^. Hence we ar gg 

for obtaining an approrfmate value c„ for , a 

also be taxen as an initial value for solving e^bton (4.4.1) 

for An Identical approach can he used fox o^ an „ 
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4*4.2. Approximated critical point_cJ of . m Section 3 - 4 , 
we studied an approximate null distribution of T. In the 
present case, let = T-q^ so that * 2’hen frtxn 

equation (3 #4 •7) we have 

(4*4.4) E(Y^) = qj_/(d-.l)^Var(Yj^)=d^ { (d-l)A4qJ/[(d-l)^(d-2)] 

with = B. Now, it is clear from Section 3*4 that-Y^/VA has a 
Student^ s t distribution with d DF* We improve it slightly by 
considering two constants a(> O) and b such that mean and variance 
of the random variable a(Y^4-b) is equal to the exact mean and 
variance of a Student's t random variable with d DF. That is^ 

E [a(Y^+b) ] = 0 and Var [a(Y^+b) ] = (^(d-2) * 

Then on using equation (4.4*4) we get 

(4-4*5) a = (d-l)/[d{(d-l)A4<j^}V2^ b = -q^/(d-l) . 

Since P [V^ > = 1 - P[tY^| < 

= 1-P [a(~C|^^^+b)<a(Yj^+b)<a(c^^^+b) ] 
= 2P [T^ > t^] , 

( 4 .4 *6 ) ^ ^ 

where t^ is the (a/2)th percentile point of Studcaat's t 
distribution with d DF. Some of the exact critical points c^. 
[solution of the equation (4.4.1) 1 and its approximated value 
c* [given in the equation (4.4*6)] are tabulated In Table 4.4*2 
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for different combinations of sanple sizes and censoring patterns 
Note that, for very small values of a and , (rj/nj^-r/n^) , > o.i, 
the approximation is not satisfacto.r^^. 


? . Similar to the approach discussed in Section 3 •4, 
we also studied a normal approximation for the critical points 
of But these are not as good as the Student's t approximated 

values for most of the cases* Hence., these values are not 
tabulated* 


4 oB • Power of the tests ^ ^ ’ 


The power of the test V for testing : 9^ ” ®2 
: 0^/02 is given by 


(4*5 *1) 


P [V > = 


_l-Pj(c„|CP) 


for 0 > 0 


for (p.< O, 


where c^ is the exact critical point/ P^(c^ip) is given by 
equation (4*2*4) and P 2 (cQ,t<p) is obtained by interchanging 
n^/n 2 /r^/r 2 by "‘^ 2 * 2*^1 equation (4*2*4) and evaluating it 
for I <P1 « 


The power functions of ^ obtained from 

equation (4 *5*1) by replacing (q/CQ^) ty (q27e^^^) 

and (O/C^^^) respectively* Some of these power values are 

CCr 

tabulated in Tables 4*5*1 and 4*5*2* 


The power of the LR. test statistic X is given by 
P [x < C- i qj T* As wo had mentioned ih Section 4*4/ the distribution 
theory of X is very ccmplicated* Hence we obtained the power 
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by applying Monte-Carlo technique# using 1000 iterations* 
However^ for r^ = r^ = 0# the exact power was obtained by using 
the pox-zer function expression as given by Paulson (I94l) [see# 
also Weinman ^ 1973]. 


For the following sets of saitple sizes# censoring patterns# 
and various values of <P we studied the conparative performance 
of all the four test statistics for a = 0*05* 


(i) nj^slO, n^^B# d=14#(r^/r2)=(0#o) and (r^/r 2 )*(l/l) in 
Table 4 *5 *1 • 

(ii) n^=n2=15# rj^=l,r2»3, d = 18 and d = 24 in Table 4*5.2* 

To high light the shape of the power functions# three such curares 


for n, 


1 . 


nj = 15, = 1 , rj = 3 and a = 24 are drawn to 


Figure 4*5*1* 

Since the difference in power values of and is very 
small, we have not drawn the power function of in Figure 4.5.'l. 
In all cases studied# it was observed that the -average^ power 
of was always greater than the "average power of 2 * 
the average waa taken of the two values corresponding to t<P and 


(P 

A normal approximation similar to thti one studied in 
section 3*5# for power function is also considered* Since 
Y = T-q, , the power of the test is giva:i by 

P [Vj_ > ^ <P ] = ^ 
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(4.5.2) -'■[Wiio] +P[„^<0] , 

inhere W. = -x^ 2) (l)^•^^■ M"* ( 'i'\ 

1 r.+l V+i )a and W -xi^; (g 


2-"-r^+r^r2+r'‘Jl*-^a 


similar to the reeulte obtained In Section S.5, we see that 
(1 = 1,2) has an asymptotic normal dlstrltajtlon with 

S(Wj^) = (IJI _ a^l*)o s Pj^o, 

Var(W^) = {A+(,^+o(l))%,„2 . 

eCWj) =((11 + c^^ha S pjo 

and VariWj) = { A+Cq^-c^l' ) £ g/^. 

Hence from ecjuation (4.5*2), we immediately get 


(4.5.3) p [v^ > da^’iK] = »(Pi/VgJ) + jC-p/ifij). 

The exact power values obtained from equation (4*5 »l) and 
its approximated values given in equation (4 .5 *3) aro tabulated 
in Table 4 *5 *3 for a = 0»05 and sele(rted values of 
and d. 


Using the power function as a base, the following conclusions 
may bo drawn. 

(a) All the three tests ^ biasedi However, the 

extent of the bias is different for different tests. 

( b) There is a considerable loss of power for all tests ev^ 
with mild censoring on the leift as is evident from 
Table 4.5 But the loss of power is negligible with 
variations in censoring on the tight* 



(c) 

(d) 

( e) 

( f ) 


The test statistic U proposed by Tiku (1931) is more 
biased than and In general the statistic U shows 

poor performance as I ( Vn^^-r/n^) l increases. Further, 
the bias of U seems to increase with increase in d. 

Table 4.5-1 shows that the test statistic V 2 is relatively 
more biased than • 

We do not have sufficient evidence to conclude that, the 
LR test statistic \ is unbiased for r^ > 0 and/or r 2 > O, 
since some simulated values are less than a as in 
Tabic 4-5.1. However, for r^ = r 2 = 0, Dubey (1973) 
established that, the LR test statistic is unbiased [see, 
also Khatri 1974]. 

Table 4 . 5 . 1 , Table 4*5.2 and Figure 4*5»1 show that there 

is very little difference in the power values of X and 

• Since the statistic X given in equation (4-3-4) is 

very complicated while the statistic given in equation 

(4*1 #3) is considerably simple, we strongly recommend the 

use of test statistic in such situations « 

The normal approximation for the power function is fairly 

good if 1 ( r^/n 

values of n^ and n 2 . However fdr large values of D, the 

approximation is not that good even with moderately large 

values of nj, and n 2 . This is well illustrated in Table 

4 -5 -3/ for the case n^^ = 15#^ n 2 = 25 and d =s 28 • For 

(r ,r ) = ( 0 , 4 ), the maximum differtaice (among all tabulated 
12 

values) between approximate exact values is 0-0274 for 
(p sz 0»30ff whereas for the iiaximura differciicc 

is 0-0566 for » -^40* 


I^~r^n 2 ) I = D (say) is small, even for small 



TABLE 4.4.1. Critical points and of the 

- — Ill — I P - (X Ou Qu 

tests ^]_/V 2 ,U and X respectively for a = 0.05* 


"l 

^2 

^1 

^2 

d 

.s . 

J 2) 

% 

4 '^ 


10 

8 

0 

0 

14 

0 .3769 

0.38 25 

0.3825 

0.2386 

10 

8 

1 

1 

14 

0.5589 

,.1 

0.5654 

0.5760 

0.0861* 

10 

10 

1 

1 

16 

0.48 27 

0.48 27 

0.48 27 

0.1104* 

15 

15 

1 


18 

0 .4049 

0.4069 

0 .5 239 

0 .1576* 

15 

15 

1 


24 

0.39 29 

0.3947 

0.5103 

0.2063* 

i 

20 

15 

4 


26 

0.4194 

0.4223 

0.4223 

0 .0589* 


^Simulated values based on 10000 sartples for each saitple size* 



TABLE 4#4^2^ Exact critical point c^*^ in tC3p row and its approximated value 
given in ecmation (4«4«6) in bottom row for a = 0«05» ^ 
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table 4»5»1» Exact power of the tests and exact or 

simulated power of X for a = OiOB, 
n 2=6 and d * 14. 





Vi 

X 



-1 

• 00 

-0 

•80 

-0 

-60 

-0 

.50 

-^0 

.40 

~0 

• 30 

~o 

•25 

-0 

• 20 

-0 

•15 

-0 

•10 

-0 

.05 

0 

.00 

0 

.05 

0 

•10 

0 

•IS 

0 

• 20 

0 

• 25 

0 

• 30 

0 

• 35 

0 

•40 

0 

.50 

0 

•60 

0 

•80 

1 

•00 


-9977 
-98 33 
•89 27 
•7619 
•55 24 
•3194 

• 2249 
•1545 
•1064 
•0755 
.0575 
•0500 

•05 23 
•066 1 
.0955 
•147 2 

• 2216 
•3358 

•460 2 
•5341 
•78 20 
•3964 
.9787 
•9957 


•998 2 
.9871 
•913 2 
•798 2 
•5987 
.357 2 

• 25 36 
•1743 
•119 2 
.08 3 2 
.0612 
•0500 

.0483 
.0567 
.078 2 
•1180 
.18 27 

• 2746 

.3879 
.5090 
.7 243 
.8631 
.9710 
.9941 


.9971 

.9795 

•8739 

.7308 

.5158 

• 2915 

• 204 2 
•1403 
.,0972 
•0100 
.0549 
.0500 

.0553 
.07 29 
•1079 
•1679 

• 2586 
.3761 

•505 2 
•6XJQ 
.8121 
.9125 
•98 21 
.9964 


.9388 

•8089 

.5580 

•4099 

• 2759 
•17 29 

•1346 

♦1046 

.0818 

.0655 

.0550 

.0500 

.0504 

.0566 

.0696 

.0905 

*1209 

.1620 

• 2147 

• 273 2 
.435 
.5856 
•8 258 
.9403 


♦9490 

•8323 

•5907 

•4404 

.3001 

.1894 

•1476 

•1143 

•0883 

.0701 

.0573 

.0500 

.0479 

.0512 

.0605 

•0169 

•1017 

•1362 

.1315 
.2376 
.3765 
•5312 
.739 2 
.9 242 


•9561 
•8498 
.6 169 
•4657 
.3205 
• 20^6 

•1588 

•1229 

.0950 

.0741 

.0594 

.0500 

.0457 
.0464 
.05 25 
•064S 
.0843 
•1125 

.1505 

•1938 

.3239 

.4723 

.7443 

.9023 


.923 

.779 

.535 

•403 

•254 

.155 

.124 

.083 

.078 

•051 

•052 

.048 

-^41 

.049 

.074 

•016 

•106 

.173 

.221 

.238 

.455 

•613 

.835 

.941 


^Simulated power based on 1000 samples for each sanple size- 
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table 4*5*2* Exact power of the tests and simulated 

p^wer of the test \ for a = 0*05/ n^=n^=l5, 
rj^=l and r2=3» 


d=18 

d=24 




<P 


• 1 -00 
■ 0.90 
■0^30 
. 0#70 
-0 #60 
-0 #50 

.0 -40 
-0-30 
'0 -20 
-0 -16 
-0 -10 
0 -00 

0 -05 
0-10 

0 -15 
0 -20 
0 -30 
0 #40 

0 .50 
0 #60 
0 .70 
0 .30 
0 .90 

1 .00 


.9960 
.9857 
.963 2 
.9128 
.3177 
.6700 

.4863 

.3073 

.1703 

.1221 

-0864 

.0500 

.0474 
.0563 
.0814 
.126 2 
• 233 2 
.5143 

.7 206 
.8597 
.9363 
.973 2 
.939 3 
.9960 


.9964 
#9871 
#9650 
#9163 
#8 234 
#6773 

.4938 
#3134 
.1741 
.1249 
#088 3 
.0500 

.0463 

.0541 

.0764 

•1180 

.2725 

.4955 

.7054 

.3503 

.9314 

.9709 

#9304 

.9959 


.9970 
.99 21 
.978 2 
.9435 
-8702 
.7421 

.5641 
.37 24 
.2130 
-1536 
.1077 
.0500 

.0338 

.0237 

.0185 

.0179 

.0331 

.0343 

#19 27 
#3593 
.5507 
.7 213 
.346 2 
.9 229 


.997 

.987 

.950 

#890 

.768 

.573 

#386 
.^0 
#112 
.103 
-072 
.05 3 

.058 

.069 

.131 

#191 

.374 

#618 

.795 

•904 

#953 

.937 

.995 

.995 


#9973 
#9919 
#975 2 
.9336 
.8461 
#6978 

.5045 

.3143 

-1710 

.1216 

.0864 

.0500 

.0435 
.0596 
.0869 
.1353 
.3097 
#5 430 

#7451 
.375 3 
.9445 
-9771 
.9910 
.9966 


#9970 
.99 24 
.9765 
.9367 
.6517 
#7058 

.5130 

.3212 

.1754 

.1248 

.CB84 

.0500 

.0472 

#0566 
.<B15 
.1269 
.2931 
#5 244 

.7310 
^670 
.9404 
.975 2 
.9902 
#9963 


-9967 

.9949 

#9861 

.9596 

.3966 

.7741 

.5910 
-3877 
• 2189 
.1567 
.109 2 
#0500 

.0336 
.0233 
.0180 
.017 2 
#0319 
.0337 

#1973 
.3763 
•5193 
.75 24 
10707 
^333 


.994 

.939 

.973 

.896 

.779 

.588 

.396 

-219 

.113 

-105 

.064 

.056 

•066 

•093 

.173 

.223 

.418 

-669 

^25 
.9 26 
.967 
.938 
.99 3 
.999 


*slittalated power based on 1000 sainples for eaoh sanple sloe. 



TABLE 4-5*3# 


app^xiinate power values of the test 
statistic V. for a = 0-05- ' ‘ “ 


("1=11 2=10, dsl6) 


(fi 

ri=r2=l 

r^=o. 

r2=l 


Exact 

Approx. 

Exact 

Approx. 

-1 .00 

.9736 

#9832 

0938 

0956 

-0 *90 

.9551 

.9568 

-9841 

0851 

~0 #30 

•9098 

.9038 

.9606 

.9577 

-0 .70 

0 293 

0140 

.909 2 

0990 

—0 #60 

.7045 

«6351 

0109 

#7960 

~0 .50 

.5 381 

#5 284 

0555 

.6475 

-0 .40 

.3593 

#367 2 

•4624 

.4721 

"0 .30 

#2110 

#2271 

.2313 

#3021 

-0.25 

•1560 

•1709 

#2094 

#2291 

-0.20 

#1143 

.1253 

#15 26 

-1675 

-0.15 

•08 44 

.0905 

.110 2 

#1186 

—0 .10 

.0647 

.0663 

.0801 

00 25 

-0 .05 

.0536 

#05 21 

0605 

0590 

0 .00 

0500 

0474 

0500 

.0473 

0 05 

.05 36 

#05 21 

1 0481 

0490 

0.10 

.0647 

.0663 

0558 

.0634 

■ 0 -15 

0844 

.0905 

.0762 

09 25 

0 .20 

.1143 

•1253 

#1145 

#1377 

0.25 

#1560 

#1709 

#1774 

.1996 

0 .30 

#2110 

#2271 

.2675 

.2776 

0 .40 

#3598 

#367 2 

#4967 

#4683 

0 .50 

#5 381 

.5 284 

#7066 

#667 2 

0 .60 

.7045 

.6851 

0456 

0 275 

0-70 

0 298 

0140 

• 0 250 

#9 273 

0 -30 

.9098 

.9038 

0649 

0755 

0 .90 

#9551 

.9568 

.9840 

09 34 

1 .00 

.9786 

.9332 

0929 

0936 


xou 


TABLE 4 *5 f 3 • Contd • 


(nj^=3.5,n^=25,d=28) 


<p 

^l=4/r2=0 

r^sO, 

jC^=4 


Exact 

Approx- 

Exact 

Approx- 

.60 

.3719 

^717 

-996 3 

•9913 

•^u *5 5 

•3199 

.3041 

.9739 

.9735 

-0.50 

•75 22 

.7185 

•9546 

.95 27 

-0.45 

*6673 

.6181 

-9141 . 

.9065 

-0.40 

•5657 

.5091 

.8443 

-8332 

-0 .35 

.4511 

.3996 

.7412 

.7306 

-0 .30 

•3314 

.2979 

.6075 

.6034 

-0.25 

-2185 

-2108 

.4588 

.4641 

-0 .20 

•1264 

*1424 

-3179 

.3291 

-0-15 

.066 2 

-0941 

.2032 

.2136 

-0.10 

.0399 

.0648 

.1222 

•1269 

-0 -05 

♦0 384 

.05 26 

.0729 

.0715 

0 .00 

.0500 

.055 2 

.0500 

.0464 

0 .05 

.0701 

-0714 

.0507 

.05 21 

0 *10 

.0984 

.1003 

.0796 

.0931 

0 .15 

• 1362 

.1418 

.1533 

.1758 

0-20 

-1847 

.1957 

.2913 

.3012 

0.25 

• 2448 

• 2614 

•4722 

-4577 

0 .30 

.3166 

.3371 

.6490 

-6 216 

0 .35 

.3985 

.4202 

.7879 

*7659 

0 .40 

.4875 

.5070 

-8814 

.8732 

0 .45 

.579 2 

*5935 1 

.9377 

-9403 

0 .50 

.6684 

•6757 

.9639 

.9753 

0 .55 

•7499 

.7500 

.9851 

.9916 

0 .60 

,f9199 

.8142 i 

.9931 

.9975 





CHAPTER V 


GENERALIZED STATISTICS FOR K RIGHT CHiFSORED SAMPLES 
Introduction * 

So far we considered the case of testing the equaJLity 

of location parameters of two exponential distrdhutions • In 

this chapter, we propose two test statistics for testing 

~ ®2 ~ *** ” ®K ~ ® against the alternative hypotheses 

H. ; 0. > max (0 .) and H„ : at least one 0^ (i = 1,2, i 

2<j<K -J ^ ^ 

is different from 0, based on K(> 3 ) independent right censored 
samples* The necessary distribution theory of the proposed test 
statistics is discussed# Some critical points and power values 
are tabulated# Finally, we compare the performance of these 
statistics with the test statistics proposed by Khatri (1974) 
and Singh (1983)# 

Let X$ (i = 1,2,-.#,K) be K Independent 

1 ik “i~®i 

sartplos frQm.E(0^,a), where n^-s^ > 1# For simplicity of notations, 
let X^ s= X^^^ (i = 1 / 2 , and - min(X^,X 2 , •• '/X^)/ that 

is, X^ is the minimum of the iWi sample and j_) is the minimum 
of all the observations# 

Khatri (1974) derived the LR tost for testing against 
H 2 # It is given by 

K ^ „ 

* E nj^(Xj.-x^ j^))/da , 


( 5 . 3 . 1 ) 



where 


(5-1.2) 


K 

2 : { I 

i=l j*il 






),d = 

1 


K 

S (n£-S£-l) 
i=l 


are same as In earlier chapters- He obtained the power function 
of the LR test and showed that 


00 

(5 -1-3) P > c|H^] = {B(d,K-l)}“^ / -/•'2(i+yrK~d+l 

e 

From equation (5-1-3), it is easy to show that dUj^/(K-l) has 
an F~distrifoution with 2 (K~l) and 2 d DF, which is denoted by 
^ 2 (K— 1 ) 2 d* also derived two union intersection test 

statistics from two different view points- These .are given by 

(5*1.4) U2 = [ma^nj^(x£-x^^ ))/n2(x2-X(^)), ...,nj,(3^~XQ))}]/do'' 

and 

(5-1-5) U3 = [max{n 2 (x 2 - 3 ^)» •-•y%:(xj^X3_)/n^(x^-X2)/ 

nj^(x^-xj(;)}]/<30'**'- 

In all three cases, the test procedure is to reject if 
Xj, > otherwise accept H , where the constant c^ (i— 1/2,3) 

is determined by solving the equation P [Uj_ > c^ IH^J = a- 
Khatri (1974) also obtained the null distributions of U 2 and U^, 
and their power functions under the assuitptions 6 ^ > 62 ^ 
which can be achieved by renaming the populations. He has 
provided some critical points of ^2 ^3* ^ pow^r 

comparison studies for H 2 and U 3 have been done- 
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Singh (1983) also dlsoussaJ the LR tost procedvire for 
testing against The test prooea.re Is egulvalent to 


rejecting If ^U) . 


"4 -- ^Sj^"i(Xi-X(jj)/{(K-i)a»} and P [u^ > o^‘‘’|H^] = a. 

He has Shown that, the distribution of 8^ is , but he 

did not study the power function of U^. Note that,U^ = dU^/(K-i) 
and the power function of has been already given by 
Khatri (1974)* 


Although there are number of tests for testing 
against H 2 based on some theoretical considerations as described 
above, x^e propose another test based on 


T 2 = f max Xj ~ min XA/a*m 

This is mainly for cortparing the performance of various tests# 
Further, a generalization of T '2 left censoring is easier 
than that of the other tests* Note that, if all samples are of 
equal sizes, then T 2 is equivalent to U 2 » 

5 *2* Test statistics and their null distrlbutians * 

Similar to the one-sided test statistics of the two- 
sample case, we propose 

(5.2.1) T^ = {Xj - rain(X 2 ,X 3 , 

for testing against Similarly, for tasting against 

H 2 # we ppoposG ' ' 
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(5*2*2) T, = { max X. - min x.}/a*. 

l<i« ^ l<i^ 1 

The test procedure is to reject if > c^, where 

^ — ^i'^o^ ~ ^ “ 1#2)/ and a is the chosen level of 

significance* For obtaining the critical points c^^ and C 2 / 

need the distributions of T. and T„ under H . These are 

1 ^ o 

obtained from the following lemmas ; 

Lemma 5 ;»2'»1 * If • * */Wjr are K independent random variables 

with pdf of given by 


( 5 . 2 . 3 ) %/w) = n^ ejqp (~nj^w)y w> 0 (i = 1, 2, •*./ K), 


then the pdf of Z = W^~min(W 2 / • * 'vJ^) is 

”n^(N~n^) exp {(N-n^)z}/N, 

n^(N-n^) exp {~n^z}/N, 


(5*2 . 4 ) f(z)=s 


z < 0 

z > O, 


where N = nj^+n2+...+i^* 

P roof . Let Y = min (W2/W24r » -./Wj^) . Then the cdf of Y is given 
by 

K ^ 

P [y < y] = 1 - n exp(-njj^y) = 1 -exp (-N 2 y), where N 2 = S n^^- 

i =2 

Consequently/ the pdf of Y is 


f(y) = N 2 e^^* (~N 2 y)/ y > 0. 

The result now follows/ on applying Leraia 3*2.1. 

Lenma 5 . 2 . 2 . Let ( 1 = 1. 2, . .^.K) be K Independsit ranaom 
variates with pdf given l:y agnation (5V2.3) . Then the pdf of 



106 


Z = max Wj ~ min Wj is 
l<i<K l<i<K ^ 


1 K K 

(5 *2*5) g(z) = i 2 E 


~n.2 K 

n^n.e H 


i=l j=l/j/^i ^ ^ hssl^tV^i, j 


{ l-expC-n. z)} , z >0, 


where N = n^+no+«* 


24- * # I 


Proof* Since (is 1,2/ •!>'./K) are K independent random 
variates with pdf f^ (w), hence, the cdf of Z is [see, David 1981, 
p . 26 J given by 


K ■» K 

(5.2*6) G(z) = E ; f^ (w) n {P^ (w+z)~F^ (w)}dw, 

i=l -«> i j=l/J>i j j 


where for w > 0, (w) = l-ej^C-njw) is the cdf of W j • 

Substituting for (w) and F^ (w), we have 


K oo -n.w K -n.w 

G(z) = E / n.c • 71 {e -e 

i=l o j=l/j/i 


-n.w -n.(w+z) 
e ^ -e ^ } dw 


K 0 ° K -n.z 

= E ; f\ {1-e ^ } dw 

i=l o j=l/jA 


K n. K 


= E 


-n .z 

{1-e ^ J 


±=1 


(5 * 2 *7 ) 


K 

E G^(z)/ 
i=l 


n. K 

where G^Cz) = IT {l-estpC-nj 


-e3{p(— H jz)} Ci = 1/2 /»-»xK) 


Now, log G.(z) = Const. + _ 


log ♦ 



Dlfforvmtiatlng w.r. to z, on both the sides of the above 
equation, we have 


SG . ( z) 


gi(2) X 

Sj^T'z) ~ ^ ®^(~njz)/{ i~exp(-njz)}] , where 


Hence, gj^(z) = -|j- S n. exp(-n.z) fl { l-es^pC-nt^z)} 

From equation (5*2*7), the pdf of Z is given by 

K 

g(z) = 2 g^Cz) 

i=l * 

which gives the required equation (5.2*6)* 

Theorem 5*2*1 * Under the null hypothesis the statistic T^ 
defined in equation (5*2*l), has the following cdf : 


(5 *2.8) P < clHj^ ] = 


-d. 


l~(N-nj^)c/d} /N, c < O 


-d. 


l-(N-n^)(l+n^c/d) /N, c > O, 


K K 

where N = S n. and d = E (n£~s^~l)* 
i=l i=l 

Proof* The random variable = (-X|^-0j_)/cr (i = 1, 2, •••/X) has 
the pdf given in equation (5*2*^)* Consequ^tly, Z = W^~^mi^ 
follows the distribution as given in equation (5*2*4). Then 
under T^ = d^W, where W = dff*/a. Now, applying Theorom 3*2*1 

we get the required equation (5*2*8)* 

For general K and n^^, the distribution theory of the 
statistic T 2 is very complicated, although it is possible uo 




follow the same approach* We therefore consider the siitplifylng 
assumption K = 3 in Theorem 5*2#2* 

K = 3< the null pdf of T 2 is given by 
3 2 

(5*2*9) fCt^) = ^ S nj_(N-nj_){(l+nj_t 2 /d)~^''^-(l+ S n^tVd)”^^} 

i=l h=l,iv6i ^ 

for t 2 > 0 , 

where N = + 02 +^ 3 * 

Proof * Note that, = (X^-e^)/cr (i = 1,2,3) has the pdf given 
in equation (5*2*3)* Prom Lenina 5*2*2, for K =s 3, the pdf of 
Z = max( W^, W 2 , W^) ~ min (Wj^,W 2 ,W 2 ) can be rewritten as 

j 

3 3 

g( 2 ) =« 2 nj(N~-nj ){ exp(-n.z)-exp(~ S n^^z)}/ z > 0 * 

^ i=l ^ ^ ^ • h=l,Vi ^ 

Since T 2 = d^W, where W = d 0 */d, we obtain the desired pdf of 
T 2 by proceeding on lines similar to that of the proof of 
Theorem 3 -.2*1* 

If n^^ = n 2 = -•* = nj^ = n, then the statistic T 2 given 
in equation ( 5 * 2 * 2 ) is equivalent to T 3 , where 

(5.2*10) T 3 = nT 2 S dU 2 * 

For equal sample sizes, we therefore use T 3 * The null 
distributioh of T 3 can be obtained as before* It can also be 
derived from the cdf of U 2 as given by Khatri (1974)* Note 
that,Khatri has used v/s for V 2 P ^ using 

Khatri^ s express ion we have 



( 5 . 2 - 11 ) 


luy 


p [T3 < CIH^] = p [U2 < c/dlH^] 

K-2 

= 1-2 (-l)^(f“J){l+( j+l)c/d}"'^, c > O 

j r£) J A — 

For studying the performance of tljese statistics^ the 
non— null distributions are derived in the next section. For 
K !>• .3, tne expressions for po\ver functions are very complicated. 
Honce^ wo have only considered the case K = 3 in detail. 

Without loss of generality, we have assumed > 62 > ©3/ since 
this can be achieved by relabeling the sample. 

5 .3 . ■ Non- null distributions of the statistics . 

The non-null distribution of is derived by using 
the following lemma : 

Lemma 5 . 3 * 1 . Let and W3 be independent random variates 

with pdf of Wj^ given by 


(5.3*1) f^ (w) = n^ exp {-n^(w- cp^)} , w > (i = l/ 2 i 3 )/ 
where > ^2 - ‘*’3' ^ ^ Wj,-mln(W2, W3) is 


(5.3.2) f(z) 


(n2+n3)bi exp{ (n2+n3)z}, ~«> < z < 0.2 

n^b2 C2xp(-nj^z)+n3b3 exp(n3z), ^2 z < 


L (~nj_z), a3 < z < 

where * 2/?)/ H = ^ . 

} 

^2 “ ^^^2 (nj^CC2-n3a3-m3<i2)/{ ' 

b3 = UjL Gxp (-n3a3 )/( 0^^113) and b^ = ng exp(n^a3)/(iaj_+H3). 


liu 


Proof. The cdf of Y = min (W2,W3) is 

[W2 < y]Ml ~ P [W3 < y]}. 
From this, the pdf of Y is given by 

Hj' exp{-n3(y-<P3)} , % 1 Y < ^2 

(n2+n3)exp{-.(n2+n3)y+(n2<P2+ii3‘P3^J / (p^ < y < 


f(y) 


OOm 


Now/ from the jpdf of Y and and making a transformation 
2 = ~ y and w = W; wo get the marginal pdf of Z as 


(5.3.3) f(z) 


; Pi(z,y)dy/ 

/A — 


—00 < z < a. 


<Pj^— z 




S p. (z,y)dy +/ PoCz,y)dy, < z < a, 

<l>2 ^ 

^ PjL^2,y)dy + ^ < z < 


<P 


3 3 

whoxe p^(z,y) =s nj^(n2+n3) eiq) ■*■ .^^^i^i "■ 


and P2^2;,y) = ^1^3 {-n^z + n3f3-(nj+n3)y} • 

Equation (5 . 3 . 3 ) now gives the requ.ired equation ( 5 .^- 2 ) on 
s impl if ication . 

Theorem 5. 3 . 1 . For > ©2 ^ non-null dlstrlbutlOTi of 

T- is given by 

P [^1 < 
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\ ^ < o 

QgCojio) + bj^ L^CCji (bj+njjc/a) 

( 5 *3 *4 ) = 

~b 3 [L^(c 2 ln 3 c/d) ~ b^(c 3 in 3 c:/d)] 

-[b^QjjC C2ln^c/d}+b4Q^{c3l3^^c/d}] , c ^ 0 , 

w> 

where Cj = daj/c, = (e^-e^Va (j _ 2 ,^), Q^(»i«)/ 

are given in equation (3.1.2) and ( i = 1 , 2 ^ 3 , 4 ) are given in 
Lemma 5 *3 .1 . 

~£22i£* bet = X^/a (i = 1,2,3). Then are independoit 

random variates with pdf given as in equation ( 5 # 3 *l) with 
<Pi = 0^/d. Tho pdf of Z = W^-min( W 2 / W 3 ) is given in Lemma 5.3.1* 
Clearly, can be written as d^W, where W = dc^/o • Now, using 

similar arguments as in Theorem 3*3.1, the cdf of T^^ upto the 
point c is the integral of the joint density of (T^,w) over the 
shaded region shown in Figure 5*3. I* 



FIGURE 5 . 3*1. Region showing the cdf of T^ upto the 

point c. 



AJL^ 


Consoqnently, 

^ [^1 1 P-j^citj^) dw = (say) for c < O 

and 

P [^1 < ciy = / (; P3dt^)dw+/ (; p.dt. )dw 

o 0 0-^1 

oo < 3 a^w C3 ^ 

-j- / (/ Pj^dt^dw-f. ; (; p„dt. )dw 

° ' ^2 ^ ^ 

00 e 

+ / (/ p^dt-)dw+/ (/ p-dt. )dw 

G3 ^ C3 dttg/W “ ^ 

= d’2'^3'^4'^5’^6'^7 (say) for c ^ O, 
whore = (n2+n3)b^ exp [~w{ l-(n2+n3)t^/d}] w'Vd! , 

P2 » ^^1^2 f’*w( l+5^3^‘t3_/d)} w/d! + n3b3 expC-wCl+n^t^/d)} w / d ? .. 

and 

P3 = nj^(b2+b^) exp {-w(l+nj^t^/d)} w^/dl- 

Simplificat ions for J^/s (j = 1/2, »'•»/?) are siinllar to that 
of of Theorem Their simplifxod foxms are as 

follows : 

J31 = bi { I~(n2+n3)c/d}”'^ 

J 2 = b^ 

J 3 =s bj^ E ^d^ *^ 2 ^ ^ "* 1 + 

^ (n 2 +n 3 )a 2 -l} 

[exp(~nj^a2){ Qjj(c 2 lO)-Q^ 2 (c 3 lO)} - {QjjCc2l 

-* Q(^(c3|nj^c/d)}]+b3 [ L^(c3!n3c/d) j 

-L^( c 2 1 113 c/d )~exp ( n3a 2 ^ ^ *^2 * ® 



'^7=^4 nj^c/d}] « 

By combining these expressions^ we get the required cdf of 
given in equation ( 5*3*4 )- 

For unequal sairple sizes, the ncn~null distribution 
of Tj is very coirplicated* It can be obtained either by proceeding 
as in Theorem 5-3*1 or by the method given by Khatri (1974)* 

For equal sample size case n^ = n^ = ng * n and 0 ^ > 02 ^ 6^ 
the non-null distribution of T^ obtainejd in a similar manner or 
derived from Khatri^ s result is given by 

( 5 . 3 .- 5 ) P [T3 < cie) ]= Q^(c3i0) - 2g3_{ (l+c/d)‘'‘^-(l+2c/d)'''^ 

~ ( 94 + 95 ) Qd^ ^3 1 +2g3Qd( *^3 * 

-(g^+g2)li^(c2lc/d)+2g^ Js^{c^\2c/A) 

+(g2“'9j^^ L^(c3lc/d^- c ^ O, 
where Cj = dTj/c, yj = n(9j^-0j)/o (j = 2,3), 

s= exp(-y 2 ~ 73 )/ 6 / 92 - (-73V2, 93 = exp( 273 ''y 2 )/ 6 / 

= exp(r3)/2, = exp (r^‘-r2)/2, gg = exp 

gj « oxp (2r2-r3)/6, and d, Q^C-l -) and L^C*!*) are given in 
equation ( 3 - 1 * 2 ) • 

One of the objectives of this chapter Is to study the 
performance of the statistic for diffeirent comblnaticxis of 



along with a con^arative study of the statistics 
1 ‘ 3 > U 3 snd for = ^^2 “ ^ 3 * Methods of obtaining the 
recjuirod critical points for all tlicsc tesc procedures are 
given in the next section* 

5 #4 • The critical points of the test statistics « 


From equation (5*2*8)/ the upper lOOa perc^t critical 
point ^ of the test statistic T-. is given by 

X/ Ou J, 


( 5 *4 - 1 ) c 


l,a 


[1 - TSilp ' 


N-n- 1 n 


/ n^^/N < l-tt/ 


whore a is the chosen level of significance. Some critical points 
of are tabulated in Table 5*4»1» 

The critical point of the test statistic Tg is 

given by P [T 3 > = ot* Prom equation (5*2*ll), it is 

clear that, c^ „ is the solution of the Ration 

O / Ou 


(5 • 4 . 2 ) 


J s=0 


_ . - . , n . rstj-ri-ip^r sections/ the sfcat:jLstics a.ii<i 

As we had mentioned in earlier 

Uj are very closely related to each otl,er. Although Khatri (1974) 
has provided a table for the critical points of he has not 

mentioned the procedure for solving eguatJon (5^.2),. Hero »o 

suggest Newton-Raphson method for obtaining by solving the^ 

irr^-t-r:, /->Tntical nolnt as the initiaj. 

equation (5 -l -a) with an approximate critical po 



value. Note that, for a not too small, {l+(j+i)o /d]"^ 

decreases rapidly as j increases. Haice,the appro^te critical 
point may be taken as the solution of the equation Cs .4 . 2 ) 

corresponding to the term j = o. Qjnsequently, 

(5.4.3) ^ 3 ,a = ' 

It is clear from the equation (5*4.2)/ that for K = 2/ ' givc 5 n 

3 /CC 

in equation (5*4.3) is the exact critical point* 


Some exact and approximate critical points are taljulated 
in Table 5 .4.2 for a = 0.05 and for some selected values of K 
and d . Prom the extensive study made in this direction/ we 
conclude that/ except for large K and very small d/ the approxi** 
mation is reasonably good* Note that/ c^^^/d is ^qual to the 
critical point of U 2 tabulated by Khatri (1974) » 


For the statistic Khatri (1974) has tabulated the 
critical points for a = 0.05/ 0-01 and for some selected values 
of nj^/n 2 /n 2 and d» For equal sanple size case/ the critical 
point is the solution of 


j =0 ^ J=° 


= Kot. 


As for T 3 , the solution of this equation for terms corresponding 
to j SB 0 / namely 


( 5 *4 <»5 ) 


(K-.l)(K+2)(l4u)*'^ 2C& 



gives a fairly good 
for largo values of 


approximation to tho exact critical point 
d* This could be used as an initial value 


for solving equation isble 5 , 4 -i gives exact and 

approximate critical points of U 3 for a = 0=05 and selected 
values of K and d. 


As Khatrl (1974) and Singh Ci983) have shown that 

has an ^^-distribution, the critical points of can 

bo obtained easily. 

Some critical points of these test statistics 
U 3 and are tabulated in Table 5-4.1 for K = 3 , n.,=n 2 =^=n^ 
OC a 0.05 and a s= O.IO. These are used for studying the power 
function of these tests. 

5 *5 • Power of the tests . 


The power of the test T^ is P [T^hc,^ c*®]' where 
Cl^a exact critical point given 5 ji equation (5'»4iyl)* 

By making use of the non-null distribution of T^ given by 
equation (5 •3<4), some power values of > ©2 ^ ®3 
evaluated and tabulated in Tables 5 «6 5^5 •‘2 and 5<>5r>3 for 

flc Bs 0 .CS and different combinations of 

Table 5-6.1 shows that^ for fixed values of n 2 /n 2 and 
d, the power of the test increases vexy rapidly as n^ increases. 
But this is not the case for changes in ^ 2 ^ ^3 fixed 

nj^, as is seen in Tables 5.5.2 and 5.5.3. 



The power of the test is obtained by using 
equation (5. 3 *5). Since the derivation of the power functions 
of and Tg involved lengthy calculations/ the siimilated 
power values of these tests are also tabulated along with the 
exact values in Table 5»5*4« This serves as a check for 
theoretical e 3 <p res s ions « 

The power functions expressions of the tests and U. 
provided by Khatri (1974) are extremely conplicated even for 
** case* Hence/ only the simulated power values are 
obtained by using 1000 iterations* Some of these values are 
tabulated in Table 5*5*4* The nature of is entirely different 
from the remaining statistics* It can be used for testing the 
equality of 6 ^'s (i = l/ 2 /*«*,K) against a specified alternative 
I> nax (© 2 / 63 # * (with suitable relabelling if necessary)* 
It can bo seen from Table 5*5*4/ that the power of is 
considerably higher than that of other three statistics* 

Since the power values for all the tests are provided 
only for difficult to compare the 

performance of these statistics* However/ some conclusicais 
can be drawn from Table 5*5 *4* Both tests T^ and perform 
equally well • The tost Uj porfotms slightly better than Tj 
when (e^-e^' while Tj performs better than O 3 if 

(e^-Sj) is small. If ( 63 - 93 ) i^ sraall.then in better than T 3 
while for ( 63 - 62 ) largo, the reverse is the case. Similar 
conclusions are expected for other values of K. The LR tost 
statistic D4 is recommaided for testing against a general 



llti 


alternative hypothesis, since its critical points are easy to 
evaluate from the F~distrlbution even for unei^ual sample sizes « 
Against a specified alternative like > max ie , the 
statistic is recommended* 
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TABLE 5 <»4 *1 * Exact critical 

points of 

the tests 

^3 


and for K = 

3 and n^ = 

: n 2 = = n. 








n = 11, d = 

30 

n = 21, d 

= 60 

X CJSua 

Ot « 0 *05 

a = 0.10 

a = 0.05 

0^ sr O •iQ 


0 .2459 

0 .1780 

0.1261 

f 

0.0913 j 

^3 

3.9036 

3.1125 

3.7878 

3.0402 1 

j 

"3 

0.1236 

0.0977 

0 .0600 

0 .0477 


2*5 2 

2*04 

2 -.45 

1.99 
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table 5*4*2- Exact critical point Cj^^ of in top row 

and its approximated value c* ^ as given in 
ecru at ion (5-4-3^ in bottom row for a = 0*06* 


\ 

d^s^ 

2 

2 

6*9443 

6 *9443 

3 

5 *1433 
5*1433 

4 

4*4590 

4-4590 

5 

4*1028 
4*10 28 

6 

3. *685 3 
3*885 3 

7 

3-7389 

3*7389 

8 

3*6337 

3*6337 

9 

3*5546 

3*5546 

10 

3*49 28 
3*49 28 

15 

3 *3158 
3.3158 

20 

3*2317 

3*2317 

50 

3 *08 7 3 
3*08 7 3 

100 

3*0411 

3*0411 


3 

4 

9 .6651 
10.6491 

11-4509 

13.4919 

6*8972 

7.2599 

8.0316 

8 -7446 

5 -8589 

6 .0595 

6.7533 

7-1326 

5 *3222 

5 .4564 

6*0937 

6*3397 

4 .9959 
5*0959 

5-6931 

5*8716 

4 *7770 

4 *8567 

5*4245 

5*5638 

4*6 201 
4*6867 

5 -2321 

5 .346 2 

4*5022 

4*5597 

5*0876 

5 *1845 

4 .4104 

4 *4613 

4*9751 
5 -0597 

4 *1479 
4 *18 21 

4-6537 
4 *7076 

4*0235 

4.0510 

4*5016 
4 *5436 

3*8106 
3*8 284 

4*2415 

4*2667 

3.7427 
3 .7573 

4 *1535 
4*1793 


5 

6 

12*7893 
15 -8885 

13-8618 
18 .0000 

8.8785 
9-8 266 

9 .5566 
10 .9 248 

7.4181 

7-9628 

7 *9494 
8*6491 

6*6648 

7*0112 

7*1203 

7-5594 

6 *2073 
6*4547 

6*6166 

6 -9 266 

5-9006 

6.0909 

6 *2789 
6*5149 

5*6809 

5*8349 

6 *0369 
6*2262 

5*5159 

5 .6453 

5 *855 2 
6.0129 

5*3875 

5 *499 2 

5-7137 
5 -^489 

5*0205 

5*0893 

5 .3094 
5 .^0'^ 

' 4 *8469 
4-8991 

5 *1181 
5 *1735 

4*5501 

4*5793 

4*7909 
4 ^ 239 

4 *4554 
4 .4795 

4*6867 

4-7129 


8 . 10 , 


15-6267 

21*6643 

16 -8016 
24*8328 

10*6096 

12*5775 

11 -4168 
13-9 386 

8 *7736 

9 *759 2 

9 *4053 
10.6514 

7*8 258 

8 *4337 

8 *366 2 

9 *126 2; 

7*2497 

7 *6722 

7-7340 
8*2 57 

6 -*8631 
7*1804 

7 *3095':> 
7 •6937 

6 .5860 
6-8373 

7-0050 

7*3106 

6 *3777 
6.5348 

6 .7761 
7-0261 

6 -2156 
6*3913 

6*5978 

6.8004 

SrnTSm 

5*85^ 

6 *C^74 
6*2052 

5-6322 
5 *6057 

5*84.55 

5-92951 

5 *1565 
5 *1941 

5 -4314 
5 <47 22 

5 .0368 
5*0653 

5*2993 

5 *3302 




TABLE 5 *4 -Si Exact critical point of in top row and its 

apprbxdjnated value given by equation (5*4*5) 
in bottom row for a a 0»05« 


\ K 

N 

2 

3 

4 

2 

3 *47 21 

3 *47 21 

4.5099 

4.7736 

5.1699 

5 .708 2 

3 

1#7144 

1 #7144 

2.1559 

2*2183 

2.4305 . 
2*5569 

4 

1#1147 

1#1147 

1*3776 

1*4028 

1.5399 

1*5900 

5 

0-B206 
0.8 205 

1*0031 

1 *0164 

1 *115 2 
1*1411 

6 

0*6475 

0 .6475 

0*7858 

0 .7940 

0.8702 

0.8860 


0*5 341 
0*5 341 

0 *6447 

0 *6503 

0.7120 

0.7226 


0 .454 2 
0*454 2 

0.5461 

0*5501 

0.6013 
0.6093 " 

9 

0*3949 

0*3949 

0*4733 

0*4764 

0^207 

0 *5 265 

10 

0*3493 

0 *349 3 

0*4175 

0.4200 

0.4587 

0 .46 33 

15 

0*2211 

0 *2211 

0.2623 

0.2633 

0.2369 

0.2339 

20 

0 *16 16 
0*1616 

0.1910 

0 *19 16 

0*2035 

0*2096 

50 

0 .0617 
0*0617 

0*0725 

0 *07 26 

0*0733 

0.0791 

100 

0 .0304 
0 WD304 

0 .0356 
0*0357 

0.0387 

0*0333 


5 

6 

8 

10 

5 .1317 
6.4833 

5 .26 22 

7 .1650 

6.0097 

8 .3541 

6.2515 

9 .39 23 

2.4610 
2*8 259 

2.5454 

3 .0548 

2.8057 

3.4395 

2*9938 

3*7622 

1.5687 

1.7356 

1*6^0 

1*8574 

1*7717 

2*0584 

1*8933 

2*2237 

1.1390 

1.2368 

1*1812 

1*3162 

1*2803 

1*4457 

1*3661 

1*5508 

0 .8900 
0*9560 

0.9 230 
1*0136 

0.9973 

1*1070 

1*06 25 

1 .18 22 

0.7287 

0*7772 

0.7556 
0*8 220 

0.8147 

0.8942 

0.8668 
0.95 20 

0.6161 

0*6539 

0.6339 

0.6904 

0*6977 

0*7438 

0 . 730 s 

0.7955 

0*5333 

0^6640 

0(65^ 

0.6946 

0*5944 

0*6435 

0.6311 

0.6324 

p*4699 
’ 0 .4956 

0 .4369 
0.5 219 

0.5 231 
0*5639 

0 .5549 

0 .5971 

0.2939 

0.3073 

0 *3044 

0.3231 

0.3260 

0.3473 

0.3450 

0.3663 

0.2136 

0.2229 

0*2211 

0 *2337 

0*2^4 

0*2505 

0 .2499 
0.2638 

0.0307 

0.0133 

0.0335 
0 .0376 

0^90 

0*0935 

0.0938 ! 
O' .093 2 

0*0396 

0.0411 

0*0409 

0.0429 

0*0436 

0*0457 

0.0459 

0..0479 
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TABLE SmSil* Exact power of the test for a = 0.05, 
n^=15, d = 30> K=3 and 0g=O. 



.00 

.00 

.0500 

.0500 

.0500 

.0500 

.0500 

.0500 

.05 

.00 

.0641 

.03 24 

.1059 

.1359 

.1745 

.2240 

.10 

-00 

.08 24 

.1359 

.2241 

.3663 

.5561 

.7207 

.15 

.00 

.1058 

-2241 

.46 24 

.7543 

.9013 

.9497 

.20 

.00 

•1359 

-3674 

.7696 

.9490 

-9827 

.9913 

.25 

.00 

.1745 

.57 S 

.9409 

.9910 

.9970 

.9985 

.30 

.00 

.2241 

.78 23 

.9886 

.9984 

,9995 

.9997 

.05 

.05 

.0584 

.0690 

.CS24 

.099 2 

.1206 

.1477 

-10 

.05 

.0749 

.1137 

.1743 

.2683 

.3949 

.5180 

•15 

.05 

.0963 

.1875 

.3630 

-5989 

.7715 

.8663 

.20 

.05 

*1236 

.3081 

.6474 

.8769 

.9533 

.976 2 

.25 

.05 

.1588 

.4891 

.8763 

.9758 

.9918 

.9959 

.30 

.05 

*2038 

.699 2 

.9712 

.9958 

.9986 

.999 3 

.35 

>.05 

*2616 

.8660 

.9946 

.999 3 

.9998 

.9999 

flO 

•10 

.07 23 

.1074 

.1632 

.2509 

.3694 

.43 20 

.15 

*10 

.09 S 

.1771 

.3395 

.55 26 

.6954 

.7706 

.20 

*10 

*1191 

• 2909 

.6004 

.6035 

.8919 

.9368 

.25 

#10 

.15 30 

.4611 

*8186 

.9417 

.9779 

.9888 

.30 

#10 . 

#1964 

.6596 

.9407 

.9886 

.9961 

.9980 

*35 

.10 

.25 21 

.8 268 

.9864 

.9980 

.9993 

.9997 

.40 

.10 

.3228 

.9301 

.9974 

.9997 

.9999 

.9999 

.20 

.20 

.1169 

.:B45 

.5863 

.7778 

.8503 

•88 36 

.25 

• 20 

.1500 

.4506 

.7912 

.8963 

.9320 

.9488 

.30 

.20 

.19 27 

.6425 

.90 29 

.9561 

.9759 

.9859 

.^5 

.20 

.2473 

.8018 

.9591 

.9870 

.9951 

.9975 

.40 

.20 

.3166 

.90 27 

.9867 

.9975 

.9991 

.9996 

-45 

.20 

.402"» 

.9572 

.9970 

.9996 

.9999 

.9999 

.30 

.30 

*19 21 

.6411 

.8999 

.9504 

•9666 

.9740 

.35 

.30 

.2466 

.7995 

.95 29 

.9770 

.9848 

.9886 

.40 

.30 

.3153 

.8939 

.9783 

.9902 

.9946 

.9969 

.45 

.30 

.4012 

.9517 

.9909 

.9971 

.9989 

.9994 

-40 

.40 

.3157 

.3936 

.9777 

.9889 

.99^ 

-9942 

•50 

*50 

.5006 

.9768 

.99 50 

.9975 

.9983 

-9937 

.60 

.60 

•7112 

.99 48 

.9989 

.9994 

•9996 

•9991 

.70 

.70 

.87 23 

.9988 

.9998 

.9999 

.9999 

•9999 
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TABLE 5 .5 •2* ^act power of the test for a = 0.05^ 
r, . = 1 c; . 3 and Bj = 0 i 


= 15 ^ 



( n^ » 5 ) 


10 


i0500 

.1053 

.2238 

.4427 

#6874 

.846? 

.9 273 
.9656 
-98 38 
.99 23 

.18 27 
..3649 
.588 3 
.7674 
.88 21 

.9435 

.9732 

.9874 

•5732 

•7388 

.8456 
.9143 
.9566 
.979 2 
-8401 

.9038 

.9432 

-9685 

.9412 

.9784 


15 


*0500 

#1058 

.2239 

.4483 

.7134 

.3620 
.9556 
.98 36 
.9940 
.9978 

.1973 
•3978 
.6481 
.8 319 
.9300 

#9732 
.9901 
*9963 
.6422 
*8 206 

.9155 
.96 21 
.9844 
.9940 
.9141 

.9595 
.9811 
.9915 
.98 CB 
.9957 


20 


.0500 
.1058 
• 2239 
.45 20 
*7283 

.9021 
.9701 
#9913 
.9975 
.999 3 

•205 3 
.4161 
.68 28 
.8638 
.9547 

#9860 

.9959 

.9988 

.6803 

.8639 

.9485 
.9815 
.99 38 
.9981 
.948 2 

.9808 
.99 30 
.9975 
.99 29 
.9990 


25 


.0500 

.1058 

#2239 

.4538 

.7375 

.9142 

.9780 

.9949 

.9989 

.9997 

.2102 
.4 273 
.7043 
.8911 
#968 2 

.9920 
.9981 
.9996 
.703 2 
.8890 

.9655 
.9901 
.9973 
.999 3 
*9654 

.9899 
.9971 
.999 2 
.9971 
.9998 


30 


.0500 

.1058 

.2240 

•4550 

.7434 

.9 219 
.98 26 
.9967 
.9994 
.9999 

•2134 

-4346 

-7184 

.9053 

.9761 

.9950 

.9991 

#9998 

.7179 

.9043 

.9749 

*9942 

.9987 

•9997 

.9749 

.9941 

.9987 

.9997 

.9987 

.9999 
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TABLE 5*5-2 Contd; 


( n 2 = 30) 


®1 

®2\ 

5 

10 

15 

20 

25 

30 

.00 

.00 

.0500 

.0500 

*0500 

.0500 

*0500 

.0500 

.05 

.00 

.1058 

*1058 

*1058 

.1059 

.1059 

.1059 

*10 

.00 

.2240 

^2240 

*2240 

.2240 

.2240 

.2240 

*15 

.00 

.4550 

.4558 

.456 3 

.4567 

.4570 

.457 2 

-20 

.00 

*7434 

.7474 

.750 2 

.75 22 

.7536 

.7547 

*25 

.00 

*9 219 

•9 269 

.9303 

.93 26 

.9343 

.9354 

• 30 

.00 

•98 26 

.9854 

.9871 

.9883 

.9890 

-.9896 

*35 

.00 

*9967 

.9977 

*998 2 

.9985 

.9987 

.9989 

*40 

.00 

.9994 

.999*7 

.9998 

*9999 

.9999 

.9999 

*10 

*10 

.0995 

*1315 

*1531 

.168 2 

•1790 

.1371 

*15 

*10 

.2070 

.2722 

*316 2 

.3469 

.3690 

.3853 

*20 

*10 

.387 3 

*4914 

.5603 

.6074 

♦6406 

.6645 

*25 

•10 

.6190 

.7209 

.7858 

.3 283 

.8570 

.8768 

.30 

*io 

.834 2 

.8903 

♦9 250 

.9463 

*9593 

.9637 

.35 

•10 

.9514 

.9711 

.98 21 

.938 3 

.99 21 

.9944 

.40 

*10 

*9894 

.9945 

.9969 

.993 2 

.9989 

.999 3 

*45 

*10 

.9930 

.9991 

.9996 

*9993 

.9999 

.9999 

.20 

.20 

*3118 

*4573 

.5444 

.5998 

.6369 

♦6627 

.25 

.20 

.4685 

.65 34 

*7545 

.8135 

.8498 

.8732 

.30 

*20 

.613 2 

*7966 

.88 26 

.9 267 

.9505 

.964 2 

*35 

.20 

.7677 

.8954 

.9498 

.974 2 

.9858 

.9915 

.40 

*20 

.8994 

.9597 

.9831 

.99 25 

.9965 

-998 2 

.45 

*20 

.9705 

.9894 

.9960 

.9984 

.9993 

.9997 

.30 

.30 

.57 24 

.7841 

.8791 

.9 257 

.9502 

.9641 

.36 

.30 

.6765 

.8705 

.94 28 

.97 22 

.985 2 

.9913 

.40 

•30 

.768 3 

.9 250 

.9736 

.9899 

.9957 

.9980 

.45 

.30 

.8591 

.9615 

.9888 

.9965 

.9988 

.9996 

.40 

.40 

-7406 

.9 204 

.97 28 

.9897 

.9957 

.9980 
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TABLE 

5 •S ♦S • 

Exact 

power of 

the test 

T^ for a 

= 0.05, 

0, 

Hi = 

30^ Up = 

20, Hg = 

15, K = 3 

and 8,^ = 








V 

5 

10 

15 

20 

25 

30 

CD 

CD 

Z 







-00 

-00 

.0500 

.0500 

.0500 

.0500 

-0500 

-0500 

-05 

•00 

.2129 

.2213 

-2231 

.2237 

-2239 

.2240 

•10 

-00 

.546 3 

-6347 

.6736 

-6959 

.7104 

-7 207 

•15 

-00 

-818 2 

-9074 

-9317 

•9416 

.9467 

-9497 

• 20 

•00 

•9438 

-98 23 

-9880 

-9898 

-9907 

-9913 

-25 

*00 

.9854 

.9969 

.9979 

-998 2 

.9984 

-9985 

•30 

.00 

•9966 

-9995 

-9996 

-9997 

-9997 

.9997 

•10 

•10 

•3773 

-4309 

.4541 

•4673 

-4759 

•48 20 

*15 

•10 

.6486 

-7 250 

-749 2 

.7603 

.7666 

-7706 

• 20, 

•10 

•8465 

.9064 

-9 224 

.9 297 

.9340 

-9 368 

• 25 

•10 

-9475 

•9787 

.9847 

-9870 

-9881 

.9888 

-30 

•10 

•985 2 

.9960 

.9973 

.9977 

.9979 

.9980 

• 35 

•10 

•9964 

•999 3 

*9995 

-9996 

.9996 

.9997 

•40 

•3,0 

•999 2 

#9999 

-9999 

.9999 

.9999 

.9999 

• 20 

• 20 

•3088 

-8609 

-8734 

-8787 

.8817 

•8836 

• 25 

-20 

•9097 

♦9380 

.9440 

•9465 

-9479 

.9483 

• 30 

• 20 

.96 35 

.9791 

-98 27 

•9843 

.9853 

•9859 1 

• 35 

• 20 

-9879 

-995 2 

-9966 

.9971 

.9973 

•9975 j 

•40 

• 20 

-9966 

•9991 

.9994 

•9995 

.9995 

-9996 

•45 

-20 

.999 2 

-9998 

-9999 

-9999 

.9999 

.9999 

• 30 

-30 

.9551 

.9689 

.9717 

.97 29 

-9736 

• -9740 

^35 

•30 

-9795 

.986 2 

.9875 

.9381 

.9884 

♦98S6 

#4 0 

»30 

#9918 

-995^ 

-9961 

.9965 

.9967 


#45 

•30 

•997^ 

-9989 

-999 2 

-9994 

.9994 

♦9994 

•40 

-40 

.9899 

.99 31 

.99 ^7 

-9940 

.9941 

♦9942 

• 50 

•50 

-9977 

•9985 

•9986 

-9987 

.9937 

♦9987 

•60 

•60 

•9995 

.9997 

.9997 

.9997 

.9997 

♦9997 

•70 

•70 

-9999 

-9999 

•9999 

.9999 

.9999 

♦9999 
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table 

5.5 .4 

• Exact 

_and simulate powers of the 

tests 

and T^, 

and simulated 

powers of 

the tests 

Ug and 

for 


a = 0 

.05, 83 = 

(n = 

0 and n^^ 

= 11 , d = 

="2=^3 =n* 

30 ) 



CD 

®2 

h 


^3 

3 



.00 

.00 

•050 

.066 

.050 

•053 

.049 

•053 

.05 

.00 

.087 

•099 

•056 

.054 

.058 

.064 

.10 

•00 

•150 

•147 

•076 

.078 

.075 

•032 , 

•15 

.00 

• 260 

• 271 

•HI 

•1^ 

.129 

.156 

.20 

•00 

•441 

•469 

•173 

.200 

.211 

• 226 

«25 

.00 

•669 

•674 

.279 

6 

.319 

• 266 

.?0 

.00 

.850 

.846 

.439 

.439 

.49 2 

.392 

.^5 

.00 

.944 

•940 

.6 29 

.6 29 

.65 3 

<•601 1 

^4 0 

.00 

•981 

.980 

.795 

.733 

•629 

.654 

•50 

.00 

-997 

•998 

.961 

.960 

.97 2 

.890 

.10 

•10 

•117 

.112 

.08 2 

•07 2 

•063 

•115 

•15 

•10 

.202 

. 206 

.107 

•109 

•100 

•165 

• 20 

•10 

•345 

•342 

•153 

•149 

*141 

.236 

• 25 

•10 

•541 

•518 

• 2^2 

• 247 

.254 

•319 

• 30 

-10 

•731 

•734 

.358 

•371 

•411 

-1461 1 

•35 

•10 

•870 

.869 

•5 22 

•5 28 

•5 29 

•6 30 

.40 

•10 

•947 

.947 

•691 

.667 

•731 

.660 

•50 

•10 

•99 3 

•99 3 

Sll 

•917 

.940 

..858 

.20 

•20 

.334 

.338 

.207 

• 225 

•169 

•4 29 

.25 

• 20 

•5 21 

•516 

• 271 

• 277 

.236 

•5 26 

i 

• 30 

• 20 

•699 

.688 

•376 

.376 

•376 

•6 29 

.35 

• 20 

.8 27 

•8 25 

•519 

•617 

•5 20 

♦760 

.40 

•20 

•908 

•906 

.668 

.674 

•719 

•769 

•50 

• 20 

.98 2 

.934 

.88 2 

.889 

.915 

.884 

.30 

.30 

.696 

.678 

.469 

.46 3 

.386 

*76 2 


[ 


^Simulated powers baseid on 1000 samples for each sartple size 
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TABLE 5*5.4. Contd. 


(n = 21, d = 60) 


®1 


^1 

rn'iJ- 

1 

^3 

T* 

3 

“3 

< 

.00 

.00 

.050 

.042 

•050 

.048 

.051 

•043 

.05 

.00 

•143 

•135 

.073 

*064 

.066 

.072 

•10 

.00 

.407 

•416 

.160 

•174 

.175 

•183 

.15 

.00 

.849 

.853 

•403 

.39 2 

.433 

.360 

• 20 

.00 

.981 

•980 

•804 

•810 

.866 

.660 

» 25 

•00 

.998 

.999 

.969 

.976 

.989 

.906 

.05 

.05 

•112 

-111 

.079 

•073 

.075 

.090 

•10 

#05 

.318 

.319 

•143 

.154 

.157 

*180 

•15 

.05 

.717 

.716 

.3 28 

.329 

.328 

.335 

• 20 

.05 

.946 

.947 

.68 2 

•675 

.679 

•670 

• 25 

.05 

*99 3 

.993 

.9 24 

.9 27 

.932 

■#9 23 

• 10 

•10 

.307 

.313 

.19 3 

.199 

.146 

.400 

*15 

♦10 

•686 

•686 

.351 

.349 

.347 

.6 28 

*20 

•10 

.900 

.913 

•661 

.661 

.699 

.783 j 

• 25 

*10 

.981 

.986 

.88 2 

.890 

-908 

.898 1 

.30 

•10 

•998 

.997 

.974 

.98 2 

.987 

.958 

*20 

«20 

*888 

#864 

•756 

.749 

.701 

.898 

.25 

.20 

*961 

.960 

.68^ 

.882 

•839 

*944 

•30 

*20 

-988 

•935 

.957 

.950 

.956 

.966 

#30 

.30 

.936 

.939 

.969 

.975 

.974 

.993 


*Siinu.lated powers iDSised on 1000 sainples for each, ssn^le size* 




CHAPTER VI 


GENERALIZED STATISTICS FOR THE EQUAL SAMPLE 
CASE v\HEN ONE OBSERVATION IS MISSING ON THE LEFT 


6*1* Introducbion and test sta-tisulcs * 

In Chapter V, tests for the equality of location 
parameters of K (> 3) populations are discussed, when the 
smallest observation is available in each sample- Here the 
problem is extended to the case v/hen the smallest observation 
is missing but the second smallest observation is available 
in each sample of. size n, that is. 


X 


(i) v-Ci) v(i) 


(i = 1,2, 


are the available observations from the ith population with 

(n~Sj_) > 2* 

Similar to T^ of Chapter V, the test V^ defined by 

(6*1.1) V = ~ min 

^ 2<i<K 

is proposed for testing H^ : 0^^ = © 2 = ••• = ©k = © against 


H.. : 0. > max (8 j), where 
^ ^ 2<i<K 

K ^“^i 


(6.1.2) do* = 2 { 2 

i=l j=2 i ^=1 


(n-'S,-2)» 


The generalization of statistics T 2 and U^ of 
Chapter V are V 2 


V^ and V^ respectively. These for equal sample 
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size case arc given by 


(6.1 .3) 


and 


V. 


= { max 
l<iS^ 



min 

l<i« ^ 


(6.1.4) ^3 = [ ^), (X^^’ ^ )} ]/0*. 

2<j<IC ^ ^ 

Both of these are proposed for testing against : atleast 
one is different from 6. Corrpared to two-sample case (K = 2)y 
the LR test is much more complicated oven for K s= 3. Even the 
derivation of ML estimates under is far from simple and 
requires a full study in itself. Consequently^ we have not 
studied the LR test in this case and have left it as an open 
problem. 


The test procedure is to reject if where 

the constants Cj^ are obtained from solving the equations 


(6.1.5) P [Vi > = a, 

whore a is the chosen level of significance. The required null 
distribution of these statistics are derived in the next section. 

6.2. Distribution theory • 

Lemma 6.2*1. Lot ( i = 1^2, .-»/K) bo K i.i.d. random variates 
with pdf 

(6.2.1) f(y) = n(n-l) [exp{-(n~l)y}-oxp(~ny)]^ Y > 0/ 
then Z = {Y^-min(Y 2 ,Y 2 / ...^Yj^)} has the pdf given by 
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( 6 . 2 .4 ) 


.exp [-{(n-l)(K-i)+j}z 2 ] 


• [e555{«(n-l)z^}-e3<p(-n2^)j , z^^Z2 > 0. 


Make a transformation z = z^-z^, Z 2 = z^r then the range of the 
transformed variables are 22 > max ( 0 ^~z), -00 < 2 < «,, 

Thus, the marginal pdf of Z is 


(6*2. 5 ) fiz) 


S fCz^, z+z^) dz 

~z 


2 # io'r'^ 2 ^ ^^2 ^ z < 0 


; f(z.,^ z+z«) dz-, , 2 > 

L. o i Z - 

where f ( ., .) is given as in equation (6.2.4). On siirplification, 
the equation (6*2.5) gives the required density function of Z. 


Lemma 6.2*2. Let (i = 1,2, .«.,K) be K i.i.d. random variates 
with pdf given as in equation (6-2-1), then the pdf of 


Z = max (Y, ) - min (Y.) is 

(6.2.6) f(z) = K S^(-l)^''J‘"^(^T^)n^’*‘^(n*-l)^‘’^{ 

j=0 J 

K"i-2 r • ..n~3iz./ -,N.—(n~l)z 

• { 1-e } ^ {n(K~j-l)e +(n--l)je' 

- ( nK-^n- j ) e^ ^ , z > 0 . 

Proof > The cdf of Z is (David, 1981, p* 26) 

00 ‘K’ 

( 6 . 2 . 7 ) G(z) =k; f(y) {P(y+z) - P(y)} dy. 
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where f( *) and F( ») are the pdf and cdf of given in eguations 
(6.2-1) and (6.2*3) respectively. Substituting for f( .) and 
F( .) in equation (6.2*7)/ expanding the term £ F(y+ 2 )''F(y)} ^ 

as a binomial sum and integrating w-r^ to y we get 


K-1 

(6*2*8) G(z)= S 0 .( 2 ) 

J=0 ^ 

where 


(6.2.9) Gj( 2) = K(~l)^~^'"'^(^T^)n^'‘’^(n-l)^“^’ 

{ l-e^ n-1 ) Zj j j 

• (nk-i- jTr^-jT- • 

Taking logarithm on both sides of equation (6.2*9) and different- 
iating with respect to z, we get 


(6 . 2 * 10 ) 


3Gj(2) 

~dz 




G j ( z) . 


Differentiating both sides of equation (6'«2.8) v/ith respect to 
Z/ and substituting for 3Gj(z)/3z from ec[uation (6 » 2*10), we 
get the required pdf of Z given in equation (6-2*6). 


Theorem 6*2*1. The null distribution of the statistic 
defined in equation (6.1*l) is given 


( 6 .2*11) 


f(v^) 


fl(Vi) / < 0 

fjCVi) , > 0, 


where 
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K-2 


J —0 


1^2^ 


- fl-J5V'*l'''^'^^V(J2j3)L 

2 

f2(v2) = (K-1) _I1^ (-l)-^'(^T2)n^''J(n-l)J+2^ 

- f 1+nv^/d}’”^'^'^^ V( ^2^3^^ ' 

and (1 = 1, are given in equation (6-2»2)» 

Proof * Note that, ~ (i = Iy2/»«*,K) has the 

distribution given as in equation (6»2»l)» Hence, under 
= dVW, where 2W = 2d0*/d has a distribution and 

Z = - min (Y 2 / Y^, • »*,Yj^) has the pdf given in equation (6.2*2)- 

By writing the jpdf of Z and W, and making a transfortnation 
Vj_ = dz/w, w = w, and integrating w.r. to w, we get the required 
pdf of Vji^ given in equation (6.2.11) * 


For K = 2, equation (6.2*11) simplifies to 


fCvj,) = 


o 




which agrees with the null distribution of T, given in equation 
( 3 . 2 . 7 ) for r^ = r 2 = 1 and n^ = n^ = n case. 

For K = 3, the null distribution of is given by 



134 


(6 *2-12) f(v^)a 


® 2^ 2^ n-l )vj/d} l-C 2n~l )v^/d} 

+ B^£i-2nvj/d}“‘^^j ^ < 0 


L®1 [%^ l+nvj^/d}''^^^] ,Vj^ > O, 

where = 2n(n-l)/{ 3( 3n-l ) ( 3n~2)} , = n^(3n-l), 

B 3 = 3n(n-l)( 2 n~l), B^ = (n-1 ) ^( 3n~2), 

Bg s= n(5i>'3) and Bg = (n-l) (5n—2) • 

It is easy to establish the following relations smong B^*s : 

B^ *■ B^ + B^ * Bg Bg * 4n*" 2 f 


B, 


Si { 


®3 ®4 1 

+ = T , 


and 


1 ^ 2 n -2 2 n~l ^ 2 n 


- # 


2 

= 


These relations can be used for showing that f(v^) given in 
equation ( 6 * 2 * 12 ) is continuous at = 0 and it is indeed a pdf 

Since the null distribution of statistics V 2 and V 3 are 
complicated/ only the case K = 3 is considered and discussed in 
the following theorems : 


Th eorem 6*2*2* Under H / the distribution of V 2 for K s= 3 is 


— d 1 \ '-dr-l 

(6*2*13) f(v 2 ) = [-B 2 £ l+2(n-l)vy'd} ''■^+B 3 £ l+( 2n~l)v^d} 

l+2nv/d}~'^~^+Bg£ l+(n~l)v/d}"'‘^"^ 

*" Bg£l+nv^d} ^ / ^2 “ 



135 


where (i e are given In equation ( 6 » 2 * 12 ). 

Iroof- Note that, (i = 1 , 2 , 3 ) are i.i^, 

random variates with coirmon pdf given in equation (6*2*l)» Now, 
from Lemma 6 * 2-2, for K = 3 , the pdf of 

Z = max (y^,y 2 ,Y 3 ) - min (Y^,Y 2 ,y 3 ) 

can be written as 

f(z)=3B^ [-B2e2'n-l>^+B3S< 211-1) z_B^ 52 nz^B^g(n-l)z_ 3 ^ 5 nzj^ ^ ^ 

Under ^qi^2 ~ where W = d 0 */cr. Nov;, making the transformati 

V 2 = dz/w, w = w, and integrating w»r«to w, we obtain the 
required pdf of V 2 given in equation (6 #2 #13)* 

Theory 6*2»'^« The null cdf of the statistic V 3 for K = ^ is 

(6.2-14) P [V 3 < ciHq] = HA^{H.nc/d}”'^-A 2 {l+(n-l)c/d}”'^ 

t A^{ l-(- 2 nc/d}~^~A 4 { l+( 20 - 1 ) c/d} '"^+A^{l+ 2 (n'*l) c/d} 

+ Ag{ l+3nc/d}‘'‘^-A^{ i+( 3n-l )c/d}‘'‘^-AQ£ l+( 3n-2)c/d}“'^ 

+ AgC l+3(n~l)c/d}" , c > 0, 

vrhere A^=:(n-l) [l+(n-l)/( 2 n-l)~(n-l) V^“’^V(^’^2)+2n^(n-l)/(3n-l)] , 

A 2 =n|]l-f-n/( 2n~l)"*(n-l ) ^/( 3n~l )-n V 3+2n(n~l )^/(3n— 2 ) ] , 
A^zzCn-l) Vc 3(3n-'l)} , A^=2n^(n~l) V? (3n-l) (3n~2)}, 

Ag =:n V{ 3 ( 3 n~ 2 ) } , Ag =( n- 1 ) ^ [l+ 2 ( n- 1 )/ 3- 2 n V C Sn-l ) ] , 
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A7=2n(n~l)2 [l/(2n-.l)~iv'(3n-2)+(n-l)/On-.l)J , 
AQ« 2 n^(n~l) [l/(2n-l)-rv^(3n~2)+(n-.l)/(3n-l)] , 
and Ag = [1-2^3 + 2(n~l)^/(3n-2)] . 

S£S£^* (i = 1/2,3)* Note that, are 

3, •i^d* random variates with pdf f(y) and cdf F(y) as given in 
equations (6*2*1) and (6*2*3) respectively* 


Let z = max { (Y^-Y^) , (Y^-Y^), (Y^-Y^), (Y 3 -Y^)} • Thus, 

P [Z < zlH^] = P [Yj~Y 2 < 2 ,Y^-Y 3 < 2 / Y^-Y^ < z, Y 3 -Y^ < z] 

= P [Y^ > Yj-z, Y 3 > Y^-z, Y^ < Y^+z, Y 3 < Yj+z] 

00 

= / P [y^-z< Y 2 <Yj_+Z/Y^~z < Y 3 < Y^+zlYj_=y] f(y) dy 

•—00 

00 

(6*2*15) = / Q(z,y) dy, 

•^CO 

where Q( 2 ,y) = f(y) p .(y+z)-p.(y-.z) ]• Since f(y) = 0 for 

3=2 

y < 0, hence, 

z 

(6.2*16) P [z < zlH^I = I Q(z/y)dy + / Q(z,y) dy. 

O 2 

Note that, the first term of the equation (6*2*16) reduces to 

; f(y){P(y+z)} ^ dy. Now substituting for. f(y) and F(y), we 

o 

get the simplified form of the equation (6*2^6) as 

/ . - -A (n-l)z^ -2nz ^ ~(2n-l)z 

(6*2*17) P [Z < zlH^] = 1 +A^e -^ 2 ^ +^ 3 ® 

a > o. 
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Under _ dZ/W, where 2W = 2da*/a has a distrilxition » 

Thus, 

00 

(6*2*18) P [v^ £^18^]=/ P[Z< cW/dlW=:w] e"*^w*^”’^dw/(d--l) I • 

o 

Substituting from equation (6*2*17) and simplifying the resulting 
ejspression, we get the required null cdf of V^* 

6 *3* Moments of the statistics under for K * 3* 

^--d~l 

Note that/ the pdf of Vj_ contains factors like (1-av; 
and ( l+av)""*^"^# This allows the evaluation of momoits of by 
using the following lerrma : 

Lemma 6 *^ *1 * For a > 0 and 0 < h < d 


(i) 

J v^( l-av)”'^”^ dv 

= (-1)^ B(h+l,d-h)/a^‘^^ 


^OO 


and 



( ii) 

; v^(l+av)*''^~^ dv 

n 

= B(h+l,d--h)/a^'*'^# 

Proof * 

Making a substitution t = (l-av) we have 


o 1 . (d-h~l) 

/ dv = / a 




= (-1)*' BCh+l^d-hVa^+l. 

Part (ii) of the leiana can be proved by making the aubrtitutlon 
V s= -u in part ( i) • 
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Now^ from Leirtma ei"**! and equation (6«2*12)^ the 
hjU} moment of about zero for K = 3 and h < d is given by 

E{v^) = B^;B(h+i,d-h)d*'« [(-!)»' B^(2n-2)*’+U{.i)h B3/(2n-l)*'+l 

+ (-1)^ B^(jn)''+l+E^/(n-l)^+UBg/n''+l] . 

Expressions for the moments of and can be written down 
in a similar manner# 

R emark 6 « 3 *1 • It is easy to see from the pdf of and V 
g iv en in equation s (6»2*12) and ( 6 • 2 #1 3 ) respectively, that 

E(V2) = 3 E(V^). 

This can also be justified from the fact that 
V 2 = Wj + VJ2 + 

where W. = min . (i = 1^2,3) and = V^# 

6 r 4 • Critical points of the tests for K » 3 « 

The upper lOCa percent critical point c^ of the test 
statistic is the solution of equation (6 #1*5)* From equation 
(6.2#12), we have F,. (o) = 4 and h®ice c^ is either the solution 

o£ 

(6.4.1) B^ tl+(il-l)Vdr'’ - -I fl+ncj^/d!"^]. a 

or of 

(6.4.2) 

* £ l-2ncj./d}"^ ] = a 
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according as a < yj, or a > y, respectively. Pro„ equation 
(6.2.13), Cj is the solution of equation 

(6.4.3) 3B^ [- ^ ^ £l+(2n-l)cy'di~* 

B 

~ to Cl+2n cyd} ^ { l+(n~l)cyd}~^ 

Bg 

rT ^ 1+ncyd} ] = a* 

Similarly, the critical point c^ is obtained from the 
equation (6 -2 •14)* 

Some critical points of V^,V 2 and are tabulated in 
Tables 6*4«1/ 6*4*2 and 6*4*3 respectively for a = 0*05 and K = 3* 

These are calculated by using Newton Raphson method with 
approximate critical points as the initial values* Approximate 
values are obtained by using the fact that the null distribution 
of { Vj^-E(V^)}/{ Var(Vj^)} is approximately normal. This gives 
the approximate critical points as 

^app* = 4£Var(Vi)}^/^ + E(V^) (i = 1,2,3), 

where c^ is the upper lOOapeitent of N(0, l). Some numerical 
calculations show that upper lCX)a/2 percent point of N(0, l) 
distribution gives closer approximation for V 2 and than the 

100 cx percent point* This may be due to the fact that the 
distribution of is in the interval (-<»,~), whereas the 
distribution of '^2 ^3 confined to the interval CO,^). 
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Some exact and approximate critical points for all the 
three tests are tabulated in Tables 6 . 4%4 for a = 0.05 and K = 3 - 
It is clear from the Table 6i4.4, that the approximate critical 
points are reasonably good for large values of d. 

6-5. Performance of the tests * 

It does not appear simple to evaluate the non-null 
distributions of these statistics. Consequently, we use Monte- 
Carlo techniques for the calculation of power of these tests. 

Some of these values based on 1000 iterations are tabulated in 
Table 6*5*1 for a = 0*05 and K = 3 * 

On the basis of these calculations, the statistic is 
recorrmended for testing against a specified alternative like 
Hi : Gi > max ( 62 / 63 ) • conclusion is similar to the use of 

of Chapter V against testing against the alternative 

H 2 / it is observed that V 2 performs better if 6^-^2 vejry small, 
otherwise the performance of is better* 

In Table 6*5.2, the power values of the test T^ (of 
Chapter V) and are tabulated for a s= 0 * 05 , K = 3, nj^=n2=n2=n=16 
and d = 6(12)42* These calculations show that there is a 
considerable loss of power due to censoring of the smallest 
observation* This highlights the Importance of smallest 
observations for testing the equality of location parameters. 
Similar conclusions were drawn by Greenberg and Sarhan (1962) 
for the estimation of location parameters* 
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TABLE 6-4:*l* Exact cx'itical points of the test 
for a = b »05 and K = 3.' 


3n-6 


1.39 38 
1 *07 21 
0.37 22 
0*7357 
I 0 *6 36 3 
0*5607 
0*5012 
0*4532 
0*4135 
0*3803 
0 *35 20 
0 .3 277 
0 *3065 
O *2878 
0*2714 
0*2567 
0 *2435 


0*9366 
0*80 26 


0.6770 

0*6340 




0 *5855 

0*5483 

0*5394 



0.5159 

0*4831 

0*4753 



0*4612 

0 *4319 

0*4249 



0 *4170 

0*3905 

0 *3841 

0 *3739 


0 *3805 

0 * 356 3 

0.3506 

0*^12 

0*3394 

0 *3499 

0*3277 

0*3224 

0 *3138 

0*3107 

0 *3239 

0 *3033 

0*2984 

0*2904 

0 .:^66 

0*3015 

0*26 23 

0 *2778 

0*2703 

0*^59 

0*23 20 

0 *2640 

0*2598 

0 * 25 23 

0*2480 

0*2648 

0*2480 

0-2440 

0 * 2375 

0*2324 

0*2497 

0*2338 

0 *2300 

0*2239 

0#2187 

0 *2361 

0 *2211 

0*2176 

0*2118 

0*2065 

{ 

0 *2240 

0 * 2098 

0*2064 

0*2009 

0*1955 
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TABLE 6<'4»2» Exact critical points of the test V 2 



for a = 0*05 

an d K = 

3. 



d 

6 

9 

15 

18 

27 

3n-6 

n 







4 

2.0681 






5 

1 *5888 

1.419 3 





6 

1.2919 

1 .1540 





7 

1*089 3 

0.97 29 

0 .889 2 




8 

0*9419 

0.8413 

0 .7689 

0.7518 



9 

0 *8 299 

0 .7412 

0.6774 

0*66 23 



10 

0 *7417 

0 .66 24 

0 .6054 

0.59 20 


i 

1 

11 

0.6706 

0 *5989 

0.5473 

0.535 2 

0.5155 

i 

1 

1 2 

0*6119 

0 *5465 

0*4994 

0*488 3 

0.4704 

0 .4669 ; 

I 

13 

0.56 27 

0 *5025 

0.4593 

0*4491 

0*4326 

0*4268 

14 

0*5 209 

0.4651 

0-4 251 

0 *4157 

0*4004 

0.3^30 

15 

0 .4848 

0.4^29 

0.3957 

0 .3869 

0*37 27 

0.3642 

16 

0 .45 34 

0 .4049 

0*3701 

0 .3618 

0.3486 

0.3394 

17 

0.4259 

0.3803 

0 .3476 

0.3398 

0.3274 

0 .3177 

18 

0 .4015 

0.3585 

0.^276 

0.3204 

0.3086 

0.3987 

19 

0 .3797 

0 .3391 

0.3099 

0.3030 

0.2919 

0.2818 

20 

0.360 2 

0.3217 

0.2940 

0.2874 

0.2769 

0 *2667 




— -- 

^ . - 
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table 6 4i4«3» Exact critical points of the test 
for "q = 0»05 ajid K = 3 » 


d 

6 

9 

15 

18 

27 

3i>-6 

n 








4 

1.9390 






5 

1 .4895 

1.3349 





6 

1.2111 

1.085 2 





7 

1.0 211 

0.9149 

0 .8 387 




8 

0.88 30 

0.7911 

0.7 25 2 

0.7097 



9 

0.7779 

0.6970 

0.6 389 

0.6 252 



10 

0.695 3 

0.6 229 

0.5710 

0 .5588 


t 

1 

11 

0.6 236 

0.563 2 

0.516 2 

0.505 2 

0*4873 

! 

i 

12 

0.5736 

0 .5139 

0.4710 

0.4610 

0 .4447 

0.4415 

i 

13 

0 .5 275 

0.47 26 

0.433 2 

0.4239 

0.4089 

0..4036 

14 

0 .488 2 

0.4374 

0 .4009 

0 *39 23 

0 .3735 

0.3718 

15 

0.4544 

0.4071 

0.373 2 

0.365 2 

0 .35 23 

0.3446 

16 

0.4250 

0 .3808 

0 .3490 

0 .3415 

« 

0.3^5 

0.3 211 

17 

0 .399 2 

0 .3576 

0.3278 

0.3208 

0 .3094 

0.3007 

18 

0.3763 

0 .3371 

0 .3090 

0.3024 

0.2917 

0.23 27 1 

19 

0.3559 

0 .3139 

0 .29 23 

0.^60 

0.2759 

0.2663 

20 

0 .3376 

0.3025 

0.277 2 

0*2713 

0.2617 

0 * 25 25 
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TABLE 6 *4 *4 • Comparison of exact and approximate critical 
points of and for a=0^5 and K = 3# 


n 

A 



'^2 


^3 


CZ 

Exact 

Approx* 

Exact 

Approx. 

Exact 

Approx- 

10 

6 

•5012 

•4853 

•7417 

•7740 

•6953 

•12SQ 

10 

12 

•44 26 

•4156 

.6.26 2 

.6 371 

•5399 

*5991 

10 


*4 249 

•39 68 

•59 20 

.6004 

.5588 

•565 2 



*416 3 

-3881 

*5756 

•5833 

.5439 

•5494 



*2435 

• 2358 

.3602 

.3758 

*3376 

*35 24 ^ 

20 


•2150 

• 2019 

• 3041 

*3094 

-.2864 

.2909 

20 


• 2064 

*19 

•2874 

• 2915 

1 .2713 

*2744 

20 

24 

-20 22 

•1885 

• 2795 

-28 3 2 

• 2641 

• 2667 

20 

30 

*1998 

•1861 

• 2748 

• 2784 

-2599 

• 26 23 

20 

36 

•198 2 

•1845 

• 2718 

• 275 3 

.2571 

.2594 

20 

42 

•1970 

*18 34 

.2696 

• 2731 

• 2551 


20 

48 

•196 2 

*18 25 

.268 0 

• 2715 

• 25 36 

•2559 

20 


•1955 

•1819 

-2667 

• 2702 

-25 25 

>254:7 
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TABLE 6»5*1* Simulated power values of the tests V. 

- - ri JL <ftt 




and 

Vj for a * 

0 .05^ K « 3 and 

%■ 

0 

It 

ro 

CD 





(n 

= ll) 





®2 

d = 15 

d = 27 



^2 

^3 

^1 

^2 

^3 .. 

.0 

.0 

•061 

#05 2 

.054 

.061 

.045 

.045 

•1 

•0 

•114 

.058 

.056 

.113 

-066 

.065 

-2 

.0 

.245 

*101 

.105 

.248 

.106 

.101 


#0 

.468 

.233 

.249 

.504 

.235 

.257 

.4 

.0 

.674 

.361 

.413 

.730 

.407 

#450 

#5 

.0 

.844 

.558 

•611 

.897 

.638 

.691 

-6 

.0 

.951 

.739 

•lie 

♦966 

.818 

-864 

•7 

•0 

.98 2 

.864 

.89 2 

#090 

#0 28 

.941 

.8 

•0 . 

.999 

.949 

.970 

.998 

.978 

.986 

•9 

.0 

•999 

.978 

.989 

.999 

.995 

.996 

•2 

*2 

•147 

•108 

♦076 

.149 

.119 

.083 

•4 

• 2 

*5 37 

.^5 

.309 

.59 2 

.326 

.34 2 

•6 

.2 

.380 

.608 

.65 2 

.913 

.679 

.734 

*8 

1 

1 

•2 ! 

.99 2 

.883 

*9 20 

.99 3 

.951 

.964 

•4 

*4: 

.5 29 

#369 

.304 

•SSI 

*409 

.345 

•6 

■#'4: 

; *844 

.611 

.656 

#868 

.677 

#711 

•8 

#4 

.97 2 

.871 

.893 

.980 

.911 

^.0 35 

*6 

#6 

.8'50 

.710 

•<1655 

.871 

.770 

.716 

•8 


-960 

^81 

.897 

.968 

.9 26 

.934 

^8 

»8 

.965 

.913 

.904 

.078 

.940 

.9 25 



— - - 


■ ““ 
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TABLE 6.5 Contd. 


(n s! 16) 


01 ®2 


• 2 *0 
• #0 
•4 -0 
•5 »0 

•6 *0 
* 2 *2 
•4 *2 

*S *2 
•4 *4 

•6 *4 

•6 *6 
•8 *6 
.8 *8 


d = 18 

d = 42 

''1 ''2 

''1 '^3 


.054 

*166 

-442 

.787 

.949 

.994 

.999 
.3 22 
.868 
.996 
.842 
.99 2 
.98 3 
.998 

.999 


.046 

.079 

.209 

.468 

.751 

*9 28 

.988 

*215 

.638 

.964 

.7 29 

.944 

.957 

.987 

.996 


.050 

.086 

.2^9 

.519 

.795 

.949 

.99 2 

.170 

.669 

.977 

.669 

♦96 2 

.947 

.991 

.999 


.058 

.170 

.48 2 

.838 

.973 

.997 

•999 

.336 

.891 

.998 

♦876 

^96 

.99 2 

.999 

.999 


.050 

•030 

.2^ 

.540 

.837 

.975 

.996 

• 244 

.717 

.984 

.787 

.974 

.980 

.995 

.999 


.056 

.085 

.264 

.591 

♦873 

.987 

.998 

♦186 

.754 

.990 

.746 

.979 

.972 

.996 

.999 
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Appendix a 


X 


JOLUTION OF THE ML EQUATIONS UNDER 
THE HYPOTHESIS 0^=6 2=6 

For a given set of ordered sample observations 

i -■ i i yr2+2 ^ "* ^ 

dn dlgoiritlim is prdsontdd to compute the ML estimates of 
location parameter 0 and scale parameter cr under the hypothesis 
In Section 2-6, for and for different 

valxies of rj| and the ML estimates 0 and o of 0 and o are 
given as follows s 


( a) 

Case % 

rj^=0,r2=0 



§ r: 


S = p/d*, 

(b) 

Case s 

r^ > 0, 

r2=° 



A 

0 = 


log 

(l+rj^/f) and a = P/d*, 

( c) 

Case ; 


> 0 


( i) 

if a 

< b, then S = 

: y -a log and a 

r2+i 

( ii) 

if a 

> b, then 0 = 

; X . and a is the solution 

r^+1 

(A,l) 

G3<P 

11 

<o 

1 + 


(d) 

Case : 

r^ > 0, 


■ 0 


A 

0 = 


log 




15 2 


and a is the solution of the equation 


(A. 2 ) exp (Q/a) = {1 y __ .,^1^ . 

-r j, 


P~d"a 


d*d+fQ~P 


where p 






. ^ X _ + r, V^+s^y 

i=r^+l “ ^ j=r2+l ^ ^ ” 2 “® 


2 2 


r^+1" 


a = Q/log(l+r 2 /f ) and b = (P-fQ)/d^. 

As remarked in Section 2.6^ equations (A.i) and (A , 2 ) 
are solved by Newton-Raphson method with (2P-fQ)/2d* as an 
initial value* 

The subroutine ESTMAT" computes the ML estimates THETA 
and SIGMA of 0 and a respectively for given values of n^^n^* irpy 
and S 2 f and vectors X and Y. 5'or this subprogram, the 
rc)quirc:sd accuracy AC, and the expected number of iterations NI 
are supplied from the main calling program* The failure 
indicator FI takes the value 1, if the iteration procedure does 
not converge in NI stc^s* 

LANGUAGE 


Fortran 10 
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STRUCTURE 


SUBROUTINE 


ESTMAT(n1,N2,R1,R2,S1^ 


s 2^ X, Y, N I> AC, THETA, S TGMA , F I ) 


pa -TTai net e rs 


N1 

Integer 

N2 

Integer 

R1 

Integer 

R2 

Integer 

SI 

Integer 

S2 

Integer 

X 

Real 

Y 

Real 

NI 

Integer 


AC 

Real 

THETA 

Real 

SIGMA 

Real 

FI 

Integer 


input 

input 

Input 

input ; 

input ; 

input j 

input vector of : 
length Nl-Rl-Sl 

input vector of t 
length N2-'R2~S2 

input t 

input : 

output 5 

output : 

output ‘ 


size of the first sample 

size of the second sample 

number of sireillest 
observations missing in the 
fii'st sample 

number of smallest 
observations missing in the 
second sample 

number of largest observations 
missing in the first sample 

number of largest observations 
missing in the second sample 

available ordered observations! 
in the first sanple 

available ordered observations 
in the second sample 

upper bound for the number of : 
iterations in -t^ich the 
iteration pjcpcess is expected i 
to converge 

desired accuracy j 

ML estimate of 0 I 

i 

ML estimate of 0 | 

failure indicator 


1 if the iteration does not 
converge 

0 otherwise- 
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C THIS SUBROUTINE COMPUTES THE ML ESTIMATES OF THETA 

C AND SICMA BASED ON TYPE II DOUBLY CENSORED 

C samples FRCM TWO EXPONENTIAL DISTRIBUTIONS 

C 

SUBROUTINE ESTMAT (Nl,N2,Rl, R2/ SI, S2,X,Y,NI,AC,THETA^ 

1 SIGMA, FI) 

INTEGER Rl,R2,Sl,S2yFI 

DOUBLE PRECISION P,D,P,Q,A, B,X(999 ), Y(999 ), Z(999 ) , S( 999 ), 
ITHETA, SIGMA, SUM, AC,U1,U2,U-^,U4, FN, FP 
IF(X(R1+1) .LE.Y(R2+1) ) GO TO 40 
C INTERCHANGING THE SAMPLES 

C 

Nll=Ni;NRl=Ri;NSl=Sl 
N1=N2;R1=R2/S1=S2 
N 2=N1 1 ; R 2-NRl / s 2=NS 1 
DO 10 I=R1+1/N1-S1 
10 Z(l)=Y(l) 

DO 20 I=R2+1^N2-S2 
20 Y(l)=X(l) 

DO 30 I=RH-1/N1-S1 
30 X(l)=Z(l) 

C CALCULATION OF F,D AND Q 

C 

40 P=N1+N2-*R1-R2 

D=P~S1-*S2 

Q=Y(R2+l)-X^^^l+^^ 
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C SUM OP ALL OBSERVATIONS 

C 

SUM=0 

DO 50 I=R1+1,N1-.S1 
50 SUM=SUM+X(I) 

DO 60 I=R2+1,N2-S2 
60 SUM=SUM+Y( I) 

C CALCULATION OP P 

C 

P=SUM+S1*X(N1~S1 )+S2*y(N2-S2)-PJ^X( Rl+1 ) 

c calculation op a and b 

c 

A=Q/DLOG( 1 .+R^P) 

B=:(P-P*Q)/D 

C SEPARATION OP DIPPERENT CASES 

C 

IF(R2.EQ.0) GO TO 75 
IP( ( R1 .EQ *0 ) .and , ( A .LE ,3) ) GO TO 3 0 
C ESTIMATING SIGM\ BY NEWTON-RiiPHSON METHOD 

C WITH INITI/^ VxiLUE S(l) 

C 

S(l) =( 2.-'tD) 

1=1 

65 Ul=l .+R2^frQ/( P-D»S{T5 ) 

U2=1 *+R1*Q/(D*S( I) +Bi- Q-P) 

U3=R2i^QifrD/( (P~D*S( I) )«*2.) 

U4=R1-mQ*D/( (D*S( I)+5H‘Q~P) ^ 2 *^ 
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C THE FUNCTION OP SIGNA IS DENOTED BY " PN" 

C AND ITS DERIVATIVE IS BY"PP" 

C 

FN=DEXP( Q/S( I) ) -Ul/u 2 

FP=^DEXP(Q/S( I) )^Q/{ S( I)**2*)-(U3*U24U1*U4)/(U2 ^2-) 
S(l+l)=S(l)~FN/PP 

IP(DABS(S(I+1)~S(I)) .LE.AC) GO TO 70 
IF(I-GE,NI) GO TO 85 
I=I+I 
GO TO 65 

c FIN/Uj V.iLUE OF SIGMA AND THETA 

C 

70 SIGM\=S(l+l) 

THETA=X( Rl+1 )~SIGM/^DLOG( 1 ,+Rl*Q/ ( D*SIGMA+P*Q— P) ) 

GO TO 90 
75 SIGMA==P/D 

THETA=X( R1+1)-SIGMA5-DL0G( 1 .+R1/F) 

GO TO 90 

90 SIGMA=(P-B^Q)/D 

THETA=Y( R2+1 )-SIGMf'i*DLOG( 1 »+R2/'^) 

GO TO 90 

c ASSIGNING THE VALUE TO THE FAILURE INDICATOR 

c 

85 

9 0 RETURN 


END 
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C 

CALLING PROGRAM 
* * * -;{■ -» * * 

C 

C SOLUTION OF THE ML EQUATIONS 

INTEGER Rl, R2,S1, S2,PI 

DOUBLE PRECISION X( 100 )^Y( 100 ), AC, THETA, SIGm 
C INPUT VALUES, THE DATA X(l) AND Y(j) ARE SII€JLATED 

C VALUES FROM E(2, l) DISTRIBUTION ARRANGED IN AN ASCENDING 

C ORDER OF MAGNITUDE 

C 

D/ VTA N1 , N 2, R1 , R2, S 1 , S 2, NI, AC, ( X( I) , 1=4, 14 ) , ( Y( J ) , J=3, 12)/ 

1 17,13, 3, 2, 3,1, 20, .00001, 2-16140, 2.21918, 2*30073, 2-84808, 
1 2-91879, 3.14132, ^ *21995, 3 .34723, 3 .35 224, 3 *45513, 3 2251, 

1 2.15 214, 2*18 279, 2*29610, 2*30496, 2*41608, 2 *53415, 2 *57010, 
1 2 .9 5 275, 2*968 00,4 *32659/ 

CALL ESTM/VT( N1 , N 2, R1 , R2, SI, S 2, X, Y, NI, AC, THETA, S IGMA, F I ) 
IF(FI-l) 10,30,10 
10 PRINT 20, THETA, SIGMA 

20 F0RMAT(10X,'THE EST MiTE OF THETA IS' ,F9 .5/lOX, ^THE 
1 ESTIM/vTE OF SIGMA IS',F9-5) 

# 

GO TO 50 

30 PRINT 40, NI 

40 P0RMAT(5X,^THE ITER/ITION DOES NOT CONVERGE IN', 34, 'STEPS' ) 

50 STOP 


END 



Appendix b 


evaluation of Q^(x1s)anD L^(x!s) 

The procedures for evaluating the functions Q^( x i s ) and 
L^(xls), where d is a positive integer are presented* Function 
subprograms QD(D^X, s) and LD(D,XyS) are given for calculating 
these functions. 

Prom Section 3.1, we have 

Q^(xis) =/ dy/(d~l)l , x > 0. 

X 

Note that, this integral is convergent only for (l+s) > 0, that 
is for s > -1 . By substituting z = y(l+s), we get 

Q.(xls) = ; e"^ z^"^ dVf . 

x( l+s ) 

(B.l) -I e-x(l+s) , s > -1. 

j=0 

Similarly, for finite positive values of x, x^e have 

La( XI s) = ; 1"®’ y'^"^ ay/( a-i ) i . 

o , 

This converges for all s, but has to be treated s^arately for 
s « 1 and s / 1 cases- Thus 

x^/di / s = 1 

(B.2) L^(xis) = ^ 

. j=o 
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The FUNCTION QD(D, X/S) computes the suiWiation given 
in equation (B.i) for D = 1,2/.*./ X > 0 and S > -i- For 

= 1/2, X > 0 and -«> < S < oo^ FUNCTION LD(D,X^S) computes 

the value of L^(xls) given in equation (B»2)* 

language 

Fortran 10 

C FUNCTION SUBPROGRA-M FOR CO^IPUTING QD(XlS) 

C 

DOUBLE PRECISION FUNCTION QD(D,X,s) 

DOUBLE PRECISION SUM,X,S,X1^FAC 

INTEGER D 

SUM=0 .0 

DO 10 J=0^D-1 

XI =1 •0;l=j 

IP(J.BQ.0) GO TO 90 

xi=x*(i#+s) 

IF(J#BQ.l) GO TO 90 
DO 80 K=2,^ 
pac=dploat(k) 

0 0 Xl=Xl^«.X-K-( 1 *+2 )/FAC 

90 SUM=SUM+X1 

10 CONTItTUE 

QD=DIXP(-X*{ 1 (1 .tS )*.D) 

return 


END 



C FUNCTION SUBPROGRAM FOR CCMPUTING LD(XiS) 

C 

DOUBLE PRECISION FUNCTION LD(D^X^S) 
DOUBLE PRECISION X^ S, SUM, XI, FAC, SI 
INTEGER D 

IF(S,EQ.l) GO TO 21 
SUM=0 .0 
DO 15 J=0,D-1 
XI =1 .o;i=j 

IF(J*BQ#0) GO TO 90 

xi=»«-(i.~s) 

IF(J.EQ.I) GO TO 90 
DO 80 K=2/I 
fac=dfloat(k) 

8 0 XI =Xl-;c- X*( 1 •-S )/FAC 

90 SUM=SUM+X1 

15 CONTIIWE 

LD=( 1 •-DEXP(~X*(1*~S))*SUM)/((1*-S)**D) 

GO TO 23 

21 Sl=l*> 

DO 25 J=l/D 
PAC=DFL0AT{J) 

25 Sl:=Sl'}frV^^‘^ 

LD=S1 

23 return 

END 



